1 |
C $Header: /u/gcmpack/MITgcm/pkg/fizhi/fizhi_coms.h,v 1.7 2004/10/14 22:11:42 molod Exp $ |
2 |
C $Name: $ |
3 |
|
4 |
C The physics state uses the dynamics dimensions in the horizontal |
5 |
C and the land dimensions in the horizontal for turbulence variables |
6 |
C |
7 |
C Fizhi State Common - State variables on physics grid |
8 |
C ---------------------------------------------------------------------- |
9 |
COMMON /physics_state/ uphy,vphy,thphy,sphy, |
10 |
& ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke |
11 |
|
12 |
_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) |
14 |
_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) |
16 |
_RL ctmt(nchp,nSx,nSy), xxmt(nchp,nSx,nSy), yymt(nchp,nSx,nSy) |
17 |
_RL zetamt(nchp,nSx,nSy) |
18 |
_RL xlmt(nchp,Nrphys,nSx,nSy), khmt(nchp,Nrphys,nSx,nSy) |
19 |
_RL tke(nchp,Nrphys,nSx,nSy) |
20 |
|
21 |
C Fizhi Tendency Common - Changes in state variables due to physics |
22 |
C ---------------------------------------------------------------------- |
23 |
COMMON /physics_tendency/ duphy,dvphy,dthphy,dsphy |
24 |
|
25 |
_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) |
27 |
_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) |
29 |
|
30 |
C Fizhi Exports Common - Physics variables on other grids for export |
31 |
C ---------------------------------------------------------------------- |
32 |
COMMON /physics_exports/ guphy,gvphy,gthphy,gsphy |
33 |
|
34 |
_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) |
36 |
_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) |
38 |
|
39 |
C Fizhi Saver Common - Internal fizhi variables to be written on pickup |
40 |
C ---------------------------------------------------------------------- |
41 |
COMMON /physics_saver/ cldtot_lw,cldras_lw,cldlsp_lw,lwlz, |
42 |
& cldtot_sw,cldras_sw,cldlsp_sw,swlz, |
43 |
& qliqavesw,qliqavelw,fccavesw,fccavelw, |
44 |
& raincon,rainlsp,snowfall, |
45 |
& iras,nlwcld,nlwlz,nswcld,nswlz, |
46 |
& imstturbsw,imstturblw |
47 |
|
48 |
_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) |
50 |
_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) |
52 |
_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) |
54 |
_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) |
56 |
_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) |
58 |
_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) |
60 |
_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) |
62 |
_RL snowfall(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
63 |
INTEGER iras(nSx,nSy) |
64 |
INTEGER nlwcld(nSx,nSy),nlwlz(nSx,nSy) |
65 |
INTEGER nswcld(nSx,nSy),nswlz(nSx,nSy) |
66 |
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) |