Changeset 2424

Timestamp:

Oct 22, 2020, 11:50:20 AM (4 years ago)

Author:

adelavois

Message:

Update of DYNAMICO / Mars interface with isotopes nqperes and nqfils

File:

• trunk/ICOSA_LMDZ/src/phymars/interface_icosa_lmdz.f90 (modified) (5 diffs)

trunk/ICOSA_LMDZ/src/phymars/interface_icosa_lmdz.f90

@@ -40 +40 @@
   
   CHARACTER(len=30),SAVE,ALLOCATABLE :: tname(:) ! tracer names
+  INTEGER,SAVE :: dyn_nqperes
+  INTEGER,SAVE,ALLOCATABLE :: dyn_nqfils(:)
   REAL,SAVE :: pday ! number of ellapsed sols since Ls=0
   REAL,SAVE :: ptime ! "universal time" as fraction of sol (e.g. 0.5 for noon)
@@ -224 +226 @@
   integer :: varid ! NetCDF variable ID
   real :: tab_cntrl(100)
+  real :: time0 !! Variable in startfi.nc to determine day_ini and hour_ini
 
     CALL init_distrib_icosa_lmdz
@@ -323 +326 @@
 
       status=nf90_get_var(ncid,varid,tab_cntrl)
-      day_ini=tab_cntrl(4)
+      day_ini=tab_cntrl(3)
+!      print*,"initialize_physics_omp: day_ini",day_ini
+      hour_ini=tab_cntrl(4)
+!      print*,"initialize_physics_omp: hour_ini",hour_ini
+
+      status=nf90_inq_varid(ncid,"Time",varid)
+      if (status.ne.nf90_noerr) then
+        write(*,*)"Failed to find Time variable"
+        write(*,*)trim(nf90_strerror(status))
+        stop
+      endif
+
+      status=nf90_get_var(ncid,varid,time0)
+      time0=int(time0) !AD: test fpr 
+      day_ini = day_ini + int(time0)
+      time0   = time0 - int(time0) 
+      hour_ini = hour_ini +time0
       print*,"initialize_physics_omp: day_ini",day_ini
-      hour_ini=tab_cntrl(29)
       print*,"initialize_physics_omp: hour_ini",hour_ini
       status=nf90_close(ncid)
@@ -332 +350 @@
       ENDIF ! of !IF (is_mpi_root)
       
+! pday is from dynamics, day_ini is then calculated in phyetat0 for physics
+! iteration is nb of dyn timesteps (as an integer in start.nc) 
+! pday=nint((iteration*dt)/day_length)
+! ptime= (iteration*dt)/day_length  -  pday
+
       CALL bcast(day_ini)
       CALL bcast(hour_ini)
@@ -354 +377 @@
 !!!            latfi,lonfi,airefi,radius,g,kappa*cpp,cpp)
 
+!!!    Temporary solution in order to run physics after HDO tracer implementation
+!!!    
+   dyn_nqperes   = nqtot
+   allocate(dyn_nqfils(dyn_nqperes))
+   dyn_nqfils(:) = 0
+!!!
    CALL phys_state_var_init(klon_omp,llm,nqtot,tname, &
                        day_ini,hour_ini,day_length,physics_timestep, &
-                       radius,g,kappa*cpp,cpp)
+                       radius,g,kappa*cpp,cpp, &
+                       dyn_nqperes,dyn_nqfils)
 
    CALL ini_fillgeom(klon_omp,latfi,lonfi,airefi)