Changeset 232 for codes/icosagcm/branches

Timestamp:

07/17/14 13:35:34 (10 years ago)

Author:

ymipsl

Message:

Work around for init etat0 state => missing field in the startfi are initialize to constant value.

YM

Location:

codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd

Files:

• phyetat0_academic.F90 (added)
• physiq.F90 (modified) (2 diffs)

codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/physiq.F90

@@ -559 +559 @@
 !        read startfi (initial parameters)
 !        ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-         call phyetat0(ngrid,nlayer,"startfi.nc",0,0,nsoilmx,nq,   &
+!ym         call phyetat0(ngrid,nlayer,"startfi.nc",0,0,nsoilmx,nq,   &
+!ym               day_ini,time_phys,tsurf,tsoil,emis,q2,qsurf,   &
+!ym               cloudfrac,totcloudfrac,hice,  &
+!ym               rnat,pctsrf_sic,tslab, tsea_ice,sea_ice)
+
+         call phyetat0_academic(ngrid,nlayer,"startfi.nc",0,0,nsoilmx,nq,   &
                day_ini,time_phys,tsurf,tsoil,emis,q2,qsurf,   &
                cloudfrac,totcloudfrac,hice,  &
@@ -569 +574 @@
            write(*,*) "dynamics day: ",pday
            write(*,*) "physics day:  ",day_ini
-           stop
+           write(*,*) "taking dynamics day for physics:  ",day_ini
+           day_ini=pday
+!ym           stop
          endif