Changeset 837 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Nov 10, 2012, 4:55:04 PM (12 years ago)

Author:

aslmd

Message:

LMDZ.GENERIC. Corrected problems with allocated arrays in start2archive and newstart. Applied a workaround to make those work without tracers (-cpp NOTRAC -- perhaps there is a better solution). Checked that everything works in debug mode.

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• datareadnc.F (modified) (2 diffs)
• for_lmdz5/README.aymeric (modified) (1 diff)
• surfini.F (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/datareadnc.F

@@ -232 +232 @@
      .         iim, jjp1, rlonv, rlatu, zmea,zstd,zsig,zgam,zthe)
 
-          CALL dump2d(iip1,jjp1,zmea,'zmea')
-          CALL dump2d(iip1,jjp1,zstd,'zstd')
-          CALL dump2d(iip1,jjp1,zsig,'zsig')
-          CALL dump2d(iip1,jjp1,zgam,'zgam')
-          CALL dump2d(iip1,jjp1,zthe,'zthe')
+          !CALL dump2d(iip1,jjp1,zmea,'zmea')
+          !CALL dump2d(iip1,jjp1,zstd,'zstd')
+          !CALL dump2d(iip1,jjp1,zsig,'zsig')
+          !CALL dump2d(iip1,jjp1,zgam,'zgam')
+          !CALL dump2d(iip1,jjp1,zthe,'zthe')
 
       endif
@@ -295 +295 @@
             phisinit(i)=1000.*phisinit(i)
       ENDDO
-      CALL dump2d(iimp1,jjp1,phisinit,'Altitude en m')
+      !CALL dump2d(iimp1,jjp1,phisinit,'Altitude en m')
       CALL multscal(iimp1*jjp1,phisinit,g,phisinit)
 

trunk/LMDZ.GENERIC/libf/phystd/for_lmdz5/README.aymeric

@@ -40 +40 @@
 
 
+/home/spiga/MODELES/LMDZ.GENERIC/makegcm -d 64x48x20 -t 0 -debug -cpp NOTRAC start2archive > logcomp
+
+/home/spiga/MODELES/LMDZ.GENERIC/makegcm -d 64x48x20 -t 0 newstart
+
+/home/spiga/MODELES/LMDZ.GENERIC/makegcm -d 128x96x20 -t 0 -debug -cpp NOTRAC newstart
+
+/home/spiga/MODELES/LMDZ.GENERIC/makegcm -d 64x48x20 -t 0 -s 1 -b 3x2 gcm
+
+
 pour initialisation idealisee, qqch comme iniaqua dans phyaqua ???

trunk/LMDZ.GENERIC/libf/phystd/surfini.F

@@ -57 +57 @@
         DO ig=1,ngrid
           IF (qsurf(ig,igcm_co2_ice) .GT. 0.) THEN
+             !!!! no good in parallel
              IF(ig.GT.ngrid/2+1) THEN
                 icap=2