/[MITgcm]/MITgcm_contrib/plumes/plumes_init_vars.F
ViewVC logotype

Diff of /MITgcm_contrib/plumes/plumes_init_vars.F

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

revision 1.1 by molod, Wed May 5 18:19:35 2004 UTC revision 1.2 by molod, Thu May 13 22:20:26 2004 UTC
# Line 8  c  Notes: Line 8  c  Notes:
8  c   1) For a Cold Start -  c   1) For a Cold Start -
9  c      This routine takes the initial condition on the dynamics grid  c      This routine takes the initial condition on the dynamics grid
10  c      and interpolates to the plumes grid to initialize the state  c      and interpolates to the plumes grid to initialize the state
11  c      variables that are on both grids.  c      variables.
12  c   2) For a Restart, read the plumes pickup file  c   2) For a Restart, read the plumes pickup file
13  c  c
14  c Calls: dyn2plumes (x4)  c Calls: dyn2plume (x2)
15  c-----------------------------------------------------------------------  c-----------------------------------------------------------------------
16         implicit none         implicit none
17  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
# Line 35  c--------------------------------------- Line 35  c---------------------------------------
35         jdim1 = 1-OLy         jdim1 = 1-OLy
36         jdim2 = sNy+OLy         jdim2 = sNy+OLy
37         idimout = Nxplumes         idimout = Nxplumes
38         jdimout = Nyplumes         Ldimout = Nrplumes
39    
40        IF ( startTime.EQ.0. .AND. nIter0.EQ.0 ) THEN        IF ( startTime.EQ.0. .AND. nIter0.EQ.0 ) THEN
41        print *,' In plumes_init_vars: Cold start '        print *,' In plumes_init_vars: Cold start '
# Line 48  c Create initial fields on phys. grid - Line 48  c Create initial fields on phys. grid -
48         do bi = myBxLo(myThid), myBxHi(myThid)         do bi = myBxLo(myThid), myBxHi(myThid)
49    
50  c Create initial fields on plumes grid - interpolate from dyn. grid  c Create initial fields on plumes grid - interpolate from dyn. grid
51          call dyn2plumes(udyntemp,idim1,idim2,jdim1,jdim2,Nr,Nsx,Nsy,          call dyn2plume(udyntemp,vdyntemp,idim1,idim2,jdim1,jdim2,
52       .    1,sNx,1,sNy,bi,bj,idimout,jdimout,Nr,Nsx,Nsy,uplume)       . 1,sNx,1,sNy,Nr,Nsx,Nsy,bi,bj,1,idimout,Ldimout,uplume)
53          call dyn2plumes(vdyntemp,idim1,idim2,jdim1,jdim2,Nr,Nsx,Nsy,          call dyn2plume(theta,theta,idim1,idim2,jdim1,jdim2,
54       .    1,sNx,1,sNy,bi,bj,idimout,jdimout,Nr,Nsx,Nsy,vplume)       . 1,sNx,1,sNy,Nr,Nsx,Nsy,bi,bj,1,idimout,Ldimout,thplume)
         call dyn2plumes(theta,idim1,idim2,jdim1,jdim2,Nr,Nsx,Nsy,  
      .    1,sNx,1,sNy,bi,bj,idimout,jdimout,Nr,Nsx,Nsy,thplume)  
         call dyn2plumes(salt,idim1,idim2,jdim1,jdim2,Nr,Nsx,Nsy,  
      .    1,sNx,1,sNy,bi,bj,idimout,jdimout,Nr,Nsx,Nsy,splume)  
55    
56         enddo         enddo
57         enddo         enddo
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22