Changeset 234 for trunk/LMDZ.MARS/libf/phymars/meso_inc

Timestamp:

Jul 19, 2011, 5:25:58 PM (14 years ago)

Author:

aslmd

Message:

LMDZ.MARS
+ MESOSCALE

• An important change to merge physiq.F and inifis.F for GCM and mesoscale
• This is mostly transparent to GCM users and developers (use of MESOSCALE precompiler flags)
• Makes it easy (and mostly automatic!) for changes in GCM physics to be impacted in mesoscale physics
• A few minor changes have followed in the GCM (slope scheme, one-point diagnostic).

• Compilation + run is OK on both sides (GCM and mesoscale).
• On the mesoscale side, call_meso_physiq?.inc and call_meso_inifis?.inc have been changed accordingly.

Here is the excerpt from README file:

19/07/2011 == AS

• Finished converging meso_physiq.F and meso_inifis.F towards physiq.F and inifis.F --> see previous point 15/07/2011 --> meso_ routines no longer exist (everything is in meso_inc and transparent to GCM users) --> GCM routines include mesoscale parts within MESOSCALE precompiler commands --> MESOSCALE parts are as hidden as possible not to mess up with GCM users/developers
• Cleaned inelegant or useless #ifdef [or] #ifndef MESOSCALE in physiq and inifis so that a minimum amount of such precompiler commands is now reached [mainly related to I/O]
• Added the SF08 slope insolation model in the general physics parameterizations. Added a callslope keyword in inifis.F and callkeys.h --> This keyword is False by default / True if you use -DMESOSCALE
• Removed the obsolete call to Viking Lander 1 diagnostic Replaced it with a diagnostic for opacity at the domain center [valid for GCM and mesoscale]

Location:

trunk/LMDZ.MARS/libf/phymars/meso_inc

Files:

• meso_inc_inifisini.F (modified) (2 diffs)
• meso_inc_inifisinvar.F (modified) (1 diff)
• meso_inc_inifisvar.F (modified) (2 diffs)
• meso_inc_slope.F (deleted)
• meso_inc_var.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_inifisini.F

@@ -8 +8 @@
 c      
       ! in 'comcstfi.h'
-      daysec=wdaysec  
       omeg=womeg                
       mugaz=wmugaz  
@@ -89 +88 @@
 c It must be set now, because it is used afterwards
 c*****************************************************
-        dtphys=wdt*float(wappel_phys)
+        dtphys=wdt*float(ptimestep)
         print*,'Physical timestep (s) ',dtphys

trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_inifisinvar.F

@@ +1 @@
+     $           ,nq,wdt
      $           ,womeg,wmugaz
      $           ,wyear_day,wperiheli,waphelie,wperi_day,wobliquit 

trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_inifisvar.F

@@ -1 @@
-      INTEGER nq,wday_ini
+      INTEGER nq
+      REAL wdt
 
-      REAL womeg,wmugaz,wdaysec
+      REAL womeg,wmugaz
       REAL wyear_day,wperiheli,waphelie,wperi_day,wobliquit
       REAL wz0,wemin_turb,wlmixmin
@@ -12 +13 @@
       REAL wtheta(ngrid),wpsi(ngrid)
       REAL wvolcapa
-      REAL wdt
-      INTEGER wappel_phys

trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_var.F

@@ -10 +10 @@
       REAL output_tab2d(ngridmx,n2d)
       REAL output_tab3d(ngridmx,nlayer,n3d)
-      REAL sl_ls, sl_lct, sl_lat, sl_tau, sl_alb, sl_the, sl_psi
-      REAL sl_fl0, sl_flu
-      REAL sl_ra, sl_di0
-      REAL sky
       REAL hfx(ngridmx)    !! pour LES avec isfflx!=0
       REAL ust(ngridmx)    !! pour LES avec isfflx!=0