/[MITgcm]/MITgcm/compare01/src/set_defaults.F
ViewVC logotype

Diff of /MITgcm/compare01/src/set_defaults.F

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

revision 1.2 by cnh, Mon Jun 8 21:42:59 1998 UTC revision 1.3 by cnh, Fri Jun 12 19:33:32 1998 UTC
# Line 85  C     \--------------------------------- Line 85  C     \---------------------------------
85        INTEGER K        INTEGER K
86    
87        REAL TAUMAX, DISTY, LY        REAL TAUMAX, DISTY, LY
88          Real*4 tmpXY(Nx,Ny)
89  C  C
90  C     /-------------------------------------------------------------------------\  C     /-------------------------------------------------------------------------\
91  C     | Parameter settings for small double gyre beta-plane experiment.         |  C     | Parameter settings for small double gyre beta-plane experiment.         |
# Line 103  CENDIGNOREENDIGNOREENDIGNOREENDIGNOREEND Line 104  CENDIGNOREENDIGNOREENDIGNOREENDIGNOREEND
104    
105    
106  CSTARTIGNORESTARTIGNORESTARTIGNORESTARTIGNORESTARTIGNORESTARTIGNORESTARTIGNORE  CSTARTIGNORESTARTIGNORESTARTIGNORESTARTIGNORESTARTIGNORESTARTIGNORESTARTIGNORE
107        phiMin      = -30.        phiMin      = -80.
108        dPhi        = 1.0        dPhi        = 4.0
109        dTheta      = 1.0        dTheta      = 4.0
110  CENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNORE  CENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNORE
111    
112    
# Line 183  C     State variables Line 184  C     State variables
184  CENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNORE  CENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNORE
185    
186    
187          THS      =  0.D0
188        THS(1)   = 20.D0        THS(1)   = 20.D0
189        THS(2)   = 10.D0        THS(2)   = 10.D0
190        THS(3)   =  8.D0        THS(3)   =  8.D0
# Line 201  C     T(_I3(2,:,:)) = 22.0D0 Line 203  C     T(_I3(2,:,:)) = 22.0D0
203  C     T(_I3(3,:,:)) = 21.0D0  C     T(_I3(3,:,:)) = 21.0D0
204  C     T(_I3(4,:,:)) = 20.0D0  C     T(_I3(4,:,:)) = 20.0D0
205        S = 10.0D0            S = 10.0D0    
206        DELTAZ(1) =  500.0D0        THS( 1)       = 16.0D0
207        DELTAZ(2) =  500.0D0        THS( 2)       = 15.2D0
208        DELTAZ(3) =  500.0D0        THS( 3)       = 14.5D0
209        DELTAZ(4) =  500.0D0        THS( 4)       = 13.9D0
210          THS( 5)       = 13.3D0
211          THS( 6)       = 12.4D0
212          THS( 7)       = 11.3D0
213          THS( 8)       =  9.9D0
214          THS( 9)       =  8.4D0
215          THS(10)       =  6.7D0
216          THS(11)       =  5.2D0
217          THS(12)       =  3.8D0
218          THS(13)       =  2.9D0
219          THS(14)       =  2.3D0
220          THS(15)       =  1.8D0
221          THS(16)       =  1.5D0
222          THS(17)       =  1.1D0
223          THS(18)       =  0.8D0
224          THS(19)       =  0.66D0
225          THS(20)       =  0.63D0
226          T(_I3( 1,:,:)) = THS( 1)
227          T(_I3( 2,:,:)) = THS( 2)
228          T(_I3( 3,:,:)) = THS( 3)
229          T(_I3( 4,:,:)) = THS( 4)
230          T(_I3( 5,:,:)) = THS( 5)
231          T(_I3( 6,:,:)) = THS( 6)
232          T(_I3( 7,:,:)) = THS( 7)
233          T(_I3( 8,:,:)) = THS( 8)
234          T(_I3( 9,:,:)) = THS( 9)
235          T(_I3(10,:,:)) = THS(10)
236          T(_I3(11,:,:)) = THS(11)
237          T(_I3(12,:,:)) = THS(12)
238          T(_I3(13,:,:)) = THS(13)
239          T(_I3(14,:,:)) = THS(14)
240          T(_I3(15,:,:)) = THS(15)
241          T(_I3(16,:,:)) = THS(16)
242          T(_I3(17,:,:)) = THS(17)
243          T(_I3(18,:,:)) = THS(18)
244          T(_I3(19,:,:)) = THS(19)
245          T(_I3(20,:,:)) = THS(20)
246    
247          SSPPT( 1)      = 34.65
248          SSPPT( 2)      = 34.75
249          SSPPT( 3)      = 34.82
250          SSPPT( 4)      = 34.87
251          SSPPT( 5)      = 34.90
252          SSPPT( 6)      = 34.90
253          SSPPT( 7)      = 34.86
254          SSPPT( 8)      = 34.78
255          SSPPT( 9)      = 34.69
256          SSPPT(10)      = 34.60
257          SSPPT(11)      = 34.58
258          SSPPT(12)      = 34.62
259          SSPPT(13)      = 34.68
260          SSPPT(14)      = 34.72
261          SSPPT(15)      = 34.73
262          SSPPT(16)      = 34.74
263          SSPPT(17)      = 34.73
264          SSPPT(18)      = 34.73
265          SSPPT(19)      = 34.72
266          SSPPT(20)      = 34.72
267          S(_I3( 1,:,:)) = SSPPT( 1)
268          S(_I3( 2,:,:)) = SSPPT( 2)
269          S(_I3( 3,:,:)) = SSPPT( 3)
270          S(_I3( 4,:,:)) = SSPPT( 4)
271          S(_I3( 5,:,:)) = SSPPT( 5)
272          S(_I3( 6,:,:)) = SSPPT( 6)
273          S(_I3( 7,:,:)) = SSPPT( 7)
274          S(_I3( 8,:,:)) = SSPPT( 8)
275          S(_I3( 9,:,:)) = SSPPT( 9)
276          S(_I3(10,:,:)) = SSPPT(10)
277          S(_I3(11,:,:)) = SSPPT(11)
278          S(_I3(12,:,:)) = SSPPT(12)
279          S(_I3(13,:,:)) = SSPPT(13)
280          S(_I3(14,:,:)) = SSPPT(14)
281          S(_I3(15,:,:)) = SSPPT(15)
282          S(_I3(16,:,:)) = SSPPT(16)
283          S(_I3(17,:,:)) = SSPPT(17)
284          S(_I3(18,:,:)) = SSPPT(18)
285          S(_I3(19,:,:)) = SSPPT(19)
286          S(_I3(20,:,:)) = SSPPT(20)
287          DELTAZ( 1) =  50.D0
288          DELTAZ( 2) =  50.D0
289          DELTAZ( 3) =  55.D0
290          DELTAZ( 4) =  60.D0
291          DELTAZ( 5) =  65.D0
292          DELTAZ( 6) =  70.D0
293          DELTAZ( 7) =  80.D0
294          DELTAZ( 8) =  95.D0
295          DELTAZ( 9) = 120.D0
296          DELTAZ(10) = 155.D0
297          DELTAZ(11) = 200.D0
298          DELTAZ(12) = 260.D0
299          DELTAZ(13) = 320.D0
300          DELTAZ(14) = 400.D0
301          DELTAZ(15) = 480.D0
302          DELTAZ(16) = 570.D0
303          DELTAZ(17) = 655.D0
304          DELTAZ(18) = 725.D0
305          DELTAZ(19) = 775.D0
306          DELTAZ(20) = 815.D0
307    
308  C     Geometry  C     Geometry
309        DM        = 20. D3        DM        = 20. D3
310        DM        = 1.        DM        = 4.
311        WRITE(20,*) ' DM',DM        WRITE(20,*) ' DM',DM
312    
313    
# Line 322  C     Wall at X = Lx Line 421  C     Wall at X = Lx
421        PMASK(_I3(:,nx,:))  = LAND  ! Land in east most cell        PMASK(_I3(:,nx,:))  = LAND  ! Land in east most cell
422  C     Wall at Y = Ly  C     Wall at Y = Ly
423        PMASK(_I3(:,:,ny))  = LAND  ! Land in northern most cell        PMASK(_I3(:,:,ny))  = LAND  ! Land in northern most cell
424  C     bathySource = textBathymetry        bathySource = binaryBathymetry
425  C     Island  C     Island
426        PMASK(_I3(Nz,1,24))  = LAND        PMASK(_I3(Nz,1,24))  = LAND
427    
# Line 352  C       fu(i,j) = -tauMax Line 451  C       fu(i,j) = -tauMax
451        ENDDO        ENDDO
452        fu(4,4) = fu(4,4)*0.917D0        fu(4,4) = fu(4,4)*0.917D0
453    
454          OPEN ( 11, FILE='windx', STATUS='old' , FORM='unformatted')
455          READ (11) tmpXY
456    C     Fu = tmpXY/DELPS(1)*G*0.1*uMask(_I3(1,:,:))           ! Dynes/cm**2 -> m/s**2
457          Fu = tmpXY/DELPS(1)*G*0.1
458          CLOSE(11)
459    
460          OPEN ( 11, FILE='windy', STATUS='old' , FORM='unformatted')
461          READ (11) tmpXY
462    C     Fu = tmpXY/DELPS(1)*G*0.1*uMask(_I3(1,:,:))           ! Dynes/cm**2 -> m/s**2
463          Fv = tmpXY/DELPS(1)*G*0.1
464          CLOSE(11)
465    
466  C  C
467  CcnhDebugStarts  CcnhDebugStarts
468  Cdbg  WRITE(0,*) ' lY    = ', lY  Cdbg  WRITE(0,*) ' lY    = ', lY
# Line 359  Cdbg  WRITE(0,*) ' distY = ', distY Line 470  Cdbg  WRITE(0,*) ' distY = ', distY
470  Cdbg  WRITE(0,*) ' tauMax = ', tauMax  Cdbg  WRITE(0,*) ' tauMax = ', tauMax
471        WRITE(0,*) ' fu '        WRITE(0,*) ' fu '
472        CALL PLOT_FIELD(fu,Nx,Ny)        CALL PLOT_FIELD(fu,Nx,Ny)
473  Cdbg  STOP        WRITE(0,*) ' fv '
474          CALL PLOT_FIELD(fv,Nx,Ny)
475    C     STOP
476  CcnhDebugEnds  CcnhDebugEnds
477  C  C
478        END        END
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22