Skip to content
Snippets Groups Projects
default_desfmn.f90 39.8 KiB
Newer Older
Gaelle TANGUY's avatar
Gaelle TANGUY committed
  LCOLD  = .TRUE.
  LNUCL  = .TRUE.
Gaelle TANGUY's avatar
Gaelle TANGUY committed
  LSNOW  = .TRUE.
  LHAIL  = .FALSE.
  CPRISTINE_ICE_LIMA = 'PLAT'
  CHEVRIMED_ICE_LIMA = 'GRAU'
  XFACTNUC_DEP = 1.0  
  XFACTNUC_CON = 1.0
  NMOD_IFN = 1
  NIND_SPECIE = 1
  LMEYERS = .FALSE.
  LIFN_HOM = .TRUE.
  CIFN_SPECIES = 'PHILLIPS'
  CINT_MIXING = 'DM2'
  XIFN_CONC(:) = 100.
Gaelle TANGUY's avatar
Gaelle TANGUY committed
  NMOD_IMM = 0
  NPHILLIPS=8
ENDIF
!
!-------------------------------------------------------------------------------
!
!*      20.   SET DEFAULT VALUES FOR MODD_CH_MNHC_n
!             -------------------------------------
!
LUSECHEM            = .FALSE.
LUSECHAQ            = .FALSE.
LUSECHIC            = .FALSE.
LCH_INIT_FIELD      = .FALSE.
LCH_CONV_SCAV       = .FALSE.
LCH_CONV_LINOX      = .FALSE.
LCH_PH              = .FALSE.
LCH_RET_ICE         = .FALSE.
XCH_PHINIT          = 5.2
XRTMIN_AQ           = 5.e-8
CCHEM_INPUT_FILE    = 'EXSEG1.nam'
CCH_TDISCRETIZATION = 'SPLIT'
NCH_SUBSTEPS        = 1
LCH_TUV_ONLINE      = .FALSE.
CCH_TUV_LOOKUP      = 'PHOTO.TUV39'
CCH_TUV_CLOUDS      = 'NONE'
XCH_TUV_ALBNEW      = -1.
XCH_TUV_DOBNEW      = -1.
XCH_TUV_TUPDATE     = 600.
CCH_VEC_METHOD      = 'MAX'
NCH_VEC_LENGTH      = 50
XCH_TS1D_TSTEP      = 600.
CCH_TS1D_COMMENT    = 'no comment'
CCH_TS1D_FILENAME   = 'IO1D'
!
!-------------------------------------------------------------------------------
!
!*      21.   SET DEFAULT VALUES FOR MODD_SERIES AND MODD_SERIE_n
!             ---------------------------------------------------
!
IF (KMI == 1) THEN
  LSERIES      = .FALSE.
  LMASKLANDSEA = .FALSE.
  LWMINMAX     = .FALSE.
  LSURF        = .FALSE.
ENDIF
!
NIBOXL = 1 !+ JPHEXT
NIBOXH = 1 !+ 2*JPHEXT
NJBOXL = 1 !+ JPHEXT
NJBOXH = 1 !+ 2*JPHEXT
NKCLS  = 1 !+ JPVEXT
NKLOW  = 1 !+ JPVEXT
NKMID  = 1 !+ JPVEXT
NKUP   = 1 !+ JPVEXT
NKCLA  = 1 !+ JPVEXT
NBJSLICE = 1
NJSLICEL(:) = 1 !+ JPHEXT
NJSLICEH(:) = 1 !+ 2*JPHEXT
NFREQSERIES  = INT(XSEGLEN /(100.*XTSTEP) )
NFREQSERIES  = MAX(NFREQSERIES,1)
!
!-------------------------------------------------------------------------------
!
!*      22.   SET DEFAULT VALUES FOR MODD_TURB_CLOUD
!             --------------------------------------
!
IF (KMI == 1) THEN
  NMODEL_CLOUD = NUNDEF
  CTURBLEN_CLOUD = 'DELT'
  XCOEF_AMPL_SAT = 5.
  XCEI_MIN = 0.001E-06
  XCEI_MAX = 0.01E-06
ENDIF
!-------------------------------------------------------------------------------
!
!*      22.   SET DEFAULT VALUES FOR MODD_MEAN_FIELD
!             --------------------------------------
!
IF (KMI == 1) THEN
  LMEAN_FIELD = .FALSE.
ENDIF
!
!-------------------------------------------------------------------------------
!
!*      22.   SET DEFAULT VALUES FOR MODD_AEROSOL
!             -----------------------------------
IF (KMI == 1) THEN ! other values are defined in modd_ch_aerosol
!
! aerosol lognormal parameterization

LVARSIGI  = .FALSE.   ! switch to active pronostic dispersion for I mode
LVARSIGJ  = .FALSE.   ! switch to active pronostic dispersion for J mode
LHETEROSO4 = .FALSE.  ! switch to active sulfates heteronegeous
                      ! production
LSEDIMAERO = .FALSE.  ! switch to active aerosol sedimentation
LAERINIT   = .FALSE.  ! switch to initialize aerosol in arome
CMINERAL      = "NONE"   ! mineral equilibrium scheme
CORGANIC      = "NONE"   ! mineral equilibrium scheme
CNUCLEATION   = "NONE" ! sulfates nucleation scheme
LDEPOS_AER(:) = .FALSE.

ENDIF

!*      23.   SET DEFAULT VALUES FOR MODD_DUST and MODD_SALT
!             ----------------------------------------------
!
IF (KMI == 1) THEN ! other values initialized in modd_dust
  LDUST      = .FALSE.
  NMODE_DST  = 3
  LVARSIG    = .FALSE.
  LSEDIMDUST = .FALSE.
  LDEPOS_DST(:) = .FALSE.

  LSALT      = .FALSE.
  LVARSIG_SLT= .FALSE.
  LSEDIMSALT = .FALSE.
  LDEPOS_SLT(:)     = .FALSE.
ENDIF
!
!-------------------------------------------------------------------------------
!
!
!*      24.   SET DEFAULT VALUES FOR MODD_PASPOL           
!             ----------------------------------
!
! other values initialized in modd_paspol
!
IF (KMI == 1) THEN 
  LPASPOL    = .FALSE.
  NRELEASE   = 0
  CPPINIT(:) ='1PT'
  XPPLAT(:)  = 0.
  XPPLON (:) = 0.
  XPPMASS(:) = 0.
  XPPBOT(:)  = 0.
  XPPTOP(:)  = 0.
  CPPT1(:)   = "20010921090000"
  CPPT2(:)   = "20010921090000"
  CPPT3(:)   = "20010921091500"
  CPPT4(:)   = "20010921091500"
ENDIF
!
!-------------------------------------------------------------------------------
!
!
!*      25.   SET DEFAULT VALUES FOR MODD_CONDSAMP         
!             ----------------------------------
!
! other values initialized in modd_condsamp
!
IF (KMI == 1) THEN 
  LCONDSAMP  = .FALSE.
  NCONDSAMP  = 3
  XRADIO(:)  = 900.
  XSCAL(:)   = 1.
  XHEIGHT_BASE = 100.
  XDEPTH_BASE  = 100.
  XHEIGHT_TOP  = 100.
  XDEPTH_TOP   = 100.
ENDIF  
!-------------------------------------------------------------------------------
!
!
!*      26.   SET DEFAULT VALUES FOR MODD_LATZ_EDFLX      
!             ----------------------------------
!
IF (KMI == 1) THEN 
  LUV_FLX=.FALSE.
  XUV_FLX1=3.E+14
  XUV_FLX2=0.
  LTH_FLX=.FALSE.
  XTH_FLX=0.75
ENDIF  
#ifdef MNH_FOREFIRE
!-------------------------------------------------------------------------------
!
!*      27.   SET DEFAULT VALUES FOR MODD_FOREFIRE         
!             ----------------------------------
!
! other values initialized in modd_forefire
!
IF (KMI == 1) THEN 
  LFOREFIRE = .FALSE.
  LFFCHEM = .FALSE.
  COUPLINGRES = 100.
  NFFSCALARS = 0
ENDIF  
#endif                 
!-------------------------------------------------------------------------------
!
!*      28.   SET DEFAULT VALUES FOR MODD_BLOWSNOW AND  MODD_BLOWSNOW_n       
!             ----------------------------------------
! 
IF (KMI == 1) THEN
   LBLOWSNOW  = .FALSE.
   XALPHA_SNOW  = 3.
   XRSNOW       = 4.
   CSNOWSEDIM  = 'TABC'
END IF
LSNOWSUBL = .FALSE.
!
!
!-------------------------------------------------------------------------------
!
!*      29.   SET DEFAULT VALUES FOR MODD_VISC           
!             ----------------------------------
!
! other values initialized in modd_VISC
!
IF (KMI == 1) THEN
  LVISC    = .FALSE.
  LVISC_UVW    = .FALSE.
  LVISC_TH    = .FALSE.
  LVISC_SV    = .FALSE.
  LVISC_R    = .FALSE.
  XPRANDTL = 0.
ENDIF
!
!-------------------------------------------------------------------------------
!
!
!*      30.   SET DEFAULT VALUES FOR MODD_DRAG           
!             ----------------------------------
!
! other values initialized in modd_DRAG
!
IF (KMI == 1) THEN
  LDRAG    = .FALSE.
  LMOUNT   = .FALSE.
  NSTART = 1
  XHSTART = 0.
ENDIF
!
!-------------------------------------------------------------------------------
!
!*      31.   SET DEFAULT VALUES FOR MODD_IBM_PARAMn         
!             --------------------------------------
!                                                           
  LIBM           = .FALSE.
  LIBM_TROUBLE   = .FALSE.
  CIBM_ADV       = 'NOTHIN'
  XIBM_EPSI      = 1.E-9
  XIBM_IEPS      = 1.E+9
  NIBM_ITR       = 8  
  XIBM_RUG       = 0.01    ! (m^1.s^-0)
  XIBM_VISC      = 1.56e-5 ! (m^2.s^-1)
  XIBM_CNU       = 0.06    ! (m^0.s^-0)  

  NIBM_LAYER_P = 2
  NIBM_LAYER_Q = 2
  NIBM_LAYER_R = 2
  NIBM_LAYER_S = 2
  NIBM_LAYER_T = 2
  NIBM_LAYER_E = 2
  NIBM_LAYER_V = 2

  XIBM_RADIUS_P = 2.
  XIBM_RADIUS_Q = 2.
  XIBM_RADIUS_R = 2.
  XIBM_RADIUS_S = 2.
  XIBM_RADIUS_T = 2.
  XIBM_RADIUS_E = 2.
  XIBM_RADIUS_V = 2.

  XIBM_POWERS_P = 1.
  XIBM_POWERS_Q = 1.
  XIBM_POWERS_R = 1.
  XIBM_POWERS_S = 1.
  XIBM_POWERS_T = 1.
  XIBM_POWERS_E = 1.
  XIBM_POWERS_V = 1.

  CIBM_MODE_INTE3_P = 'LAI'
  CIBM_MODE_INTE3_Q = 'LAI'
  CIBM_MODE_INTE3_R = 'LAI'
  CIBM_MODE_INTE3_S = 'LAI'
  CIBM_MODE_INTE3_T = 'LAI'
  CIBM_MODE_INTE3_E = 'LAI'
  CIBM_MODE_INTE3_V = 'LAI'

  CIBM_MODE_INTE1_P = 'CL2'
  CIBM_MODE_INTE1_Q = 'CL2'
  CIBM_MODE_INTE1_R = 'CL2'
  CIBM_MODE_INTE1_S = 'CL2'
  CIBM_MODE_INTE1_T = 'CL2'
  CIBM_MODE_INTE1_E = 'CL2'
  CIBM_MODE_INTE1NV = 'CL2'
  CIBM_MODE_INTE1TV = 'CL2'
  CIBM_MODE_INTE1CV = 'CL2'

  CIBM_MODE_BOUND_P = 'SYM'
  CIBM_MODE_BOUND_Q = 'SYM'
  CIBM_MODE_BOUND_R = 'SYM'
  CIBM_MODE_BOUND_S = 'SYM'
  CIBM_MODE_BOUND_T = 'SYM'
  CIBM_MODE_BOUND_E = 'SYM'
  CIBM_MODE_BOUNT_V = 'ASY'
  CIBM_MODE_BOUNN_V = 'ASY'
  CIBM_MODE_BOUNC_V = 'ASY'

  XIBM_FORC_BOUND_P = 0.
  XIBM_FORC_BOUND_Q = 0.
  XIBM_FORC_BOUND_R = 0.
  XIBM_FORC_BOUND_S = 0. 
  XIBM_FORC_BOUND_T = 0.
  XIBM_FORC_BOUND_E = 0.
  XIBM_FORC_BOUNN_V = 0.
  XIBM_FORC_BOUNT_V = 0.
  XIBM_FORC_BOUNC_V = 0.

  CIBM_TYPE_BOUND_P = 'NEU'
  CIBM_TYPE_BOUND_Q = 'NEU'
  CIBM_TYPE_BOUND_R = 'NEU'
  CIBM_TYPE_BOUND_S = 'NEU'
  CIBM_TYPE_BOUND_T = 'NEU'
  CIBM_TYPE_BOUND_E = 'NEU'
  CIBM_TYPE_BOUNT_V = 'DIR'
  CIBM_TYPE_BOUNN_V = 'DIR'
  CIBM_TYPE_BOUNC_V = 'DIR'

  CIBM_FORC_BOUND_P = 'CST'
  CIBM_FORC_BOUND_Q = 'CST'
  CIBM_FORC_BOUND_R = 'CST'
  CIBM_FORC_BOUND_S = 'CST'
  CIBM_FORC_BOUND_T = 'CST'
  CIBM_FORC_BOUND_E = 'CST'
  CIBM_FORC_BOUNN_V = 'CST'
  CIBM_FORC_BOUNT_V = 'CST'
  CIBM_FORC_BOUNC_V = 'CST'
  CIBM_FORC_BOUNR_V = 'CST'

!
!-------------------------------------------------------------------------------
!
!*      32.   SET DEFAULT VALUES FOR MODD_RECYCL_PARAMn         
!             --------------------------------------
!
  LRECYCL  = .FALSE.
  LRECYCLN = .FALSE.
  LRECYCLW = .FALSE.
  LRECYCLE = .FALSE.
  LRECYCLS = .FALSE.
  XDRECYCLN  = 0.
  XTMOY      = 0.
  XTMOYCOUNT = 0.
  XNUMBELT   = 28.
  XRCOEFF    = 0.2
  XTBVTOP    = 500.
  XTBVBOT    = 300.
END SUBROUTINE DEFAULT_DESFM_n