/[MITgcm]/MITgcm_contrib/gael/pkg/smooth/smooth_ensinv.F
ViewVC logotype

Diff of /MITgcm_contrib/gael/pkg/smooth/smooth_ensinv.F

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

revision 1.1 by gforget, Tue Jun 19 18:23:18 2007 UTC revision 1.2 by gforget, Fri Oct 16 03:36:34 2009 UTC
# Line 56  c#include "adcommon.h" Line 56  c#include "adcommon.h"
56       & 'wvar_T_4_50levels.bin'       & 'wvar_T_4_50levels.bin'
57        call mdsreadfield(fnamegeneric,32,'RL',nR,        call mdsreadfield(fnamegeneric,32,'RL',nR,
58       & wthetalev,1,mythid)       & wthetalev,1,mythid)
59        _EXCH_XYZ_R8( wthetalev, mythid )        _EXCH_XYZ_RL( wthetalev, mythid )
60        write(fnamegeneric(1:80),'(1a)')        write(fnamegeneric(1:80),'(1a)')
61       & 'wvar_S_4_50levels.bin'       & 'wvar_S_4_50levels.bin'
62        call mdsreadfield(fnamegeneric,32,'RL',nR,        call mdsreadfield(fnamegeneric,32,'RL',nR,
63       & wsaltlev,1,mythid)       & wsaltlev,1,mythid)
64        _EXCH_XYZ_R8( wsaltlev, mythid )        _EXCH_XYZ_RL( wsaltlev, mythid )
65    
66  C initialize the random number generator:  C initialize the random number generator:
67        theta(1,1,1,1,1)=port_rand(1)        theta(1,1,1,1,1)=port_rand(1)
# Line 104  C define a perturbation: Line 104  C define a perturbation:
104          ENDDO          ENDDO
105         ENDDO         ENDDO
106        ENDDO        ENDDO
107        _EXCH_XYZ_R8 ( theta, myThid )        _EXCH_XYZ_RL ( theta, myThid )
108        _EXCH_XYZ_R8 ( salt, myThid )        _EXCH_XYZ_RL ( salt, myThid )
109    
110        if (ii.GT.nbRandPrev) then        if (ii.GT.nbRandPrev) then
111    
# Line 137  C scale the variance Line 137  C scale the variance
137          ENDDO          ENDDO
138         ENDDO         ENDDO
139        ENDDO        ENDDO
140        _EXCH_XYZ_R8 ( theta, myThid )        _EXCH_XYZ_RL ( theta, myThid )
141        _EXCH_XYZ_R8 ( salt, myThid )        _EXCH_XYZ_RL ( salt, myThid )
142    
143  c do the profiles_inloop operations  c do the profiles_inloop operations
144        do iloop=1,nTimeSteps        do iloop=1,nTimeSteps
# Line 208  c initialise adtheta/adsalt and do adinl Line 208  c initialise adtheta/adsalt and do adinl
208          ENDDO          ENDDO
209         ENDDO         ENDDO
210        ENDDO        ENDDO
211        _EXCH_XYZ_R8 ( adtheta, myThid )        _EXCH_XYZ_RL ( adtheta, myThid )
212        _EXCH_XYZ_R8 ( adsalt, myThid )        _EXCH_XYZ_RL ( adsalt, myThid )
213        do iloop = ntimesteps, 1, -1        do iloop = ntimesteps, 1, -1
214          mytime = starttime+float(iloop-1)*deltatclock          mytime = starttime+float(iloop-1)*deltatclock
215          call adprofiles_inloop( mytime,mythid )          call adprofiles_inloop( mytime,mythid )
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22