1 |
C $Header$ |
C $Header$ |
2 |
C $Name$ |
C $Name$ |
3 |
|
|
|
|
|
4 |
C The physics state uses the dynamics dimensions in the horizontal |
C The physics state uses the dynamics dimensions in the horizontal |
5 |
C and the land dimensions in the horizontal for turbulence variables |
C and the land dimensions in the horizontal for turbulence variables |
6 |
c |
C |
7 |
c Fizhi State Common - State variables on physics grid |
C Fizhi State Common - State variables on physics grid |
8 |
c ---------------------------------------------------------------------- |
C ---------------------------------------------------------------------- |
9 |
common /physics_state/ uphy,vphy,thphy,sphy, |
COMMON /physics_state/ uphy,vphy,thphy,sphy, |
10 |
. ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke |
& ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke |
11 |
|
|
12 |
_RL uphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL uphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
13 |
_RL vphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL vphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
14 |
_RL thphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL thphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
15 |
_RL sphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL sphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
16 |
_RL ctmt(nchp,Nsx,Nsy), xxmt(nchp,Nsx,Nsy), yymt(nchp,Nsx,Nsy) |
_RL ctmt(nchp,nSx,nSy), xxmt(nchp,nSx,nSy), yymt(nchp,nSx,nSy) |
17 |
_RL zetamt(nchp,Nsx,Nsy) |
_RL zetamt(nchp,nSx,nSy) |
18 |
_RL xlmt(nchp,Nrphys,Nsx,Nsy), khmt(nchp,Nrphys,Nsx,Nsy) |
_RL xlmt(nchp,Nrphys,nSx,nSy), khmt(nchp,Nrphys,nSx,nSy) |
19 |
_RL tke(nchp,Nrphys,Nsx,Nsy) |
_RL tke(nchp,Nrphys,nSx,nSy) |
20 |
|
|
21 |
c Fizhi Tendency Common - Changes in state variables due to physics |
C Fizhi Tendency Common - Changes in state variables due to physics |
22 |
c ---------------------------------------------------------------------- |
C ---------------------------------------------------------------------- |
23 |
common /physics_tendency/ duphy,dvphy,dthphy,dsphy |
COMMON /physics_tendency/ duphy,dvphy,dthphy,dsphy |
24 |
|
|
25 |
_RL duphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL duphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
26 |
_RL dvphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL dvphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
27 |
_RL dthphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL dthphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
28 |
_RL dsphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL dsphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
29 |
|
|
30 |
c Fizhi Exports Common - Physics variables on other grids for export |
C Fizhi Exports Common - Physics variables on other grids for export |
31 |
c ---------------------------------------------------------------------- |
C ---------------------------------------------------------------------- |
32 |
common /physics_exports/ guphy,gvphy,gthphy,gsphy |
COMMON /physics_exports/ guphy,gvphy,gthphy,gsphy |
33 |
|
|
34 |
_RL guphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,Nsx,Nsy) |
_RL guphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
35 |
_RL gvphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,Nsx,Nsy) |
_RL gvphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
36 |
_RL gthphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,Nsx,Nsy) |
_RL gthphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
37 |
_RL gsphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,Nsx,Nsy) |
_RL gsphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
38 |
|
|
39 |
c Fizhi Saver Common - Internal fizhi variables to be written on pickup |
C Fizhi Saver Common - Internal fizhi variables to be written on pickup |
40 |
c ---------------------------------------------------------------------- |
C ---------------------------------------------------------------------- |
41 |
common /physics_saver/ cldtot_lw,cldras_lw,cldlsp_lw,lwlz, |
COMMON /physics_saver/ cldtot_lw,cldras_lw,cldlsp_lw,lwlz, |
42 |
. cldtot_sw,cldras_sw,cldlsp_sw,swlz, |
& cldtot_sw,cldras_sw,cldlsp_sw,swlz, |
43 |
. qliqavesw,qliqavelw,fccavesw,fccavelw, |
& qliqavesw,qliqavelw,fccavesw,fccavelw, |
44 |
. raincon,rainlsp,snowfall, |
& raincon,rainlsp,snowfall, |
45 |
. iras,nlwcld,nlwlz,nswcld,nswlz, |
& iras,nlwcld,nlwlz,nswcld,nswlz, |
46 |
. imstturbsw,imstturblw |
& imstturbsw,imstturblw |
47 |
|
|
48 |
_RL cldtot_lw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL cldtot_lw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
49 |
_RL cldras_lw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL cldras_lw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
50 |
_RL cldlsp_lw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL cldlsp_lw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
51 |
_RL lwlz(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL lwlz(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
52 |
_RL cldtot_sw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL cldtot_sw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
53 |
_RL cldras_sw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL cldras_sw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
54 |
_RL cldlsp_sw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL cldlsp_sw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
55 |
_RL swlz(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL swlz(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
56 |
_RL qliqavesw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL qliqavesw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
57 |
_RL qliqavelw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL qliqavelw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
58 |
_RL fccavesw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL fccavesw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
59 |
_RL fccavelw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy) |
_RL fccavelw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy) |
60 |
_RL raincon(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nsx,Nsy) |
_RL raincon(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
61 |
_RL rainlsp(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nsx,Nsy) |
_RL rainlsp(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
62 |
_RL snowfall(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nsx,Nsy) |
_RL snowfall(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
63 |
integer iras(Nsx,Nsy) |
INTEGER iras(nSx,nSy) |
64 |
integer nlwcld(Nsx,Nsy),nlwlz(Nsx,Nsy) |
INTEGER nlwcld(nSx,nSy),nlwlz(nSx,nSy) |
65 |
integer nswcld(Nsx,Nsy),nswlz(Nsx,Nsy) |
INTEGER nswcld(nSx,nSy),nswlz(nSx,nSy) |
66 |
integer imstturbsw(Nsx,Nsy),imstturblw(Nsx,Nsy) |
INTEGER imstturbsw(nSx,nSy),imstturblw(nSx,nSy) |
67 |
|
|
68 |
|
C ---------------------------------------------------------------------- |
69 |
|
C turbStart :: true when doing a cold-start for turbulence |
70 |
|
COMMON /FIZHI_START/ turbStart |
71 |
|
LOGICAL turbStart(nSx,nSy) |