Changeset 4092 for trunk/LMDZ.VENUS/libf/phyvenus/compo_hedin_mod2.F90

Timestamp:

Feb 28, 2026, 6:36:58 PM (6 weeks ago)

Author:

emillour

Message:

Venus PCM:
Code cleanup: remove "hedin.h" and make variables included there be module
variables of compo_hedin_mod2.F90. Turn "nirdat.h" into module nirdata.F90
and incroporate initialization routine nir_leedat.F in it.
EM

File:

• trunk/LMDZ.VENUS/libf/phyvenus/compo_hedin_mod2.F90 (modified) (2 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/compo_hedin_mod2.F90

@@ -2 +2 @@
   !!05/2013 Laura Salmi
   !!03/2014 revision Gabriella Gilli
-  !!Calcul des vmr pour CO2, CO, O, N et N2 en s'appuyant sur les tables donnees dans l'article Hedin (1983)
+  !!Compute vmr for CO2, CO, O, N and N2 based on tables from Hedin (1983) paper
   
   implicit none      
+
+  integer, parameter :: zsize=23  ! for z from 100 to 210 km
+  integer, parameter :: musize=10 ! for sza ranging from -1 to 1
+  
+  real, save, protected :: co2_hedin(musize,zsize)
+  real, save, protected :: co_hedin(musize,zsize)
+  real, save, protected :: n2_hedin(musize,zsize)
+  real, save, protected :: o_hedin(musize,zsize)
+  real, save, protected :: n_hedin(musize,zsize)
+  
+  real, save, protected :: pres_hedin(musize,zsize)
+  real, save, protected :: mu_hedin(musize)
+  
   contains
 
         subroutine compo_hedin83_init2
-        use YOMCST_mod
+        use YOMCST_mod, only: RKBOL
         implicit none
-!#include "YOMCST.h"
-#include "hedin.h"
 
         REAL :: alpha
@@ -133 +144 @@
 !!Interpolation des profils de Hedin (1983) sur la grille du modele
 
-        use dimphy
+        use dimphy, only: klon, klev
         implicit none
-        include 'hedin.h'
+
         REAL, intent(in) :: pression(klon,klev), mu0(klon)
         integer :: i,k,iz,jmu,jz,z_ok, mu_ok