Changeset 1589 for trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd

Timestamp:

Aug 30, 2016, 9:04:28 AM (8 years ago)

Author:

emillour

Message:

Generic GCM:
Fix newstart which was broken due to recent updates of the physics/dynamics interface.
EM

File:

• trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/newstart.F (modified) (6 diffs)

trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/newstart.F

@@ -204 +204 @@
       
 ! get value of nsoilmx and allocate corresponding arrays
-      nsoilmx=18 ! default value
+      ! a default value of nsoilmx is set in comsoil_h
       call getin_p("nsoilmx",nsoilmx)
       
@@ -333 +333 @@
         enddo
         
-        ! Lmodif set to 0 to disable modifications possibility in phyeta0                           
         write(*,*) 'Reading file START'
         fichnom = 'start.nc'
@@ -339 +338 @@
      .       ps,phis,time)
 
+        CALL iniconst
+        CALL inigeom
+
+! Initialize the physics
+         CALL iniphysiq(iim,jjm,llm,
+     &                  (jjm-1)*iim+2,comm_lmdz,
+     &                  daysec,day_ini,dtphys,
+     &                  rlatu,rlatv,rlonu,rlonv,
+     &                  aire,cu,cv,rad,g,r,cpp,
+     &                  1)
+
+        ! Lmodif set to 0 to disable modifications possibility in phyeta0                           
+        Lmodif=0
         write(*,*) 'Reading file STARTFI'
         fichnom = 'startfi.nc'
@@ -441 +453 @@
 ! Load parameters from run.def file
       CALL defrun_new( 99, .TRUE. )
+! Initialize dynamics geometry and co. (which, when using start.nc may 
+! have changed because of modifications (tabi, preff, pa) above)
       CALL iniconst 
       CALL inigeom
@@ -446 +460 @@
       idum=0
 
-! Initialize the physics
+! Initialize the physics for start_archive only
+      IF (choix_1.eq.0) THEN
          CALL iniphysiq(iim,jjm,llm,
      &                  (jjm-1)*iim+2,comm_lmdz,
@@ -453 +468 @@
      &                  aire,cu,cv,rad,g,r,cpp,
      &                  1)
+      ENDIF
 
 c=======================================================================