Changeset 1710 for trunk/ICOSA_LMDZ/src

Timestamp:

May 18, 2017, 2:45:16 PM (8 years ago)

Author:

aslmd

Message:

corrected starting day bug. either no startfi.nc then starting day is read in *.def; or there is a startfi and starting day is read in startfi. AS+EM

File:

• trunk/ICOSA_LMDZ/src/phystd/interface_icosa_lmdz.f90 (modified) (4 diffs)

trunk/ICOSA_LMDZ/src/phystd/interface_icosa_lmdz.f90

@@ -197 +197 @@
   USE inifis_mod, ONLY : inifis 
 !  USE phyaqua_mod, ONLY : iniaqua 
-   
+ 
+  USE netcdf  
   
   IMPLICIT NONE
@@ -239 +240 @@
   INTEGER :: iflag_phys    
   INTEGER :: nq
+
+  !! to get starting date
+  !! --------------------
+  logical :: startphy_file
+  ! NetCDF stuff
+  integer :: status ! NetCDF return code
+  integer :: ncid ! NetCDF file ID
+  integer :: varid ! NetCDF variable ID
+  real :: tab_cntrl(100)
 
     CALL init_distrib_icosa_lmdz
@@ -339 +349 @@
 !    CALL inifis(klon_omp,nlayer,nqtot,pdayref,punjours,nday,ptimestep, &
 !            latitude,longitude,cell_area,prad,pg,pr,pcpp)
-    start_day=0
-    CALL getin('start_day',start_day)
+!$OMP MASTER
+
+    startphy_file=.true.
+    CALL getin('startphy_file',startphy_file)
+
+    IF (startphy_file) THEN
+
+      status=nf90_open('startfi.nc',NF90_NOWRITE,ncid)
+      if (status.ne.nf90_noerr) then
+        write(*,*)"Failed to open startfi.nc"
+        write(*,*)trim(nf90_strerror(status))
+        stop
+      endif
+
+      status=nf90_inq_varid(ncid,"controle",varid)
+      if (status.ne.nf90_noerr) then
+        write(*,*)"Failed to find controle variable"
+        write(*,*)trim(nf90_strerror(status))
+        stop
+      endif
+
+      status=nf90_get_var(ncid,varid,tab_cntrl)
+      start_day=tab_cntrl(3)
+      print*,"start_day YORGL",start_day
+      status=nf90_close(ncid)
+
+    ELSE
+
+      start_day=0
+      CALL getin('start_day',start_day)
+
+    ENDIF
+
     day_length=86400
     CALL getin('day_length',day_length)
+    year_length = 31557600 ! 365.25 * 86400
+    CALL getin('year_length',year_length)
     ndays=nint(itaumax*(dt/day_length))! number of days to run
     physics_timestep=dt*itau_physics
+!$OMP END MASTER
+!$OMP BARRIER
     CALL inifis(klon_omp,llm,nqtot,start_day,day_length,ndays,physics_timestep, &
             latfi,lonfi,airefi,radius,g,kappa*cpp,cpp)
@@ -373 +418 @@
     ! Initialize "calendar"
 !$OMP MASTER
-    day_length=86400
-    CALL getin('day_length',day_length)
-    year_length = 31557600 ! 365.25 * 86400
-    CALL getin('year_length',year_length)
-    start_day = 0
-    CALL getin('start_day',start_day)
     ! initialize pday and ptime
     pday = start_day