--- 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