/[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.3 by molod, Mon Jun 14 20:34:50 2004 UTC revision 1.4 by molod, Fri Jul 23 22:32:28 2004 UTC
# Line 29  c--------------------------------------- Line 29  c---------------------------------------
29  #include "gridalt_mapping.h"  #include "gridalt_mapping.h"
30  #include "fizhi_coms.h"  #include "fizhi_coms.h"
31  #include "fizhi_land_coms.h"  #include "fizhi_land_coms.h"
32    #include "fizhi_earth_coms.h"
33  #include "EEPARAMS.h"  #include "EEPARAMS.h"
34  #include "SURFACE.h"  #include "SURFACE.h"
35  #include "PARAMS.h"  #include "PARAMS.h"
# Line 166  c         Physics works top-down. so -> Line 167  c         Physics works top-down. so ->
167  c Now initialize tke, xlmt, khmt, xxmt, yymt, ctmt, zetamt,  c Now initialize tke, xlmt, khmt, xxmt, yymt, ctmt, zetamt,
168  c Now initialize land state too - tcanopy, etc...  c Now initialize land state too - tcanopy, etc...
169          do L = 1,Nrphys          do L = 1,Nrphys
170          do i = 1,nchp          do i = 1,nchptot
171           tke(i,L,bi,bj) = 0.           tke(i,L,bi,bj) = 0.
172           xlmt(i,L,bi,bj) = 0.           xlmt(i,L,bi,bj) = 0.
173           khmt(i,L,bi,bj) = 0.           khmt(i,L,bi,bj) = 0.
174          enddo          enddo
175          enddo          enddo
176  c Now initialize land state too - tcanopy, etc... ZERO FOR NOW,  c Now initialize land state too - tcanopy, etc... SET FOR NOW,
177  c                                              READ CLIM FOR REAL  c                                              READ CLIM FOR REAL
178          do i = 1,nchp          do i = 1,nchptot
179          tcanopy(i,bi,bj) = 0.          tcanopy(i,bi,bj) = 283.
180          tdeep(i,bi,bj) = 0.          tdeep(i,bi,bj) = 282.5
181          ecanopy(i,bi,bj) = 0.          ecanopy(i,bi,bj) = 2.e-2
182          swetshal(i,bi,bj) = 0.          swetshal(i,bi,bj) = 0.6
183          swetroot(i,bi,bj) = 0.          swetroot(i,bi,bj) = 0.5
184          swetdeep(i,bi,bj) = 0.          swetdeep(i,bi,bj) = 0.5
185          capac(i,bi,bj) = 0.          capac(i,bi,bj) = 0.
186          snodep(i,bi,bj) = 0.          snodep(i,bi,bj) = 0.
187          enddo          enddo
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22