/[MITgcm]/MITgcm/verification/fizhi-gridalt-hs/code/fizhi_init_vars.F
ViewVC logotype

Contents of /MITgcm/verification/fizhi-gridalt-hs/code/fizhi_init_vars.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.3 - (show annotations) (download)
Thu May 5 21:27:39 2005 UTC (18 years, 10 months ago) by molod
Branch: MAIN
CVS Tags: checkpoint58l_post, checkpoint57t_post, checkpoint57o_post, checkpoint58e_post, checkpoint57v_post, checkpoint58u_post, checkpoint58w_post, checkpoint57m_post, checkpoint57s_post, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint57k_post, checkpoint60, checkpoint61, checkpoint62, checkpoint63, checkpoint58r_post, checkpoint57i_post, checkpoint57y_post, checkpoint58n_post, checkpoint58x_post, checkpoint58t_post, checkpoint58h_post, checkpoint57y_pre, checkpoint58q_post, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint58j_post, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint59j, checkpoint57r_post, checkpoint59, checkpoint58, checkpoint57h_done, checkpoint58f_post, checkpoint57x_post, checkpoint58d_post, checkpoint58c_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint58a_post, checkpoint58i_post, checkpoint57q_post, checkpoint58g_post, checkpoint58o_post, checkpoint57z_post, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint58y_post, checkpoint58k_post, checkpoint58v_post, checkpoint58s_post, checkpoint61f, checkpoint61g, checkpoint61d, checkpoint61e, checkpoint61b, checkpoint61c, checkpoint58p_post, checkpoint61a, checkpoint61n, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61q, checkpoint57j_post, checkpoint61z, checkpoint61x, checkpoint61y, checkpoint58b_post, checkpoint57h_pre, checkpoint58m_post, checkpoint57l_post, checkpoint57h_post
Changes since 1.2: +153 -31 lines
Add this routine so as not to call fizhi_vegsurftiles routine - not needed for this exp.

1 C $Header: /u/gcmpack/MITgcm/pkg/fizhi/fizhi_init_vars.F,v 1.18 2005/05/04 22:15:17 molod Exp $
2 C $Name: $
3
4 #include "FIZHI_OPTIONS.h"
5 subroutine fizhi_init_vars (myThid)
6 c-----------------------------------------------------------------------
7 c Routine to initialise the fizhi state.
8 c
9 c Input: myThid - Process number calling this routine
10 c
11 c Notes:
12 c 1) For a Cold Start -
13 c This routine takes the initial condition on the dynamics grid
14 c and interpolates to the physics grid to initialize the state
15 c variables that are on both grids. It initializes the variables
16 c of the turbulence scheme to 0., and the land state from a model
17 c climatology.
18 c 2) For a Restart, read the fizhi pickup file
19 c 3) The velocity component physics fields are on an A-Grid
20 c
21 c Calls: dyn2phys (x4)
22 c-----------------------------------------------------------------------
23 implicit none
24 #include "SIZE.h"
25 #include "fizhi_SIZE.h"
26 #include "fizhi_land_SIZE.h"
27 #include "GRID.h"
28 #include "DYNVARS.h"
29 #include "gridalt_mapping.h"
30 #include "fizhi_coms.h"
31 #include "fizhi_land_coms.h"
32 #include "fizhi_earth_coms.h"
33 #include "EEPARAMS.h"
34 #include "SURFACE.h"
35 #include "PARAMS.h"
36 #include "chronos.h"
37
38 integer myThid
39
40 c pe on dynamics and physics grid refers to bottom edge
41 _RL pephy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys+1,nSx,nSy)
42 _RL pedyn(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr+1,nSx,nSy)
43 _RL windphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy)
44 _RL udyntemp(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
45 _RL vdyntemp(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
46 _RL tempphy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,nSx,nSy)
47
48 integer i, j, L, bi, bj, Lbotij
49 integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2
50 logical alarm
51 external alarm
52
53 im1 = 1-OLx
54 im2 = sNx+OLx
55 jm1 = 1-OLy
56 jm2 = sNy+OLy
57 idim1 = 1
58 idim2 = sNx
59 jdim1 = 1
60 jdim2 = sNy
61
62 c First Check to see if we can start a fizhi experiment at current time
63 c All Fizhi alarms must be on for the first time step of a segment
64
65 if( .not.alarm('moist') .or. .not.alarm('turb') .or.
66 . .not.alarm('radsw') .or. .not.alarm('radlw') ) then
67 print *,' Cant Start Fizhi experiment at ',nymd,' ',nhms
68 stop
69 endif
70
71 C Deal Here with Variables that are on a Fizhi Pickup or need Initialization
72
73 IF ( startTime.EQ.baseTime .AND. nIter0.EQ.0 ) THEN
74 print *,' In fizhi_init_vars: Beginning of New Experiment '
75
76 do bj = myByLo(myThid), myByHi(myThid)
77 do bi = myBxLo(myThid), myBxHi(myThid)
78
79 C Build pressures on dynamics grid
80 do j = 1,sNy
81 do i = 1,sNx
82 do L = 1,Nr
83 pedyn(i,j,L,bi,bj) = 0.
84 enddo
85 enddo
86 enddo
87 do j = 1,sNy
88 do i = 1,sNx
89 Lbotij = ksurfC(i,j,bi,bj)
90 if(Lbotij.ne.0.)
91 . pedyn(i,j,Lbotij,bi,bj) = Ro_surf(i,j,bi,bj) + etaH(i,j,bi,bj)
92 enddo
93 enddo
94 do j = 1,sNy
95 do i = 1,sNx
96 Lbotij = ksurfC(i,j,bi,bj)
97 do L = Lbotij+1,Nr+1
98 pedyn(i,j,L,bi,bj) = pedyn(i,j,L-1,bi,bj) -
99 . drF(L-1)*hfacC(i,j,L-1,bi,bj)
100 enddo
101 c Do not use a zero field as the top edge pressure for interpolation
102 if(pedyn(i,j,Nr+1,bi,bj).lt.1.e-5)
103 . pedyn(i,j,Nr+1,bi,bj) = 1.e-5
104 enddo
105 enddo
106 C Build pressures on physics grid
107 do j = 1,sNy
108 do i = 1,sNx
109 pephy(i,j,1,bi,bj)=Ro_surf(i,j,bi,bj) + etaH(i,j,bi,bj)
110 do L = 2,Nrphys+1
111 pephy(i,j,L,bi,bj)=pephy(i,j,L-1,bi,bj)-dpphys0(i,j,L-1,bi,bj)
112 enddo
113 c Do not use a zero field as the top edge pressure for interpolation
114 if(pephy(i,j,Nrphys+1,bi,bj).lt.1.e-5)
115 . pephy(i,j,Nrphys+1,bi,bj) = 1.e-5
116 enddo
117 enddo
118 c
119 c Create an initial wind magnitude field on the physics grid -
120 c Use a log wind law with z0=1cm, u*=1 cm/sec,
121 c do units and get u = .025*ln(dP*10), with dP in pa.
122 do L = 1,Nrphys
123 do j = 1,sNy
124 do i = 1,sNx
125 windphy(i,j,L,bi,bj) = 0.025 *
126 . log((pephy(i,j,1,bi,bj)-pephy(i,j,L+1,bi,bj))*10.)
127 enddo
128 enddo
129 enddo
130
131 enddo
132 enddo
133
134 c Create initial fields on phys. grid - Move Dynamics u and v to A-Grid
135 call CtoA(myThid,uvel,vvel,maskW,maskS,im1,im2,jm1,jm2,Nr,
136 . Nsx,Nsy,1,sNx,1,sNy,udyntemp,vdyntemp)
137
138 do bj = myByLo(myThid), myByHi(myThid)
139 do bi = myBxLo(myThid), myBxHi(myThid)
140
141 c Create initial fields on phys. grid - interpolate from dyn. grid
142 call dyn2phys(udyntemp,pedyn,im1,im2,jm1,jm2,Nr,Nsx,Nsy,
143 . 1,sNx,1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,1,tempphy)
144 c Note: Interpolation gives bottom-up arrays (level 1 is bottom),
145 c Physics works top-down. so -> need to flip arrays
146 do L = 1,Nrphys
147 do j = 1,sNy
148 do i = 1,sNx
149 uphy(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
150 enddo
151 enddo
152 enddo
153 call dyn2phys(vdyntemp,pedyn,im1,im2,jm1,jm2,Nr,Nsx,Nsy,
154 . 1,sNx,1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,1,tempphy)
155 do L = 1,Nrphys
156 do j = 1,sNy
157 do i = 1,sNx
158 vphy(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
159 enddo
160 enddo
161 enddo
162 call dyn2phys(theta,pedyn,im1,im2,jm1,jm2,Nr,Nsx,Nsy,
163 . 1,sNx,1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,0,tempphy)
164 do L = 1,Nrphys
165 do j = 1,sNy
166 do i = 1,sNx
167 thphy(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
168 enddo
169 enddo
170 enddo
171
172 call dyn2phys(salt,pedyn,im1,im2,jm1,jm2,Nr,Nsx,Nsy,
173 . 1,sNx,1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,0,tempphy)
174 do L = 1,Nrphys
175 do j = 1,sNy
176 do i = 1,sNx
177 sphy(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
178 enddo
179 enddo
180 enddo
181
182 c Zero out fizhi tendency arrays on the fizhi grid
183 do L = 1,Nrphys
184 do j = 1,sNy
185 do i = 1,sNx
186 duphy(i,j,L,bi,bj) = 0.
187 dvphy(i,j,L,bi,bj) = 0.
188 dthphy(i,j,L,bi,bj) = 0.
189 dsphy(i,j,L,bi,bj) = 0.
190 enddo
191 enddo
192 enddo
193
194 c Zero out fizhi tendency arrays on the dynamics grid
195 do L = 1,Nr
196 do j = jm1,jm2
197 do i = im1,im2
198 guphy(i,j,L,bi,bj) = 0.
199 gvphy(i,j,L,bi,bj) = 0.
200 gthphy(i,j,L,bi,bj) = 0.
201 gsphy(i,j,L,bi,bj) = 0.
202 enddo
203 enddo
204 enddo
205
206 c Initialize vegetation tile tke, xlmt, khmt, xxmt, yymt, ctmt, zetamt,
207 if( (nhms.eq.nhms0) .and. (nymd.eq.nymd0) ) then
208 print *,' Cold Start: Zero out Turb second moments '
209 do i = 1,nchp
210 ctmt(i,bi,bj) = 0.
211 xxmt(i,bi,bj) = 0.
212 yymt(i,bi,bj) = 0.
213 zetamt(i,bi,bj) = 0.
214 enddo
215 do L = 1,Nrphys
216 do i = 1,nchp
217 tke(i,L,bi,bj) = 0.
218 xlmt(i,L,bi,bj) = 0.
219 khmt(i,L,bi,bj) = 0.
220 enddo
221 enddo
222 else
223 print *,' Need initial Values for TKE - dont have them! '
224 stop
225 endif
226
227 c Now initialize vegetation tile land state too - tcanopy, etc...
228 c call fizhi_init_vegsurftiles( nymd,nhms, 'D', myThid )
229 c Now initialize land state too - tcanopy, etc... SET FOR NOW,
230 c READ CLIM FOR REAL
231 do i = 1,nchp
232 tcanopy(i,bi,bj) = 283.
233 tdeep(i,bi,bj) = 282.5
234 ecanopy(i,bi,bj) = 2.e-2
235 swetshal(i,bi,bj) = 0.6
236 swetroot(i,bi,bj) = 0.5
237 swetdeep(i,bi,bj) = 0.5
238 capac(i,bi,bj) = 0.
239 snodep(i,bi,bj) = 0.
240 enddo
241
242 c Now initialize fizhi arrays that will be on a pickup
243 print *,' Initialize fizhi arrays that will be on pickup '
244 imstturblw(bi,bj) = 0
245 imstturbsw(bi,bj) = 0
246 iras(bi,bj) = 0
247 nlwcld(bi,bj) = 0
248 nlwlz(bi,bj) = 0
249 nswcld(bi,bj) = 0
250 nswlz(bi,bj) = 0
251 do L = 1,Nrphys
252 do j = 1,sNy
253 do i = 1,sNx
254 swlz(i,j,L,bi,bj) = 0.
255 lwlz(i,j,L,bi,bj) = 0.
256 qliqavesw(i,j,L,bi,bj) = 0.
257 qliqavelw(i,j,L,bi,bj) = 0.
258 fccavesw(i,j,L,bi,bj) = 0.
259 fccavelw(i,j,L,bi,bj) = 0.
260 cldtot_sw(i,j,L,bi,bj) = 0.
261 cldras_sw(i,j,L,bi,bj) = 0.
262 cldlsp_sw(i,j,L,bi,bj) = 0.
263 cldtot_lw(i,j,L,bi,bj) = 0.
264 cldras_lw(i,j,L,bi,bj) = 0.
265 cldlsp_lw(i,j,L,bi,bj) = 0.
266 enddo
267 enddo
268 enddo
269 do j = 1,sNy
270 do i = 1,sNx
271 rainlsp(i,j,bi,bj) = 0.
272 raincon(i,j,bi,bj) = 0.
273 snowfall(i,j,bi,bj) = 0.
274 enddo
275 enddo
276
277 enddo
278 enddo
279
280 ELSE
281 print *,' In fizhi_init_vars: Read from restart '
282
283 C-- Read fizhi package state variables from pickup file
284
285 call fizhi_read_pickup( nIter0, myThid )
286 CALL FIZHI_READ_VEGTILES( nIter0, 'D', myThid )
287
288 ENDIF
289
290 return
291 end

  ViewVC Help
Powered by ViewVC 1.1.22