Changeset 1535 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Apr 8, 2016, 8:46:54 AM (9 years ago)

Author:

emillour

Message:

Mars:

• Fix for the 1D model initializations.

EM

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F (modified) (5 diffs)

trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F

@@ -12 +12 @@
       use slope_mod, only: theta_sl, psi_sl
       use phyredem, only: physdem0,physdem1
-      use comgeomphy, only: initcomgeomphy
+      use comgeomphy, only: initcomgeomphy, airephy, cuphy, cvphy,
+     &                      rlond, rlatd
       use planete_h, only: year_day, periheli, aphelie, peri_day,
      &                     obliquit, emin_turb, lmixmin
       use comcstfi_h, only: pi, rad, omeg, g, mugaz, rcp, r, cpp
-      use time_phylmdz_mod, only: daysec, dtphys, day_step
+      use time_phylmdz_mod, only: daysec, dtphys, day_step,
+     &                            ecritphy, iphysiq
       use dimradmars_mod, only: tauscaling,tauvis
       USE comvert_mod, ONLY: ap,bp,aps,bps,pa,preff,sig
       USE logic_mod, ONLY: hybrid
+      use regular_lonlat_mod, only: init_regular_lonlat
       USE phys_state_var_init_mod, ONLY: phys_state_var_init
       IMPLICIT NONE
@@ -117 +120 @@
 ! initialize "serial/parallel" related stuff
 !      CALL init_phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
-      CALL init_phys_lmdz(iim,jjm+1,llm,1,(/(jjm-1)*iim+2/))
+      CALL init_phys_lmdz(1,1,llm,1,(/1/))
       call initcomgeomphy
 
@@ -419 +422 @@
       write(*,*) " day_step = ",day_step
 
+      iphysiq=1 ! in 1D model physics are called every time step
+      ecritphy=day_step ! default value for ecritphy, output every time step
+
       ndt=10 ! default value for ndt
       PRINT *,'Number of sols to run ?'
@@ -477 +483 @@
 !  and allocates some arrays 
 !Mars possible matter with dtphys in input and include!!!
+! Initializations below should mimick what is done in iniphysiq for 3D GCM
+      call init_regular_lonlat(1,1,longitude,latitude,
+     &                   (/0.,0./),(/0.,0./))
       call phys_state_var_init(1,llm,nq,
      .          day0,time,daysec,dtphys,rad,g,r,cpp)
@@ -526 +535 @@
 c   mesh surface (not a very usefull quantity in 1D)
       area(1)=1.E+0
+      airephy(1)=area(1)
 
 c   surface geopotential is not used (or useful) since in 1D