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

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

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

-revision 1.16 by molod,
Mon Mar 21 20:17:01 2005 UTC
+revision 1.21 by jmc,
Wed May  6 02:57:46 2009 UTC
 Line 5 
 C $Name$
         subroutine fizhi_init_vars (myThid)
  c-----------------------------------------------------------------------
  c  Routine to initialise the fizhi state.
  c
  c  Input: myThid       - Process number calling this routine
  c
  c  Notes:
  c   1) For a Cold Start -
  c      This routine takes the initial condition on the dynamics grid
  c      and interpolates to the physics grid to initialize the state
  c      variables that are on both grids. It initializes the variables
 Line 34 
 c---------------------------------------
  #include "SURFACE.h"
  #include "PARAMS.h"
  #include "chronos.h"
+ #ifdef ALLOW_EXCH2
+ #include "W2_EXCH2_TOPOLOGY.h"
+ #include "W2_EXCH2_PARAMS.h"
+ #endif /* ALLOW_EXCH2 */
         integer myThid
+        INTEGER xySize
+ #if defined(ALLOW_EXCH2)
+        PARAMETER ( xySize = W2_ioBufferSize )
+ #else
+        PARAMETER ( xySize = Nx*Ny )
+ #endif
+        Real*8 globalArr( xySize*8 )
  c pe on dynamics and physics grid refers to bottom edge
         _RL pephy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys+1,nSx,nSy)
         _RL pedyn(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr+1,nSx,nSy)
-Line 47 
 c pe on dynamics and physics grid refers
+Line 59 
 c pe on dynamics and physics grid refers
         integer i, j, L, bi, bj, Lbotij
         integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2
+        integer xsize, ysize
         logical alarm
         external alarm
+ #if defined(ALLOW_EXCH2)
+        xsize = exch2_global_Nx
+        ysize = exch2_global_Ny
+ #else
+        xsize = Nx
+        ysize = Ny
+ #endif
         im1 = 1-OLx
         im2 = sNx+OLx
         jm1 = 1-OLy
-Line 64 
 c    All Fizhi alarms must be on for the
+Line 84 
 c    All Fizhi alarms must be on for the
        if( .not.alarm('moist') .or. .not.alarm('turb')   .or.
       .    .not.alarm('radsw') .or. .not.alarm('radlw') ) then
-        print *,' Cant Start Fizhi experiment at ',nymd,' ',nhms
+        write(15,*) ' Cant Start Fizhi experiment at ',nymd,' ',nhms
         stop
        endif
  C Deal Here with Variables that are on a Fizhi Pickup or need Initialization
-       IF ( startTime.EQ.0. .AND. nIter0.EQ.0 ) THEN
+       IF ( startTime.EQ.baseTime .AND. nIter0.EQ.0 ) THEN
        print *,' In fizhi_init_vars: Beginning of New Experiment '
         do bj = myByLo(myThid), myByHi(myThid)
-Line 87 
 C Build pressures on dynamics grid
+Line 107 
 C Build pressures on dynamics grid
          do j = 1,sNy
          do i = 1,sNx
           Lbotij = ksurfC(i,j,bi,bj)
           if(Lbotij.ne.0.)
       .    pedyn(i,j,Lbotij,bi,bj) = Ro_surf(i,j,bi,bj) + etaH(i,j,bi,bj)
          enddo
          enddo
-Line 130 
 c   do units and get u = .025*ln(dP*10),
+Line 150 
 c   do units and get u = .025*ln(dP*10),
         enddo
         enddo
  c Create initial fields on phys. grid - Move Dynamics u and v to A-Grid
         call CtoA(myThid,uvel,vvel,maskW,maskS,im1,im2,jm1,jm2,Nr,
       .                     Nsx,Nsy,1,sNx,1,sNy,udyntemp,vdyntemp)
-Line 160 
 c         Physics works top-down. so ->
+Line 180 
 c         Physics works top-down. so ->
          enddo
          enddo
          call dyn2phys(theta,pedyn,im1,im2,jm1,jm2,Nr,Nsx,Nsy,
-      . 1,sNx,1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,0,tempphy)
+      . 1,sNx,1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,2,tempphy)
          do L = 1,Nrphys
          do j = 1,sNy
          do i = 1,sNx
-Line 168 
 c         Physics works top-down. so ->
+Line 188 
 c         Physics works top-down. so ->
          enddo
          enddo
          enddo
          call dyn2phys(salt,pedyn,im1,im2,jm1,jm2,Nr,Nsx,Nsy,
       . 1,sNx,1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,0,tempphy)
          do L = 1,Nrphys
-Line 203 
 c Zero out fizhi tendency arrays on the
+Line 222 
 c Zero out fizhi tendency arrays on the
          enddo
          enddo
- c Initialize tke, xlmt, khmt, xxmt, yymt, ctmt, zetamt,
+ c Initialize vegetation tile tke, xlmt, khmt, xxmt, yymt, ctmt, zetamt,
          if( (nhms.eq.nhms0) .and. (nymd.eq.nymd0) ) then
           print *,' Cold Start: Zero out Turb second moments '
           do i = 1,nchp
-Line 224 
 c Initialize tke, xlmt, khmt, xxmt, yymt
+Line 243 
 c Initialize tke, xlmt, khmt, xxmt, yymt
           stop
          endif
- c Now initialize land state too - tcanopy, etc... SET FOR NOW,
+ c Now initialize vegetation tile land state too - tcanopy, etc...
- c                                              READ CLIM FOR REAL
+         call fizhi_init_vegsurftiles( globalArr, xsize, ysize,
-         do i = 1,nchp
+      &                                nymd,nhms, 'D', myThid )
-         tcanopy(i,bi,bj) = 283.
-         tdeep(i,bi,bj) = 282.5
-         ecanopy(i,bi,bj) = 2.e-2
-         swetshal(i,bi,bj) = 0.6
-         swetroot(i,bi,bj) = 0.5
-         swetdeep(i,bi,bj) = 0.5
-         capac(i,bi,bj) = 0.
-         snodep(i,bi,bj) = 0.
-         enddo
+ c Now initialize fizhi arrays that will be on a pickup
          print *,' Initialize fizhi arrays that will be on pickup '
          imstturblw(bi,bj) = 0
          imstturbsw(bi,bj) = 0
-Line 276 
 c
+Line 287 
 c
        ELSE
        print *,' In fizhi_init_vars: Read from restart '
  C--   Read fizhi package state variables from pickup file
         call fizhi_read_pickup( nIter0, myThid )

 Legend:



Removed from v.1.16
 


changed lines


 
Added in v.1.21
 Legend:



Removed from v.1.16
 


changed lines


 
Added in v.1.21
-Removed from v.1.16
+Added in v.1.21

	ViewVC Help
Powered by ViewVC 1.1.22