/[MITgcm]/MITgcm/pkg/fizhi/fizhi_init_vars.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/fizhi_init_vars.F

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

revision 1.8 by molod, Thu Aug 12 15:21:22 2004 UTC revision 1.12 by molod, Mon Mar 21 18:37:33 2005 UTC
# Line 33  c--------------------------------------- Line 33  c---------------------------------------
33  #include "EEPARAMS.h"  #include "EEPARAMS.h"
34  #include "SURFACE.h"  #include "SURFACE.h"
35  #include "PARAMS.h"  #include "PARAMS.h"
36    #include "chronos.h"
37    
38         integer myThid         integer myThid
39    
# Line 56  c pe on dynamics and physics grid refers Line 57  c pe on dynamics and physics grid refers
57         jdim1 = 1         jdim1 = 1
58         jdim2 = sNy         jdim2 = sNy
59    
60    C Deal Here with Variables that are on a Fizhi Pickup or need Initialization
61    
62        IF ( startTime.EQ.0. .AND. nIter0.EQ.0 ) THEN        IF ( startTime.EQ.0. .AND. nIter0.EQ.0 ) THEN
63        print *,' In fizhi_init_vars: Cold start '        print *,' In fizhi_init_vars: Beginning of New Experiment '
64    
65         do bj = myByLo(myThid), myByHi(myThid)         do bj = myByLo(myThid), myByHi(myThid)
66         do bi = myBxLo(myThid), myBxHi(myThid)         do bi = myBxLo(myThid), myBxHi(myThid)
# Line 165  c         Physics works top-down. so -> Line 168  c         Physics works top-down. so ->
168          enddo          enddo
169          enddo          enddo
170    
171  c Now initialize tke, xlmt, khmt, xxmt, yymt, ctmt, zetamt,  c Zero out fizhi tendency arrays on the fizhi grid
 c Now initialize land state too - tcanopy, etc...  
172          do L = 1,Nrphys          do L = 1,Nrphys
173          do i = 1,nchptot(bi,bj)          do j = 1,sNy
174           tke(i,L,bi,bj) = 0.          do i = 1,sNx
175           xlmt(i,L,bi,bj) = 0.           duphy(i,j,L,bi,bj) = 0.
176           khmt(i,L,bi,bj) = 0.           dvphy(i,j,L,bi,bj) = 0.
177             dthphy(i,j,L,bi,bj) = 0.
178             dsphy(i,j,L,bi,bj) = 0.
179            enddo
180            enddo
181            enddo
182    
183    c Zero out fizhi tendency arrays on the dynamics grid
184            do L = 1,Nr
185            do j = jm1,jm2
186            do i = im1,im2
187             guphy(i,j,L,bi,bj) = 0.
188             gvphy(i,j,L,bi,bj) = 0.
189             gthphy(i,j,L,bi,bj) = 0.
190             gsphy(i,j,L,bi,bj) = 0.
191          enddo          enddo
192          enddo          enddo
193            enddo
194    
195    c Initialize tke, xlmt, khmt, xxmt, yymt, ctmt, zetamt,
196            if( (nhms.eq.nhms0) .and. (nymd.eq.nymd0) ) then
197             print *,' Cold Start: Zero out Turb second moments '
198             do i = 1,nchp
199              ctmt(i,bi,bj) = 0.
200              xxmt(i,bi,bj) = 0.
201              yymt(i,bi,bj) = 0.
202              zetamt(i,bi,bj) = 0.
203             enddo
204             do L = 1,Nrphys
205             do i = 1,nchp
206              tke(i,L,bi,bj) = 0.
207              xlmt(i,L,bi,bj) = 0.
208              khmt(i,L,bi,bj) = 0.
209             enddo
210             enddo
211            else
212             print *,' Need initial Values for TKE - don't have them! '
213             stop
214            endif
215    
216  c Now initialize land state too - tcanopy, etc... SET FOR NOW,  c Now initialize land state too - tcanopy, etc... SET FOR NOW,
217  c                                              READ CLIM FOR REAL  c                                              READ CLIM FOR REAL
218          do i = 1,nchptot(bi,bj)          do i = 1,nchp
219          tcanopy(i,bi,bj) = 283.          tcanopy(i,bi,bj) = 283.
220          tdeep(i,bi,bj) = 282.5          tdeep(i,bi,bj) = 282.5
221          ecanopy(i,bi,bj) = 2.e-2          ecanopy(i,bi,bj) = 2.e-2
# Line 187  c Line 226  c
226          snodep(i,bi,bj) = 0.          snodep(i,bi,bj) = 0.
227          enddo          enddo
228    
229            print *,' Initialize fizhi arrays that will be on pickup '
230            imstturblw(bi,bj) = 0
231            imstturbsw(bi,bj) = 0
232            iras(bi,bj) = 0
233            nlwcld(bi,bj) = 0
234            nlwlz(bi,bj) = 0
235            nswcld(bi,bj) = 0
236            nswlz(bi,bj) = 0
237            do L = 1,Nrphys
238            do j = 1,sNy
239            do i = 1,sNx
240             swlz(i,j,L) = 0.
241             lwlz(i,j,L) = 0.
242             qliqavesw(i,j,L) = 0.
243             qliqavelw(i,j,L) = 0.
244             fccavesw(i,j,L) = 0.
245             fccavelw(i,j,L) = 0.
246             cldtot_sw(i,j,L) = 0.
247             cldras_sw(i,j,L) = 0.
248             cldlsp_sw(i,j,L) = 0.
249             cldtot_lw(i,j,L) = 0.
250             cldras_lw(i,j,L) = 0.
251             cldlsp_lw(i,j,L) = 0.
252             lwdt(i,j,L,bi,bj) = 0.
253             swdt(i,j,L,bi,bj) = 0.
254             turbt(i,j,L,bi,bj) = 0.
255             moistt(i,j,L,bi,bj) = 0.
256             turbq(i,j,L,1,bi,bj) = 0.
257             moistq(i,j,L,1,bi,bj) = 0.
258             turbu(i,j,L,bi,bj) = 0.
259             moistu(i,j,L,bi,bj) = 0.
260             turbv(i,j,L,bi,bj) = 0.
261             moistv(i,j,L,bi,bj) = 0.
262            enddo
263            enddo
264            enddo
265            do j = 1,sNy
266            do i = 1,sNx
267             rainlsp(i,j) = 0.
268             raincon(i,j) = 0.
269             snowfall(i,j) = 0.
270            enddo
271            enddo
272    
273         enddo         enddo
274         enddo         enddo
275    
# Line 194  c Line 277  c
277        print *,' In fizhi_init_vars: Read from restart '        print *,' In fizhi_init_vars: Read from restart '
278                                                                                                                                                                        
279  C--   Read fizhi package state variables from pickup file  C--   Read fizhi package state variables from pickup file
280    
281         call fizhi_read_pickup( nIter0, myThid )         call fizhi_read_pickup( nIter0, myThid )
282         CALL FIZHI_READ_VEGTILES( nIter0, 'D', myThid )         CALL FIZHI_READ_VEGTILES( nIter0, 'D', myThid )
283    
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22