/[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.13 by dgoldberg, Thu May 23 22:12:33 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 319  C The 0 is the "iteration" argument. The Line 325  C The 0 is the "iteration" argument. The
325    
326  ! finish initialize glen constant  ! finish initialize glen constant
327    
 ! initialize BDOT  
   
       IF ( STREAMICEBdotConfig.EQ.'FILE' ) THEN  
   
        IF ( STREAMICEBdotFile .NE. ' ' ) THEN  
         _BARRIER  
   
         CALL READ_FLD_XY_RL( STREAMICEBdotFile, ' ',  
      &      BDOT_STREAMICE, 0, myThid )  
   
        ELSE  
         WRITE(msgBuf,'(A)') 'INIT GLEN - FILENAME MISSING'  
         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
      &                    SQUEEZE_RIGHT , 1)  
        ENDIF  
   
 !      ELSE IF (STREAMICEBdotConfig.EQ.'UNIFORM' ) THEN  
 !  
 !        DO bj = myByLo(myThid), myByHi(myThid)  
 !         DO bi = myBxLo(myThid), myBxHi(myThid)  
 !          DO j=1,sNy  
 !           DO i=1,sNx  
 !            BDOT_streamice(i,j,bi,bj) = B_glen_isothermal  
 !           ENDDO  
 !          ENDDO  
 !         ENDDO  
 !        ENDDO  
   
       ELSEIF (STREAMICEBdotConfig.EQ.' ' ) THEN  
   
         DO bj = myByLo(myThid), myByHi(myThid)  
          DO bi = myBxLo(myThid), myBxHi(myThid)  
           DO j=1,sNy  
            DO i=1,sNx  
             BDOT_streamice(i,j,bi,bj) = 0. _d 0  
            ENDDO  
           ENDDO  
          ENDDO  
         ENDDO  
   
       ELSE  
   
        WRITE(msgBuf,'(A)') 'INIT MELT RATE - NOT IMPLENTED'  
        CALL PRINT_ERROR( msgBuf, myThid)  
        STOP 'ABNORMAL END: S/R STREAMICE_INIT_VAR'  
       ENDIF  
   
 ! finish initialize BDOT  
   
   
328  ! initialize basal traction  ! initialize basal traction
329    
330        IF ( STREAMICEbasalTracConfig.EQ.'FILE' ) THEN        IF ( STREAMICEbasalTracConfig.EQ.'FILE' ) THEN
# Line 445  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 480  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.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22