cuparm.f

  module cuparm_mod
!
    implicit none
!
    real,parameter :: H1=1.e0,h1d5=1.5e0,h2d5=2.5e0,h3000=3000.e0,       &
      h10e5=100000.e0,d00=0.e0,d125=.125e0,d50=.5e0,d608=.608e0,         &
      g=9.8e0,cp=1004.6e0,capa=0.28589641e0,rog=287.04/9.8,              &
      elwv=2.50e6,elivw=2.72e6,row=1.e3,epsq=2.e-12,                     &
      a2=17.2693882e0,a3=273.16e0,a4=35.86e0,                            &
      t0=273.16e0,t1=274.16e0,pq0=379.90516e0,stresh=1.10e0,             &
      stabs=1.0e0,stabd=.90e0,stabfc=1.00e0,dttop=0.0e0,                 &
!---VVVVV
      rhf=0.10,epsup=1.00,epsdn=1.05,epsth=0.0,                          &
      pbm=13000.,pqm=20000.,pno=1000.,pone=2500.,zsh=2000.,              &
      pfrz=15000.,pshu=45000.,                                           &
 
!    &, RHF=0.20,EPSUP=0.93,EPSDN=1.00,EPSTH=0.3
!    &, RHF=0.20,EPSUP=1.00,EPSDN=1.00,EPSTH=0.3
!AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
!
!    AUGUST '91: SCHEME HAVING THE OPTION OF USING DIFFERENT FAST AND
!    SLOW PROFILES FOR SEA AND FOR LAND POINTS; AND ALSO THE "SEA" AND
!    THE "LAND" SCHEME EVERYWHERE.  OVER LAND PROFILES DEPART FROM THE
!    FAST (DRY) PROFILES ONLY FOR PRECIPITATION/TIME STEP >
!    A PRESCRIBED VALUE (CURRENTLY, IN THE VERSION #3 DONE WEDNESDAY
!    18 SEPTEMBER, 1/4 INCH/24 H).  USE OF VARIOUS SWITCHES AS FOLLOWS.
!
!       THE "OLD" ("HARD", =ZAVISA OCT.1990) LAND SCHEME WITH FIXED
!       LAND PROFILES IS RUN BY
!       * SETTING OCT90=.TRUE. IN THE FIRST EXECUTABLE LINE FOLLOWING
!            THESE COMMENTS (THIS REACTIVATES EFI=H1 OVER LAND IF
!            .NOT.UNIS, AND CLDEFI(K)=EFIMN AS THE LEFTOVER CLDEFI
!            VALUE AT SWAP POINTS);
!       * DEFINING FAST LAND PROFILES SAME AS FAST SEA PROFILES (OR BY
!            CHOOSING ANOTHER SET OF LAND PROFILES ZAVISA USED AT
!            EARLIER TIME);
!       * SETTING FSL=1.; AND
!       * DEFINING STEFI (STARTING EFI) EQUAL TO AVGEFI.
!            (THE LAST THREE POINTS ARE HANDLED BY SWITCHING AROUND THE
!            "CFM" COMMENTS AT TWO PLACES)
!
!       THE "OLD,OLD" (APPR. ORIGINAL BETTS) SCHEME IS RUN BY
!       * SPECIFYING UNIL=.TRUE.;
!       * SETTING FSL=1.;
!       * SETTING OCT90=.TRUE.
!            (WITH THESE SETTINGS FAST LAND PROFILES ONLY ARE USED).
!                                                                     FM
     fss=.85e0,efimn=.20e0,efmnt=.70e0,fcc=.50,fcp=h1-fcc,              &
!
!         IN THIS VERSION 3.5, OVER LAND AND FOR THE FAST PROFILES, DSPB
!         IS PRESCRIBED TO BE 25 PERCENT DRIER THAN THE FAST SEA VALUE
!         (IN ROUGH AGREEMENT WITH BINDER, QJ, IN PRESS) WHILE DSP0 AND
!         DSPT ARE EACH 20 PERCENT DRIER THAN THE CORRESPONDING FAST
!         SEA VALUES.                WITH FSL=.875 THIS MAKES THE
!         AVERAGE OF THE FAST AND THE SLOW LAND PROFILES SOMEWHAT DRIER
!         THAN THE OCT90 FIXED LAND PROFILES.                         FM
!
     dspbfl=-4843.75e0,dsp0fl=-7050.00e0,dsptfl=-2250.0e0,fsl=.850e0,   &
!***   ACTIVATE THE FOLLOWING LINE IF OCT90=.TRUE. (AND COMMENT OUT THE
!***   PRECEDING LINE):
!    DSPBFL=-3875.E0,DSP0FL=-5875.E0,DSPTFL=-1875.E0,FSL=1.0E0,         &
     dspbfs=-3875.e0,dsp0fs=-5875.e0,dsptfs=-1875.e0,                   &
     dspbsl=dspbfl*fsl,dsp0sl=dsp0fl*fsl,dsptsl=dsptfl*fsl,             &
     dspbss=dspbfs*fss,dsp0ss=dsp0fs*fss,dsptss=dsptfs*fss,             &
!*** NEW CONVECTION SCHEME WITH CROSSING DSP PROFILES ******************
!+-  &, UNIS=.FALSE.,EFIMN=.71E0,EFMNT=.71,FCC=0.5,FCP=H1-FCC
!+-  &, DSPBL=-3875.E0,DSP0L=-5875.E0,DSPTL=-1875.E0
!+-  &, DSPBS=-2875.E0,DSP0S=-5125.E0,DSPTS=-4875.E0
!+-  &, DSPBF=-4375.E0,DSP0F=-4375.E0,DSPTF=-1000.E0
!*** BETTS CONVECTION SCHEME *******************************************
!    &, UNIS=.FALSE.,EFIMN=.9999E0,EFMNT=.9999E0,FCC=.50,FCP=H1-FCC
!    &, DSPBL=-3875.E0,DSP0L=-5875.E0,DSPTL=-1875.E0
!    &, DSPBF=-3875.E0,DSP0F=-5875.E0,DSPTF=-1875.E0
!    &, DSPBS=-3875.E0,DSP0S=-5875.E0,DSPTS=-1875.E0
!***********************************************************************
     trel=3000.,epsntp=.0010e0,efifc=5.0e0,                             &
     avgefi=(efimn+1.e0)*.5e0,dspc=-3000.e0,epsp=1.e-7,                 &
     stefi=1.e0,                                                        &
!*** ACTIVATE THE FOLLOWING LINE AND COMMENT OUT THE PRECEDING LINE IF
!*** OCT90=.TRUE.
!    &, STEFI=AVGEFI
     slopbl=(dspbfl-dspbsl)/(h1-efimn),                                 &
     slop0l=(dsp0fl-dsp0sl)/(h1-efimn),                                 &
     sloptl=(dsptfl-dsptsl)/(h1-efimn),                                 &
     slopbs=(dspbfs-dspbss)/(h1-efimn),                                 &
     slop0s=(dsp0fs-dsp0ss)/(h1-efimn),                                 &
     slopts=(dsptfs-dsptss)/(h1-efimn),                                 &
     slope=(h1   -efmnt)/(h1-efimn)
   real, parameter ::                                                   &
     A23M4L=a2*(a3-a4)*elwv,                                            &
     elocp=elivw/cp,cprlg=cp/(row*g*elwv),rcp=h1/cp
   logical,parameter ::                                                 &
     UNIS=.false.,unil=.false.,oct90=.false.                           
  end module cuparm_mod