Changeset 1246 for trunk/LMDZ.MARS/libf/phymars/conf_phys.F

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/conf_phys.F (modified) (4 diffs)

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

@@ -1 @@
-      SUBROUTINE conf_phys(nq)
+      SUBROUTINE conf_phys(ngrid,nlayer,nq)
  
 !=======================================================================
@@ -38 +38 @@
       use surfdat_h, only: albedo_h2o_ice, inert_h2o_ice,
      &                     frost_albedo_threshold
-      use yomaer_h,only: tauvis
       use control_mod, only: ecritphy
       use planete_h
       USE comcstfi_h, only: daysec,dtphys
+      use dimradmars_mod, only: naerkind, name_iaer,
+     &                      ini_scatterers,tauvis
 
       IMPLICIT NONE
@@ -52 +53 @@
 !#include "surfdat.h"
 !#include "dimradmars.h"
-!#include "yomaer.h"
 #include "datafile.h"
 !#include "slope.h"
 #include "microphys.h"
 !#include "tracer.h"
-! naerkind is set in scatterers.h (built when compiling with makegcm -s #)
-#include"scatterers.h"
-
-      INTEGER,INTENT(IN) :: nq
-      INTEGER ig,ierr
+
+      INTEGER,INTENT(IN) :: ngrid,nlayer,nq
+      INTEGER ig,ierr,j
  
       CHARACTER ch1*12
@@ -536 +534 @@
 
 
+! SCATTERERS
+         write(*,*) "how many scatterers?"
+         naerkind=1 ! default value
+         call getin("naerkind",naerkind)
+         write(*,*)" naerkind = ",naerkind
+
+! Test of incompatibility
+c        Logical tests for radiatively active water-ice clouds:
+         IF ( (activice.AND.(.NOT.water)).OR.
+     &        (activice.AND.(naerkind.LT.2)) ) THEN
+           WRITE(*,*) 'If activice is TRUE, water has to be set'
+           WRITE(*,*) 'to TRUE, and "naerkind" must be at least'
+           WRITE(*,*) 'equal to 2.'
+           CALL ABORT
+         ELSE IF ( (.NOT.activice).AND.(naerkind.GT.1) ) THEN
+           WRITE(*,*) 'naerkind is greater than unity, but'
+           WRITE(*,*) 'activice has not been set to .true.'
+           WRITE(*,*) 'in callphys.def; this is not logical!'
+           CALL ABORT
+         ENDIF
+
+!------------------------------------------
+!------------------------------------------
+! once naerkind is known allocate arrays
+! -- we do it here and not in phys_var_init
+! -- because we need to know naerkind first
+         CALL ini_scatterers(ngrid,nlayer)
+!------------------------------------------
+!------------------------------------------
+
+
+c        Please name the different scatterers here ----------------
+         name_iaer(1) = "dust_conrath"   !! default choice is good old Conrath profile
+         IF (doubleq.AND.active) name_iaer(1) = "dust_doubleq" !! two-moment scheme
+         if (nq.gt.1) then
+          ! trick to avoid problems compiling with 1 tracer
+          ! and picky compilers who know name_iaer(2) is out of bounds
+          j=2
+         IF (water.AND.activice) name_iaer(j) = "h2o_ice"      !! radiatively-active clouds
+         IF (submicron.AND.active) name_iaer(j) = "dust_submicron" !! JBM experimental stuff
+         endif ! of if (nq.gt.1)
+c        ----------------------------------------------------------
+
 ! THERMOSPHERE