UPP/cuparm_8f_source.html

  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.,   &


!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, &            !< Deficit Saturation Pressure (DSP) at cloud base - fast land ?

!    DSP0FL=-5875.E0, &            !< Deficit Saturation Pressure (DSP) at freezing level - fast land ?

!    DSPTFL=-1875.E0, &            !< Deficit Saturation Pressure (DSP) at cloud top - fast land ?

!    FSL=1.0E0, &                  !< Factor to obtain "slow" profiles over land ?

     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