Changeset 1246 for trunk/LMDZ.MARS/libf/phymars/dimradmars_mod.F90

Timestamp:

May 11, 2014, 2:37:58 PM (11 years ago)

Author:

aslmd

Message:

LMDZ.MARS. Made number of scatterers a free dimension not in need to be prescribe at compiling time. Instead it must be set in callphys.def. See README for further information about this commit.

File:

• trunk/LMDZ.MARS/libf/phymars/dimradmars_mod.F90 (modified) (3 diffs)

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

@@ -16 +16 @@
 ! Number of kind of tracer radiative properties
 ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-! naerkind is set in scatterers.h (built when compiling with makegcm -s #)
-
-#include"scatterers.h"
-! NB: May have to change value of nsizemax below when changing scatterers
+! naerkind is set by reading callphys.def
+! -- see conf_phys
+! -- value of nsizemax below is comfortably high
+!    but beware in case you add a lot of scatterers
+  INTEGER, SAVE :: naerkind
+
+  ! AS: previously in aerkind.h
+  character*20, SAVE, ALLOCATABLE :: name_iaer(:)  ! name of the scatterers
+  integer iaer_dust_conrath ! Typical dust profiles using a
+                            ! Conrath type analytical equation
+  integer iaer_dust_doubleq ! Dust profile is given by the
+                            ! mass mixing ratio of the two-
+                            ! moment scheme method (doubleq)
+  integer iaer_dust_submicron ! Dust profile is given by a
+                              ! submicron population of dust
+                              ! particles
+  integer iaer_h2o_ice ! Water ice particles
+
+  ! AS: was in aeropacity
+  INTEGER,SAVE,ALLOCATABLE :: iaerdust(:)
+
+  ! AS: was in suaer
+  CHARACTER(LEN=30), SAVE, ALLOCATABLE :: file_id(:,:)
 
 ! Reference wavelengths used to compute reference optical depth (m)
 ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-
-  real,save :: longrefir(naerkind),longrefvis(naerkind)
+  REAL,SAVE,ALLOCATABLE :: longrefir(:),longrefvis(:)
   
 ! Definition of spectral intervals at thermal infrared wavelengths (LW)
@@ -86 +104 @@
   REAL,SAVE,ALLOCATABLE :: nueffdust(:,:) ! Dust effective variance
 
+!! ------------------------------------------------------
+!! AS: what was previously in yomaer
+!   Shortwave
+!   ~~~~~~~~~
+! 
+! tauvis: dust optical depth at reference wavelength  ("longrefvis" set
+! in dimradmars_mod : typically longrefvis = 0.67E-6 m, as measured by Viking )
+
+! For the "naerkind" kind of aerosol radiative properties : 
+! QVISsQREF  :  Qext / Qext("longrefvis")   <--- For both solar bands
+! omegavis   :  sinle scattering albedo     <--- For both solar bands
+! gvis       :  assymetry factor            <--- For both solar bands
+! 
+!   Longwave
+!   ~~~~~~~~
+! 
+! For the "naerkind" kind of aerosol radiative properties : 
+! QIRsQREF :  Qext / Qext("longrefvis")     <--- For the nir bandes IR
+! omegaIR  :  mean single scattering albedo <--- For the nir bandes IR
+! gIR      :  mean assymetry factor         <--- For the nir bandes IR
+! 
+  real,save :: tauvis
+  real,save,allocatable :: QVISsQREF(:,:,:)
+  real,save,allocatable :: omegavis(:,:,:)
+  real,save,allocatable :: gvis(:,:,:)
+  real,save,allocatable :: QIRsQREF(:,:,:)
+  real,save,allocatable :: omegaIR(:,:,:)
+  real,save,allocatable :: gIR(:,:,:)
+! Actual number of grain size classes in each domain for a
+!   given aerosol:
+  integer,save,allocatable :: nsize(:,:)
+! Particle size axis (depend on the kind of aerosol and the
+!   radiation domain)
+  real,save,allocatable :: radiustab(:,:,:)
+! Extinction coefficient at reference wavelengths;
+!   These wavelengths are defined in dimradmars_mod, and called
+!   longrefvis and longrefir.
+  real,save,allocatable :: QREFvis(:,:)
+  real,save,allocatable :: QREFir(:,:)
+  real,save,allocatable :: omegaREFvis(:,:)
+  real,save,allocatable :: omegaREFir(:,:)
+!! ------------------------------------------------------
+
 contains
-  
+
   subroutine ini_dimradmars_mod(ngrid,nlayer)
   
@@ -107 +168 @@
    allocate(fluxrad(ngrid))
    allocate(tauscaling(ngrid))
+   allocate(nueffdust(ngrid,nlayer))
+
+  end subroutine ini_dimradmars_mod
+ 
+  subroutine ini_scatterers(ngrid,nlayer)
+
+  implicit none
+
+  integer,intent(in) :: ngrid ! number of atmospheric columns
+  integer,intent(in) :: nlayer ! number of atmospheric layers
+
+   allocate(name_iaer(naerkind))
    allocate(aerosol(ngrid,nlayer,naerkind))
-   allocate(nueffdust(ngrid,nlayer))
-
-  end subroutine ini_dimradmars_mod
-  
+   allocate(longrefir(naerkind))
+   allocate(longrefvis(naerkind))
+   allocate(iaerdust(naerkind))
+   allocate(file_id(naerkind,2))
+
+   allocate(QVISsQREF(nsun,naerkind,nsizemax))
+   allocate(omegavis(nsun,naerkind,nsizemax))
+   allocate(gvis(nsun,naerkind,nsizemax))
+   allocate(QIRsQREF(nir,naerkind,nsizemax))
+   allocate(omegaIR(nir,naerkind,nsizemax))
+   allocate(gIR(nir,naerkind,nsizemax))
+   allocate(nsize(naerkind,2))
+   allocate(radiustab(naerkind,2,nsizemax))
+   allocate(QREFvis(naerkind,nsizemax))
+   allocate(QREFir(naerkind,nsizemax))
+   allocate(omegaREFvis(naerkind,nsizemax))
+   allocate(omegaREFir(naerkind,nsizemax))
+
+  end subroutine ini_scatterers
+
 end module dimradmars_mod