/[MITgcm]/MITgcm_contrib/dgoldberg/streamice/streamice_init_varia.F
ViewVC logotype

Diff of /MITgcm_contrib/dgoldberg/streamice/streamice_init_varia.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.12 by dgoldberg, Thu Mar 7 15:23:19 2013 UTC revision 1.14 by dgoldberg, Sat Jun 8 22:15:34 2013 UTC
# Line 65  C     ZERO OUT FLOATING POINT ARRAYS Line 65  C     ZERO OUT FLOATING POINT ARRAYS
65            C_basal_friction(i,j,bi,bj) = C_basal_fric_const            C_basal_friction(i,j,bi,bj) = C_basal_fric_const
66            B_glen(i,j,bi,bj) = B_glen_isothermal            B_glen(i,j,bi,bj) = B_glen_isothermal
67            H_streamice_prev(i,j,bi,bj) = 0. _d 0            H_streamice_prev(i,j,bi,bj) = 0. _d 0
68    #ifdef ALLOW_STREAMICE_2DTRACER
69    #ifdef STREAMICE_TRACER_AB
70             GAD_trac_2d (i,j,bi,bj) = 0. _d 0
71    !         GAD_trac_2dNm1 (i,j,bi,bj) = 0. _d 0
72    #endif
73    #endif
74  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
75            ru_old_si(i,j,bi,bj) = 0. _d 0            ru_old_si(i,j,bi,bj) = 0. _d 0
76            rv_old_si(i,j,bi,bj) = 0. _d 0            rv_old_si(i,j,bi,bj) = 0. _d 0
77            zu_old_si(i,j,bi,bj) = 0. _d 0            zu_old_si(i,j,bi,bj) = 0. _d 0
78            zv_old_si(i,j,bi,bj) = 0. _d 0            zv_old_si(i,j,bi,bj) = 0. _d 0
79            h_after_uflux_SI(i,j,bi,bj) = 0. _d 0  !          h_after_uflux_SI(i,j,bi,bj) = 0. _d 0
80  #endif  #endif
81  #ifdef USE_ALT_RLOW  #ifdef USE_ALT_RLOW
82            R_low_si(i,j,bi,bj) = 0. _d 0            R_low_si(i,j,bi,bj) = 0. _d 0
# Line 395  C The 0 is the "iteration" argument. The Line 401  C The 0 is the "iteration" argument. The
401    
402  ! finish initialize basal trac  ! finish initialize basal trac
403    
404    #ifdef ALLOW_STREAMICE_2DTRACER
405    
406          IF ( STREAMICETRAC2DINITFILE .NE. ' ' ) THEN
407             _BARRIER
408    
409            CALL READ_FLD_XY_RL( STREAMICETRAC2dInitFile, ' ',
410         &      trac2d, 0, myThid )
411    
412           ELSE
413            WRITE(msgBuf,'(A)') 'TRAC2dInit - NO FILE SPECIFIED'
414            CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
415         &                    SQUEEZE_RIGHT , 1)
416            DO bj = myByLo(myThid), myByHi(myThid)
417             DO bi = myBxLo(myThid), myBxHi(myThid)
418              DO j=1,sNy
419               DO i=1,sNx
420                trac2d(i,j,bi,bj) = 0.0
421               ENDDO
422              ENDDO
423             ENDDO
424            ENDDO
425    
426           ENDIF
427    
428          _EXCH_XY_RL (trac2d, myThid)
429    
430    
431    #endif /*STREAMICE_ALLOW_2DTRACER*/
432    
433        CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )        CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )
434    
435        _EXCH_XY_RL(H_streamice, myThid )        _EXCH_XY_RL(H_streamice, myThid )
# Line 430  C The 0 is the "iteration" argument. The Line 465  C The 0 is the "iteration" argument. The
465       & B_glen, 0, myThid )       & B_glen, 0, myThid )
466        CALL WRITE_FLD_XY_RL ( "H_streamice", "init",        CALL WRITE_FLD_XY_RL ( "H_streamice", "init",
467       & H_streamIce, 0, myThid )       & H_streamIce, 0, myThid )
468    #ifdef ALLOW_STREAMICE_2DTRACER
469          CALL WRITE_FLD_XY_RL ( "2DTracer", "init",
470         & trac2d, 0, myThid )
471    #endif
472        CALL WRITE_FLD_XY_RL ( "area_shelf_streamice", "init",        CALL WRITE_FLD_XY_RL ( "area_shelf_streamice", "init",
473       & area_shelf_streamice, 0, myThid )       & area_shelf_streamice, 0, myThid )
474        CALL WRITE_FLD_XY_RL ( "STREAMICE_hmask", "init",        CALL WRITE_FLD_XY_RL ( "STREAMICE_hmask", "init",
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22