40 |
C fit exactly into the input (dynamics) grid levels |
C fit exactly into the input (dynamics) grid levels |
41 |
C*********************************************************************** |
C*********************************************************************** |
42 |
implicit none |
implicit none |
43 |
|
cinterp1 #include "PACKAGES_CONFIG.h" |
44 |
#include "CPP_OPTIONS.h" |
#include "CPP_OPTIONS.h" |
45 |
|
|
46 |
integer im1, im2, jm1, jm2, lmdyn, lmphy, Nsx, Nsy, flg |
integer im1, im2, jm1, jm2, lmdyn, lmphy, Nsx, Nsy, flg |
57 |
_RL pphy(im1:im2,jm1:jm2,lmphy) |
_RL pphy(im1:im2,jm1:jm2,lmphy) |
58 |
_RL dpkedyn, dpkephy, windsum, qd |
_RL dpkedyn, dpkephy, windsum, qd |
59 |
integer i,j,L,Lout1,Lout2,Lphy |
integer i,j,L,Lout1,Lout2,Lphy |
60 |
_RL getcon, kappa |
cinterp1 _RL kappa |
61 |
|
#ifdef ALLOW_FIZHI |
62 |
|
cinterp1 _RL getcon |
63 |
|
#else |
64 |
|
cinterp1 #include 'SIZE.h' |
65 |
|
cinterp1 #include 'EEPARAMS.h' |
66 |
|
cinterp1 #include 'PARAMS.h' |
67 |
|
#endif |
68 |
|
|
69 |
kappa = getcon('KAPPA') |
#ifdef ALLOW_FIZHI |
70 |
|
cinterp1 kappa = getcon('KAPPA') |
71 |
|
#else |
72 |
|
cinterp1 kappa = atm_kappa |
73 |
|
#endif |
74 |
|
|
75 |
C define physics grid mid level pressures |
C define physics grid mid level pressures |
76 |
do Lphy = 1,lmphy |
do Lphy = 1,lmphy |