Changeset 2597

Timestamp:

Dec 17, 2021, 10:28:22 AM (3 years ago)

Author:

emillour

Message:

Mars GCM:
Move "zls" (solar longitude) from being a local variable of physiq
to module comsaison_h. This is for convenience (saving zls and making
it available to the "plasma dynamical core").
Also cleaned up and commented comsaison_h in the process.
JYC+EM

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/comsaison_h.F90 (modified) (1 diff)
• libf/phymars/physiq_mod.F (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -3549 +3549 @@
 - missing rho at launch altitude
 - changed inputs, because R and Cp are needed to compute rho and BV
+
+== 17/12/2021 == JYC+EM
+Move "zls" (solar longitude) from being a local variable of physiq
+to module comsaison_h. This is for convenience (saving zls and making
+it available to the "plasma dynamical core").
+Also cleaned up and commented comsaison_h in the process.
+

trunk/LMDZ.MARS/libf/phymars/comsaison_h.F90

@@ -3 +3 @@
 implicit none
 
-  logical,save :: callsais
-  integer,save :: isaison
-  real,save :: dist_sol
-  real,save :: declin
-  real,save,allocatable :: mu0(:)
-  real,save,allocatable :: fract(:)
+  real,save :: dist_sol ! Sun-Mars distance (AU)
+  real,save :: declin ! declination (rad)
+  real,save :: zls ! solar longitude (rad)
+  real,save,allocatable :: mu0(:) ! cosine of solar zenith angle
+  real,save,allocatable :: fract(:) ! day fraction over the time interval (0: night, 1: day)
   real,save,allocatable :: local_time(:) ! local solar time as fraction of day (0,1)
 
-!$OMP THREADPRIVATE(callsais,isaison,dist_sol,declin,mu0,fract,local_time)
+!$OMP THREADPRIVATE(dist_sol,declin,zls,mu0,fract,local_time)
 
 contains

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -54 +54 @@
      &                     capcal, fluxgrd, qsurf,
      &                     hmons,summit,base,watercap,watercaptag
-      use comsaison_h, only: dist_sol, declin, mu0, fract, local_time
+      use comsaison_h, only: dist_sol, declin, zls, 
+     &                       mu0, fract, local_time
       use slope_mod, only: theta_sl, psi_sl
       use conc_mod, only: rnew, cpnew, mmean
@@ -316 +317 @@
       REAL tau(ngrid,naerkind)     ! Column dust optical depth at each point
                                    ! AS: TBD: this one should be in a module !
-      REAL zls                       !  solar longitude (rad)
       REAL zday                      ! date (time since Ls=0, in martian days)
       REAL zzlay(ngrid,nlayer)     ! altitude at the middle of the layers