Changeset 3572 for trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90

Timestamp:

Jan 13, 2025, 2:52:33 PM (3 weeks ago)

Author:

debatzbr

Message:

Remove generic_aerosols and generic_condensation, along with their related variables (useless). RENAME THE VARIABLE AEROHAZE TO OPTICHAZE.

File:

• trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90 (modified) (4 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90

@@ -776 +776 @@
      if (is_master)write(*,*)trim(rname)//&
      "Radiatively active haze ?"
-     aerohaze=.false. ! default value
-     call getin_p("aerohaze",aerohaze)
-     if (is_master)write(*,*)trim(rname)//&
-     "aerohaze = ",aerohaze
+     optichaze=.false. ! default value
+     call getin_p("optichaze",optichaze)
+     if (is_master)write(*,*)trim(rname)//&
+     "optichaze = ",optichaze
 
      if (is_master)write(*,*)trim(rname)//&
@@ -1198 +1198 @@
 
 !=================================
-! TWOLAY scheme and NH3 cloudare left for retrocompatibility only,
-! You should now use N-LAYER scheme (see below).
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": Radiatively active two-layer aerosols?"
-     aeroback2lay=.false.     ! default value
-     call getin_p("aeroback2lay",aeroback2lay)
-     if (is_master) write(*,*)trim(rname)//": aeroback2lay = ",aeroback2lay
-
-     if (aeroback2lay.and.is_master) then
-       print*,'Warning : The TWOLAY AEROSOL scheme is deprecated and buggy...'
-       print*,'You should use the generic n-layer scheme (see aeronlay).'
-     endif
-
-     if (is_master) write(*,*)trim(rname)//": Radiatively active ammonia cloud?"
-     aeronh3=.false.     ! default value
-     call getin_p("aeronh3",aeronh3)
-     if (is_master) write(*,*)trim(rname)//": aeronh3 = ",aeronh3
-
-     if (aeronh3.and.is_master) then
-       print*,'Warning : You are using specific NH3 cloud scheme ...'
-       print*,'You should use the generic n-layer scheme (see aeronlay).'
-     endif
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": TWOLAY AEROSOL: total optical depth "//&
-       "in the tropospheric layer (visible)"
-     obs_tau_col_tropo=8.D0
-     call getin_p("obs_tau_col_tropo",obs_tau_col_tropo)
-     if (is_master) write(*,*)trim(rname)//&
-       ": obs_tau_col_tropo = ",obs_tau_col_tropo
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": TWOLAY AEROSOL: total optical depth "//&
-       "in the stratospheric layer (visible)"
-     obs_tau_col_strato=0.08D0
-     call getin_p("obs_tau_col_strato",obs_tau_col_strato)
-     if (is_master) write(*,*)trim(rname)//&
-       ": obs_tau_col_strato = ",obs_tau_col_strato
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": TWOLAY AEROSOL: optprop_back2lay_vis?"
-     optprop_back2lay_vis = 'optprop_saturn_vis_n20.dat'
-     call getin_p("optprop_back2lay_vis",optprop_back2lay_vis)
-     if (is_master) write(*,*)trim(rname)//&
-       ": optprop_back2lay_vis = ",trim(optprop_back2lay_vis)
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": TWOLAY AEROSOL: optprop_back2lay_ir?"
-     optprop_back2lay_ir = 'optprop_saturn_ir_n20.dat'
-     call getin_p("optprop_back2lay_ir",optprop_back2lay_ir)
-     if (is_master) write(*,*)trim(rname)//&
-       ": optprop_back2lay_ir = ",trim(optprop_back2lay_ir)
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": TWOLAY AEROSOL: pres_bottom_tropo? in pa"
-     pres_bottom_tropo=66000.0
-     call getin_p("pres_bottom_tropo",pres_bottom_tropo)
-     if (is_master) write(*,*)trim(rname)//&
-       ": pres_bottom_tropo = ",pres_bottom_tropo
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": TWOLAY AEROSOL: pres_top_tropo? in pa"
-     pres_top_tropo=18000.0
-     call getin_p("pres_top_tropo",pres_top_tropo)
-     if (is_master) write(*,*)trim(rname)//&
-       ": pres_top_tropo = ",pres_top_tropo
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": TWOLAY AEROSOL: pres_bottom_strato? in pa"
-     pres_bottom_strato=2000.0
-     call getin_p("pres_bottom_strato",pres_bottom_strato)
-     if (is_master) write(*,*)trim(rname)//&
-       ": pres_bottom_strato = ",pres_bottom_strato
-
-     ! Sanity check
-     if (pres_bottom_strato .gt. pres_top_tropo) then
-       if(is_master) then
-       print*,'Error : TWOLAY AEROSOL, Please ensure that in callphys.def'
-       print*,'you have pres_top_tropo > pres_bottom_strato !'
-       endif
-       call abort_physic(rname," pres_top_tropo > pres_bottom_strato!",1)
-     endif
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": TWOLAY AEROSOL: pres_top_strato? in pa"
-     pres_top_strato=100.0
-     call getin_p("pres_top_strato",pres_top_strato)
-     if (is_master) write(*,*)trim(rname)//&
-       ": pres_top_strato = ",pres_top_strato
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": TWOLAY AEROSOL: particle size in the ", &
-       "tropospheric layer, in meters"
-     size_tropo=2.e-6
-     call getin_p("size_tropo",size_tropo)
-     if (is_master) write(*,*)trim(rname)//": size_tropo = ",size_tropo
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": TWOLAY AEROSOL: particle size in the ", &
-       "stratospheric layer, in meters"
-     size_strato=1.e-7
-     call getin_p("size_strato",size_strato)
-     if (is_master) write(*,*)trim(rname)//": size_strato = ",size_strato
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": NH3 (thin) cloud: total optical depth"
-     tau_nh3_cloud=7.D0
-     call getin_p("tau_nh3_cloud",tau_nh3_cloud)
-     if (is_master) write(*,*)trim(rname)//": tau_nh3_cloud = ",tau_nh3_cloud
-
-     if (is_master) write(*,*)trim(rname)//": NH3 (thin) cloud pressure level"
-     pres_nh3_cloud=7.D0
-     call getin_p("pres_nh3_cloud",pres_nh3_cloud)
-     if (is_master) write(*,*)trim(rname)//": pres_nh3_cloud = ",pres_nh3_cloud
-
-     if (is_master) write(*,*)trim(rname)//": NH3 (thin) cloud: particle sizes"
-     size_nh3_cloud=3.e-6
-     call getin_p("size_nh3_cloud",size_nh3_cloud)
-     if (is_master) write(*,*)trim(rname)//": size_nh3_cloud = ",size_nh3_cloud
-
-!=================================
-! Generic N-LAYER aerosol scheme
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": Radiatively active generic n-layer aerosols?"
-     aeronlay=.false.     ! default value
-     call getin_p("aeronlay",aeronlay)
-     if (is_master) write(*,*)trim(rname)//": aeronlay = ",aeronlay
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": Number of generic aerosols layers?"
-     nlayaero=1     ! default value
-     call getin_p("nlayaero",nlayaero)
-     ! Avoid to allocate arrays of size 0
-     if (aeronlay .and. nlayaero.lt.1) then
-       if (is_master) then
-       print*, " You are trying to set no generic aerosols..."
-       print*, " Set aeronlay=.false. instead ! I abort."
-       endif
-       call abort_physic(rname,"no generic aerosols...",1)
-     endif
-     if (.not. aeronlay) nlayaero=1
-     if (is_master) write(*,*)trim(rname)//": nlayaero = ",nlayaero
-
-     ! This is necessary, we just set the number of aerosol layers
-     IF(.NOT.ALLOCATED(aeronlay_tauref))      ALLOCATE(aeronlay_tauref(nlayaero))
-     IF(.NOT.ALLOCATED(aeronlay_lamref))      ALLOCATE(aeronlay_lamref(nlayaero))
-     IF(.NOT.ALLOCATED(aeronlay_choice))      ALLOCATE(aeronlay_choice(nlayaero))
-     IF(.NOT.ALLOCATED(aeronlay_pbot))        ALLOCATE(aeronlay_pbot(nlayaero))
-     IF(.NOT.ALLOCATED(aeronlay_ptop))        ALLOCATE(aeronlay_ptop(nlayaero))
-     IF(.NOT.ALLOCATED(aeronlay_sclhght))     ALLOCATE(aeronlay_sclhght(nlayaero))
-     IF(.NOT.ALLOCATED(aeronlay_size))        ALLOCATE(aeronlay_size(nlayaero))
-     IF(.NOT.ALLOCATED(aeronlay_nueff))       ALLOCATE(aeronlay_nueff(nlayaero))
-     IF(.NOT.ALLOCATED(optprop_aeronlay_ir))  ALLOCATE(optprop_aeronlay_ir(nlayaero))
-     IF(.NOT.ALLOCATED(optprop_aeronlay_vis)) ALLOCATE(optprop_aeronlay_vis(nlayaero))
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": Generic n-layer aerosols: Optical depth at reference wavelenght"
-     aeronlay_tauref=1.0E-1
-     call getin_p("aeronlay_tauref",aeronlay_tauref)
-     if (is_master) write(*,*)trim(rname)//&
-       ": aeronlay_tauref = ",aeronlay_tauref
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": Generic n-layer aerosols: Reference wavelenght for optical depths (m)"
-     aeronlay_lamref=0.6E-6
-     call getin_p("aeronlay_lamref",aeronlay_lamref)
-     if (is_master) write(*,*)trim(rname)//&
-       ": aeronlay_lamref = ",aeronlay_lamref
-
-     if (is_master) then
-       write(*,*)trim(rname)//&
-       ": Generic n-layer aerosols: Vertical profile choice : "
-       write(*,*)trim(rname)//&
-       "             (1) Tau btwn ptop and pbot follows atm. scale height"
-       write(*,*)trim(rname)//&
-       "         or  (2) Tau above pbot follows its own scale height"
-     endif
-     aeronlay_choice=1
-     call getin_p("aeronlay_choice",aeronlay_choice)
-     if (is_master) write(*,*)trim(rname)//&
-       ": aeronlay_choice = ",aeronlay_choice
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": Generic n-layer aerosols: bottom pressures (Pa)"
-     aeronlay_pbot=2000.0
-     call getin_p("aeronlay_pbot",aeronlay_pbot)
-     if (is_master) write(*,*)trim(rname)//": aeronlay_pbot = ",aeronlay_pbot
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": Generic n-layer aerosols: (if choice=1) Top pressures (Pa) "
-     aeronlay_ptop=300000.0
-     call getin_p("aeronlay_ptop",aeronlay_ptop)
-     if (is_master) write(*,*)trim(rname)//": aeronlay_ptop = ",aeronlay_ptop
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": Generic n-layer aerosols: (if choice=2) Scale height / atm. scale height"
-     aeronlay_sclhght=0.2
-     call getin_p("aeronlay_sclhght",aeronlay_sclhght)
-     if (is_master) write(*,*)trim(rname)//&
-       ": aeronlay_sclhght = ",aeronlay_sclhght
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": Generic n-layer aerosols: particles effective radii (m)"
-     aeronlay_size=1.e-6
-     call getin_p("aeronlay_size",aeronlay_size)
-     if (is_master) write(*,*)trim(rname)//": aeronlay_size = ",aeronlay_size
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": Generic n-layer aerosols: particles radii effective variance"
-     aeronlay_nueff=0.1
-     call getin_p("aeronlay_nueff",aeronlay_nueff)
-     if (is_master) write(*,*)trim(rname)//": aeronlay_nueff = ",aeronlay_nueff
-
-     if (is_master) write(*,*)trim(rname)//&
-       ": Generic n-layer aerosols: VIS optical properties file"
-     optprop_aeronlay_vis = 'optprop_saturn_vis_n20.dat'
-     call getin_p("optprop_aeronlay_vis",optprop_aeronlay_vis)
-     if (is_master) write(*,*)trim(rname)//&
-       ": optprop_aeronlay_vis = ",optprop_aeronlay_vis
-
-     if (is_master) write(*,*)trim(rname)//&
-     ": Generic n-layer aerosols: IR optical properties file"
-     optprop_aeronlay_ir = 'optprop_saturn_ir_n20.dat'
-     call getin_p("optprop_aeronlay_ir",optprop_aeronlay_ir)
-     if (is_master) write(*,*)trim(rname)//&
-       ": optprop_aeronlay_ir = ",optprop_aeronlay_ir
-
-
-!=================================
 
      if (is_master) write(*,*)trim(rname)//&
@@ -1460 +1229 @@
      call getin_p("sedimentation",sedimentation)
      if (is_master) write(*,*)trim(rname)//": sedimentation = ",sedimentation
-
-     if (is_master) write(*,*)trim(rname)//": Generic Condensation of tracers ?"
-     generic_condensation=.false. !default value
-     call getin_p("generic_condensation",generic_condensation)
-     if (is_master) write(*,*)trim(rname)//": generic_condensation = ",generic_condensation
 
      if (is_master) write(*,*)trim(rname)//": Spectral Dependant albedo ?"
@@ -1533 +1297 @@
          call abort_physic(rname, 'if paleo is set, fast must be true', 1)
      endif
-     if ((haze_proffix.or.haze_radproffix).and..not.aerohaze) then
-         call abort_physic(rname, 'for now, haze/rad proffix only works w aerohaze=T', 1)
+     if ((haze_proffix.or.haze_radproffix).and..not.optichaze) then
+         call abort_physic(rname, 'for now, haze/rad_proffix only works with optichaze=T', 1)
      endif
       if (carbox.and.condcosurf.and.no_n2frost) then