|
subroutine, public | myjpbl (ntsd, me, dt_phs, epsl, epsq2, ht, stdh, dz, del, pmid, pinh, th, t, exner, q, cwm, u, v, tsk, qsfc, chklowq, thz0, qz0, uz0, vz0, xland, sice, snow, q2, exch_h, ustar, z0, el_myj, pblh, kpbl, ct, akhs, akms, elflx, mixht, thlm, qlm, rublten, rvblten, rthblten, rqblten, rqcblten, dusfc, dvsfc, dtsfc, dqsfc, xkzo, xkzmo, ict, ids, ide, jds, jde, ims, ime, jms, jme, its, ite, jts, jte, lm) |
|
subroutine | mixlen (lmh, rsqdt, u, v, thv, the, q2, epsl, epsq2, z, p, ps, rxner, gm, gh, el, pblh, lpbl, lmxl, ct, mixht, i, j, lm) |
|
subroutine | prodq2 (ntsd, me, lmh, dtturbl, ustar, gm, gh, el, q2, epsl, epsq2, i, j, lm) |
|
subroutine | difcof (ntsd, me, lmh, lmxl, gm, gh, el, t, q2, z, akm, akh, i, j, lm, print_diag, kpbl) |
|
subroutine | vdifq (lmh, dtdif, q2, el, z, i, j, lm) |
|
subroutine | vdifh (dtdif, lmh, thz0, qz0, rkhs, chklowq, ct, th, q, cwm, rkh, z, rho, i, j, lm) |
|
subroutine | vdifv (lmh, dtdif, uz0, vz0, rkms, d, u, v, rkm, z, rho, i, j, lm) |
|
subroutine, public | myjpbl_init (ids, ide, jds, jde, lm, ims, ime, jms, jme, its, ite, jts, jte) |
|
subroutine | herftbl |
|
real(kind=kfpt) function | herf (x) |
|
subroutine | tablept |
|
subroutine | tablett |
|
subroutine | spline (jtbl, nold, xold, yold, y2, nnew, xnew, ynew, p, q) |
|
|
real(kind=kfpt), parameter | a2 =17.2693882 |
|
real(kind=kfpt), parameter | a3 =273.15 |
|
real(kind=kfpt), parameter | a4 =35.86 |
|
real(kind=kfpt), parameter | cp =1004.6 |
|
real(kind=kfpt), parameter | eliv =2.850e6 |
|
real(kind=kfpt), parameter | elwv =2.501e6 |
|
real(kind=kfpt), parameter | r_v =461.6 |
|
real(kind=kfpt), parameter | epsq =1.e-12 |
|
real(kind=kfpt), parameter | g =9.8060226 |
|
real(kind=kfpt), parameter | pq0 =379.90516 |
|
real(kind=kfpt), parameter | r_d =287.04 |
|
real(kind=kfpt), parameter | ep_1 =R_V/R_D-1. |
|
real(kind=kfpt), parameter | p608 =R_V/R_D-1. |
|
real(kind=kfpt), parameter | pi =3.141592653589793 |
|
real(kind=kfpt), parameter | rhowater =1000. |
|
real(kind=kfpt), parameter | stbolt =5.67051E-8 |
|
real(kind=kfpt), parameter | cappa =R_D/CP |
|
real(kind=kfpt), parameter | eliwv =2.683e6 |
|
real(kind=kfpt), parameter | conw =1./G |
|
real(kind=kfpt), parameter | cont =CP/G |
|
real(kind=kfpt), parameter | conq =ELWV/G |
|
real(kind=kfpt), parameter | elevfc =0.6 |
|
real(kind=kfpt), parameter | vkarman =0.4 |
|
real(kind=kfpt), parameter | xls =ELIV |
|
real(kind=kfpt), parameter | xlv =ELWV |
|
real(kind=kfpt), parameter | rlivwv =XLS/XLV |
|
real(kind=kfpt), parameter | elocp =2.72E6/CP |
|
real(kind=kfpt), parameter | eps1 =1.E-12 |
|
real(kind=kfpt), parameter | eps2 =0. |
|
real(kind=kfpt), parameter | epsru =1.E-7 |
|
real(kind=kfpt), parameter | epsrs =1.E-7 |
|
real(kind=kfpt), parameter | epstrb =1.E-24 |
|
real(kind=kfpt), parameter | fh =1.10 |
|
real(kind=kfpt), parameter | alph =0.30 |
|
real(kind=kfpt), parameter | beta =1./273. |
|
real(kind=kfpt), parameter | el0max =1000. |
|
real(kind=kfpt), parameter | el0min =1. |
|
real(kind=kfpt), parameter | elfc =1. |
|
real(kind=kfpt), parameter | gam1 =0.2222222222222222222 |
|
real(kind=kfpt), parameter | a1 =0.659888514560862645 |
|
real(kind=kfpt), parameter | a2x =0.6574209922667784586 |
|
real(kind=kfpt), parameter | b1 =11.87799326209552761 |
|
real(kind=kfpt), parameter | b2 =7.226971804046074028 |
|
real(kind=kfpt), parameter | c1 =0.000830955950095854396 |
|
real(kind=kfpt), parameter | elz0 =0. |
|
real(kind=kfpt), parameter | esq =5.0 |
|
real(kind=kfpt), parameter | seafc =0.98 |
|
real(kind=kfpt), parameter | pq0sea =PQ0*SEAFC |
|
real(kind=kfpt), parameter | btg =BETA*G |
|
real(kind=kfpt), parameter | esqhf =0.5*5.0 |
|
real(kind=kfpt), parameter | rb1 =1./B1 |
|
real(kind=kfpt), parameter | adnh = 9.*A1*A2X*A2X*(12.*A1+3.*B2)*BTG*BTG |
|
real(kind=kfpt), parameter | adnm =18.*A1*A1*A2X*(B2-3.*A2X)*BTG |
|
real(kind=kfpt), parameter | anmh =-9.*A1*A2X*A2X*BTG*BTG |
|
real(kind=kfpt), parameter | anmm =-3.*A1*A2X*(3.*A2X+3.*B2*C1+18.*A1*C1-B2)*BTG |
|
real(kind=kfpt), parameter | bdnh = 3.*A2X*(7.*A1+B2)*BTG |
|
real(kind=kfpt), parameter | bdnm = 6.*A1*A1 |
|
real(kind=kfpt), parameter | beqh = A2X*B1*BTG+3.*A2X*(7.*A1+B2)*BTG |
|
real(kind=kfpt), parameter | beqm =-A1*B1*(1.-3.*C1)+6.*A1*A1 |
|
real(kind=kfpt), parameter | bnmh =-A2X*BTG |
|
real(kind=kfpt), parameter | bnmm =A1*(1.-3.*C1) |
|
real(kind=kfpt), parameter | bshh =9.*A1*A2X*A2X*BTG |
|
real(kind=kfpt), parameter | bshm =18.*A1*A1*A2X*C1 |
|
real(kind=kfpt), parameter | bsmh =-3.*A1*A2X*(3.*A2X+3.*B2*C1+12.*A1*C1-B2)*BTG |
|
real(kind=kfpt), parameter | cesh =A2X |
|
real(kind=kfpt), parameter | cesm =A1*(1.-3.*C1) |
|
real(kind=kfpt), parameter | cnv =EP_1*G/BTG |
|
real(kind=kfpt), parameter | aeqh =9.*A1*A2X*A2X*B1*BTG*BTG +9.*A1*A2X*A2X*(12.*A1+3.*B2)*BTG*BTG |
|
real(kind=kfpt), parameter | aeqm =3.*A1*A2X*B1*(3.*A2X+3.*B2*C1+18.*A1*C1-B2) *BTG+18.*A1*A1*A2X*(B2-3.*A2X)*BTG |
|
real(kind=kfpt), parameter | requ =-AEQH/AEQM |
|
real(kind=kfpt), parameter | epsgh =1.E-9 |
|
real(kind=kfpt), parameter | epsgm =REQU*EPSGH |
|
real(kind=kfpt), parameter | ubryl =(18.*REQU*A1*A1*A2X*B2*C1*BTG +9.*A1*A2X*A2X*B2*BTG*BTG) /(REQU*ADNM+ADNH) |
|
real(kind=kfpt), parameter | ubry =(1.+EPSRS)*UBRYL |
|
real(kind=kfpt), parameter | ubry3 =3.*UBRY |
|
real(kind=kfpt), parameter | aubh =27.*A1*A2X*A2X*B2*BTG*BTG-ADNH*UBRY3 |
|
real(kind=kfpt), parameter | aubm =54.*A1*A1*A2X*B2*C1*BTG -ADNM*UBRY3 |
|
real(kind=kfpt), parameter | bubh =(9.*A1*A2X+3.*A2X*B2)*BTG-BDNH*UBRY3 |
|
real(kind=kfpt), parameter | bubm =18.*A1*A1*C1 -BDNM*UBRY3 |
|
real(kind=kfpt), parameter | cubr =1. - UBRY3 |
|
real(kind=kfpt), parameter | rcubr =1./CUBR |
|
integer(kind=kint), parameter | itbl =401 |
|
integer(kind=kint), parameter | jtbl =1201 |
|
integer(kind=kint), parameter | kerfm =301 |
|
integer(kind=kint), parameter | kerfm2 =KERFM-2 |
|
real(kind=kfpt), parameter | pl =2500. |
|
real(kind=kfpt), parameter | ph =105000. |
|
real(kind=kfpt), parameter | thl =210. |
|
real(kind=kfpt), parameter | thh =365. |
|
real(kind=kfpt), parameter | xemin =0. |
|
real(kind=kfpt), parameter | xemax =3. |
|
real(kind=kfpt), save, private | rdp |
|
real(kind=kfpt), save, private | rdq |
|
real(kind=kfpt), save, private | rdth |
|
real(kind=kfpt), save, private | rdthe |
|
real(kind=kfpt), save, private | rdxe |
|
real(kind=kfpt), dimension(1:itbl), save, private | sthe |
|
real(kind=kfpt), dimension(1:itbl), save, private | the0 |
|
real(kind=kfpt), dimension(1:jtbl), save, private | qs0 |
|
real(kind=kfpt), dimension(1:jtbl), save, private | sqs |
|
real(kind=kfpt), dimension(1:kerfm), save, private | herff |
|
real(kind=kfpt), dimension(1:itbl, 1:jtbl), save, private | ptbl |
|
real(kind=kfpt), dimension(1:jtbl, 1:itbl), save, private | ttbl |
|