Changeset 1824 for trunk/LMDZ.COMMON/libf/dyn3d_common

Timestamp:

Nov 21, 2017, 4:03:44 PM (7 years ago)

Author:

emillour

Message:

Common dynamics:

• enable possiblity to store multiple time steps in the restart.nc file (flag "ecritstart" gives the frequency, in dynamical steps).
• fixed dynredem_mod.F90 to correctly write multiple time steps.
• fixed computation of JH_cur in the mars case where "hour_ini" contains the initial time of day read from the start.nc file
• minor fix in dynetat0.F90

RY

Location:

trunk/LMDZ.COMMON/libf/dyn3d_common

Files:

• dynetat0.F90 (modified) (1 diff)
• dynredem_mod.F90 (modified) (2 diffs)

trunk/LMDZ.COMMON/libf/dyn3d_common/dynetat0.F90

@@ -307 +307 @@
     s1='value of '//TRIM(str1)//' ='
     s2=' read in starting file differs from parametrized '//TRIM(str2)//' ='
-    WRITE(msg,'(10x,a,i4,2x,a,i4)'),s1,n1,s2,n2
+    WRITE(msg,'(10x,a,i4,2x,a,i4)')TRIM(s1),n1,TRIM(s2),n2
     CALL ABORT_gcm(TRIM(modname),TRIM(msg),1)
   END IF

trunk/LMDZ.COMMON/libf/dyn3d_common/dynredem_mod.F90

@@ -34 +34 @@
   INTEGER :: start(4), count(4)
 !===============================================================================
-  start(:)=[1,1,1,1]; count(:)=[iip1,jjp1,ll,nb]
+  IF (ll.eq.1) THEN
+    start(:)=[1,1,nb,1]
+    count(:)=[iip1,jjp1,1,1]
+  ELSE
+    start(:)=[1,1,1,nb]
+    count(:)=[iip1,jjp1,ll,1]
+  ENDIF
   CALL err(NF90_INQ_VARID(ncid,id,nvarid),"inq",id)
   CALL err(NF90_PUT_VAR(ncid,nvarid,var,start,count),"put",id)
@@ -60 +66 @@
   INTEGER :: start(4), count(4)
 !===============================================================================
-  start(:)=[1,1,1,1]; count(:)=[iip1,jjm,ll,nb]
+  IF (ll.eq.1) THEN
+    start(:)=[1,1,nb,1]
+    count(:)=[iip1,jjm,1,1]
+  ELSE
+    start(:)=[1,1,1,nb]
+    count(:)=[iip1,jjm,ll,1]
+  ENDIF
   CALL err(NF90_INQ_VARID(ncid,id,nvarid),"inq",id)
   CALL err(NF90_PUT_VAR(ncid,nvarid,var,start,count),"put",id)