/[MITgcm]/MITgcm_contrib/shelfice_remeshing/CLEAN/code/shelfice_thermodynamics.F
ViewVC logotype

Diff of /MITgcm_contrib/shelfice_remeshing/CLEAN/code/shelfice_thermodynamics.F

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

revision 1.3 by dgoldberg, Mon Jan 18 10:36:41 2016 UTC revision 1.4 by dgoldberg, Tue Jan 19 16:46:55 2016 UTC
# Line 731  C     endif (not) useISOMIPTD Line 731  C     endif (not) useISOMIPTD
731  C--  Calculate new loading anomaly (in case the ice-shelf mass was updated)  C--  Calculate new loading anomaly (in case the ice-shelf mass was updated)
732  #ifndef ALLOW_AUTODIFF  #ifndef ALLOW_AUTODIFF
733    
 c       DO bj = myByLo(myThid), myByHi(myThid)  
 c        DO bi = myBxLo(myThid), myBxHi(myThid)  
 c         DO j = 1-OLy, sNy+OLy  
 c          DO i = 1-OLx, sNx+OLx  
 c           IF (R_shelfice(i,j,bi,bj) .LT. 0) THEN    
 c             ETACOUNT=ETACOUNT+1  
 c             GMSL=GMSL+etaN(i,j,bi,bj)  
 c           ENDIF  
 c          ENDDO  
 c         ENDDO  
 c        ENDDO  
 c       ENDDO  
 c           GMSL=GMSL/ETACOUNT  
 c  
 c       DO bj = myByLo(myThid), myByHi(myThid)  
 c        DO bi = myBxLo(myThid), myBxHi(myThid)  
 c         DO j = 1-OLy, sNy+OLy  
 c          DO i = 1-OLx, sNx+OLx  
 c         IF (shelficeMass(i,j,bi,bj).GT. 1030*(990+GMSL)) THEN  
 c             shelficeMass(i,j,bi,bj) = 1030*(990+GMSL)  
 c         ENDIF  
 c  
 c           shelficeLoadAnomaly(i,j,bi,bj) = gravity  
 c     &      *( shelficeMass(i,j,bi,bj) + rhoConst*Ro_surf(i,j,bi,bj) )  
 c  
 c          ENDDO  
 c         ENDDO  
 c        ENDDO  
 c       ENDDO  
 c      ENDIF  
   
734    
735         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
736          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
# Line 789  c      ENDIF Line 758  c      ENDIF
758           shelficeLoadAnomaly(i,j,bi,bj) = gravity           shelficeLoadAnomaly(i,j,bi,bj) = gravity
759       &      *( mass_eff + rhoConst*Ro_surf(i,j,bi,bj) )       &      *( mass_eff + rhoConst*Ro_surf(i,j,bi,bj) )
760    
761               EFFMASS(i,j,bi,bj)=mass_eff
762    
763  #endif /* ALLOW_STREAMICE */  #endif /* ALLOW_STREAMICE */
764    
# Line 828  C     Friction velocity Line 797  C     Friction velocity
797        ENDIF        ENDIF
798          CALL DIAGNOSTICS_FILL(R_shelfice,'SHI_Rshelfice',          CALL DIAGNOSTICS_FILL(R_shelfice,'SHI_Rshelfice',
799       &      0,1,0,1,1,myThid)       &      0,1,0,1,1,myThid)
800            CALL DIAGNOSTICS_FILL(EFFMASS,'SHI_MassEff',
801         &      0,1,0,1,1,myThid)
802  #endif /* ALLOW_DIAGNOSTICS */          
803    #endif
804    
805    
806  #endif /* ALLOW_SHELFICE */  #endif /* ALLOW_SHELFICE */
807        RETURN        RETURN
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22