Parent Directory
|
Revision Log
|
Revision Graph
|
Patch
--- MITgcm_contrib/plumes/plumes_init_vars.F 2004/05/05 18:19:35 1.1
+++ MITgcm_contrib/plumes/plumes_init_vars.F 2004/05/13 22:20:26 1.2
@@ -8,10 +8,10 @@
c 1) For a Cold Start -
c This routine takes the initial condition on the dynamics grid
c and interpolates to the plumes grid to initialize the state
-c variables that are on both grids.
+c variables.
c 2) For a Restart, read the plumes pickup file
c
-c Calls: dyn2plumes (x4)
+c Calls: dyn2plume (x2)
c-----------------------------------------------------------------------
implicit none
#include "CPP_OPTIONS.h"
@@ -35,7 +35,7 @@
jdim1 = 1-OLy
jdim2 = sNy+OLy
idimout = Nxplumes
- jdimout = Nyplumes
+ Ldimout = Nrplumes
IF ( startTime.EQ.0. .AND. nIter0.EQ.0 ) THEN
print *,' In plumes_init_vars: Cold start '
@@ -48,14 +48,10 @@
do bi = myBxLo(myThid), myBxHi(myThid)
c Create initial fields on plumes grid - interpolate from dyn. grid
- call dyn2plumes(udyntemp,idim1,idim2,jdim1,jdim2,Nr,Nsx,Nsy,
- . 1,sNx,1,sNy,bi,bj,idimout,jdimout,Nr,Nsx,Nsy,uplume)
- call dyn2plumes(vdyntemp,idim1,idim2,jdim1,jdim2,Nr,Nsx,Nsy,
- . 1,sNx,1,sNy,bi,bj,idimout,jdimout,Nr,Nsx,Nsy,vplume)
- 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)
+ call dyn2plume(udyntemp,vdyntemp,idim1,idim2,jdim1,jdim2,
+ . 1,sNx,1,sNy,Nr,Nsx,Nsy,bi,bj,1,idimout,Ldimout,uplume)
+ call dyn2plume(theta,theta,idim1,idim2,jdim1,jdim2,
+ . 1,sNx,1,sNy,Nr,Nsx,Nsy,bi,bj,1,idimout,Ldimout,thplume)
enddo
enddo
| ViewVC Help | |
| Powered by ViewVC 1.1.22 |