Changeset 3378

Timestamp:

Jun 13, 2024, 2:52:13 PM (6 months ago)

Author:

afalco

Message:

Pluto PCM: changed default params to match aerosol data. AF

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• inifis_mod.F90 (modified) (3 diffs)
• radinc_h.F90 (modified) (7 diffs)

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

@@ -629 +629 @@
     ! use Pluto.old routines
      if (is_master) write(*,*) trim(rname)//": use vdifc from old Pluto ?"
-     oldplutovdifc=.true. ! default value
+     oldplutovdifc=.false. ! default value
      call getin_p("oldplutovdifc",oldplutovdifc)
      if (is_master) write(*,*) trim(rname)//": oldplutovdifc = ",oldplutovdifc
@@ -635 +635 @@
      if (is_master) write(*,*) trim(rname)//&
        ": call pluto.old correlated-k radiative transfer ?"
-     oldplutocorrk=.true. ! default value
+     oldplutocorrk=.false. ! default value
      call getin_p("oldplutocorrk",oldplutocorrk)
      if (is_master) write(*,*) trim(rname)//": oldplutocorrk = ",oldplutocorrk
@@ -641 +641 @@
      if (is_master) write(*,*) trim(rname)//&
        ": call pluto.old sedimentation ?"
-     oldplutosedim=.true. ! default value
+     oldplutosedim=.false. ! default value
      call getin_p("oldplutosedim",oldplutosedim)
      if (is_master) write(*,*) trim(rname)//": oldplutosedim = ",oldplutosedim

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

@@ -11 +11 @@
 !======================================================================
 !
-!     RADINC.H 
+!     RADINC.H
 !
 !     Includes for the radiation code; RADIATION LAYERS, LEVELS,
 !     number of spectral intervals. . .
-! 
+!
 !======================================================================
 
@@ -27 +27 @@
 
 !     L_NLEVRAD corresponds to the surface - i.e., the GCM Level that
-!     is at the surface.  PLEV(L_NLEVRAD) = P(J,I)+PTROP, 
+!     is at the surface.  PLEV(L_NLEVRAD) = P(J,I)+PTROP,
 !     PLEV(2) = PTROP, PLEV(1) = ptop
 
 !     L_NLAYRAD is the number of radiation code layers
 !     L_NLEVRAD is the number of radiation code levels.  Level N is the
-!               top of layer N. 
+!               top of layer N.
 !
 !     L_NSPECTI is the number of IR spectral intervals
@@ -39 +39 @@
 !               GAUSS POINT 17 (aka the last one) is the special case
 !
-!     L_NPREF   is the number of reference pressures that the 
+!     L_NPREF   is the number of reference pressures that the
 !               k-coefficients are calculated on
 !     L_PINT    is the number of Lagrange interpolated reference
@@ -74 +74 @@
       real,    parameter :: L_TAUMAX  = 35
 
-      ! For Planck function integration: 
+      ! For Planck function integration:
       ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
       ! Integration boundary temperatures are NTstart/NTfac and Ntstop/NTfac
@@ -86 +86 @@
 
       ! Maximum number of grain size classes for aerosol convolution:
-      ! This must correspond to size of largest dataset used for aerosol 
+      ! This must correspond to size of largest dataset used for aerosol
       ! optical properties in datagcm folder.
-      integer, parameter :: nsizemax = 60
+      integer, parameter :: nsizemax = 200
 
       character(len=100),save :: corrkdir
@@ -105 +105 @@
   real*8, intent(in) :: tplanckmax
   real*8, intent(in) :: dtplanck
-  
+
   L_NLAYRAD = nbp_lev
   L_LEVELS  = 2*(nbp_lev-1)+3
@@ -113 +113 @@
   NTstart = int(tplanckmin * NTfac)
   NTstop  = int(tplanckmax * NTfac)
- 
+
   end subroutine