Changeset 23082

Timestamp:

08/08/18 11:33:50 (7 years ago)

Author:

dmenemen

Message:

cpl_issm now checks for issm mask (9999's) and sets

R_shelfIce=0; this change does not modify Archive4002.arch

cpl_issm now updates R_shelfIce at every coupling time step instead

of just during statTime; this change modifies Archive4002.arch

Location:

issm/trunk-jpl/test

Files:

• Archives/Archive4002.arch (modified) ( previous)
• MITgcm/code/cpl_issm.F (modified) (4 diffs)

issm/trunk-jpl/test/MITgcm/code/cpl_issm.F

@@ -51 +51 @@
       _R8 xfer_array(Nx,Ny)
       _R8 local(1:sNx,1:sNy,nSx,nSy)
+      CHARACTER*(MAX_LEN_MBUF) suff
 
 C Initialization steps I1, I2, and I3:
@@ -182 +183 @@
 C Recurring step C4 except during Termination:
 C  C4. ISSM sends IceModelTime=(N-1)*CouplingTime base to MITgcm
-CDM   IF( myTime .NE. endtime .AND.
-      IF( myTime .EQ. startTime .AND.
+      IF( myTime .NE. endtime .AND.
      &     MOD(myTime,CouplingTime) .LT. deltaT/2. ) THEN
-         CALL WRITE_FLD_XY_RS( 'R_shelfIce_1',' ',R_shelfIce,-1,myThid)
+         WRITE(suff,'(I10.10)') myIter
+         CALL WRITE_FLD_XY_RS( 'R_shelfIce1_',suff,R_shelfIce,-1,myThid)
          IF( myProcId .EQ. 0 ) THEN
             _BEGIN_MASTER( myThid )         
@@ -199 +200 @@
                DO j=1,sNy
                   DO i=1,sNx
-                     IF( ABS(local(i,j,bi,bj)+
-     &                    0.89190188 ).GT.1E-5 ) THEN
+                     IF( local(i,j,bi,bj).LT.9998 ) THEN
                         R_shelfIce(i,j,bi,bj) = local(i,j,bi,bj)
+                     ELSE
+                        R_shelfIce(i,j,bi,bj) = 0. _d 0
                      ENDIF
                   ENDDO
@@ -209 +211 @@
 C- fill in the overlap (+ BARRIER):
          _EXCH_XY_RS( R_shelfIce, myThid )
-         CALL WRITE_FLD_XY_RS( 'R_shelfIce_2',' ',R_shelfIce,-1,myThid)
+         CALL WRITE_FLD_XY_RS( 'R_shelfIce2_',suff,R_shelfIce,-1,myThid)
       ENDIF
 C End recurring step C4.