85 |
INTEGER K |
INTEGER K |
86 |
|
|
87 |
REAL TAUMAX, DISTY, LY |
REAL TAUMAX, DISTY, LY |
88 |
Real*4 tmpXY(Nx,Ny) |
Real*4 tmpXY(Nx,Ny) |
89 |
|
Real*4 tmpXYZ( _I3(Nz,Nx,Ny) ) |
90 |
C |
C |
91 |
C /-------------------------------------------------------------------------\ |
C /-------------------------------------------------------------------------\ |
92 |
C | Parameter settings for small double gyre beta-plane experiment. | |
C | Parameter settings for small double gyre beta-plane experiment. | |
97 |
CP = 3900.D0 |
CP = 3900.D0 |
98 |
CENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNORE |
CENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNORE |
99 |
|
|
100 |
delT = 1200.D0 |
delT = 3600.D0 |
101 |
freeSurfFac = 1.D0 |
freeSurfFac = 1.D0 |
102 |
startTime = 0. |
startTime = 0. |
103 |
endTime = startTime+10.*delT |
endTime = startTime+10.*delT |
111 |
CENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNORE |
CENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNOREENDIGNORE |
112 |
|
|
113 |
|
|
114 |
a2MomXY = 4.D2 |
a2MomXY = 5.0D5 |
115 |
a2MomZ = 1.D-2*G*RONIL*G*RONIL |
a2MomZ = 1.D-3*G*RONIL*G*RONIL |
116 |
|
C a2MomXY = 0.0D5 |
117 |
|
C a2MomZ = 0.D-3*G*RONIL*G*RONIL |
118 |
a4MomXY = 0.D0 |
a4MomXY = 0.D0 |
119 |
a2TempXY = 4.D2 |
a2TempXY = 1.D3 |
120 |
a2TempZ = 1.D-2*G*RONIL*G*RONIL |
a2TempZ = 3.D-5*G*RONIL*G*RONIL |
121 |
|
a2SaltXY = 1.D3 |
122 |
|
a2SaltZ = 3.D-5*G*RONIL*G*RONIL |
123 |
|
C a2SaltXY = 0.D3 |
124 |
|
C a2SaltZ = 0.D-5*G*RONIL*G*RONIL |
125 |
C a2TempZ = 1.D-5*G*RONIL*G*RONIL |
C a2TempZ = 1.D-5*G*RONIL*G*RONIL |
126 |
a4TempXy = 0.D0 |
a4TempXy = 0.D0 |
127 |
|
a4SaltXy = 0.D0 |
128 |
|
|
129 |
|
|
130 |
CSTARTIGNORESTARTIGNORESTARTIGNORESTARTIGNORESTARTIGNORESTARTIGNORESTARTIGNORE |
CSTARTIGNORESTARTIGNORESTARTIGNORESTARTIGNORESTARTIGNORESTARTIGNORESTARTIGNORE |
138 |
WRITE(20,*) ' a2TempZ',a2TempZ |
WRITE(20,*) ' a2TempZ',a2TempZ |
139 |
C********* GMGS ******************* |
C********* GMGS ******************* |
140 |
C Ratio of time step |
C Ratio of time step |
141 |
asyncfac = 1.0 |
asyncfac = 30.0 |
142 |
C Scale dependent mixing scheme parameters |
C Scale dependent mixing scheme parameters |
143 |
difref = a2tempXY |
difref = a2tempXY |
144 |
difmin = 0.1 * a2tempXY |
difmin = 0.1 * a2tempXY |
405 |
|
|
406 |
CBEGINIGNOREBEGINIGNOREBEGINIGNOREBEGINIGNOREBEGINIGNOREBEGINIGNOREBEGINIGNORE |
CBEGINIGNOREBEGINIGNOREBEGINIGNOREBEGINIGNOREBEGINIGNOREBEGINIGNOREBEGINIGNORE |
407 |
sBeta = 0. |
sBeta = 0. |
408 |
|
sBeta = 7.4E-4 |
409 |
WRITE(20,*) ' sBeta',sBeta |
WRITE(20,*) ' sBeta',sBeta |
410 |
tAlpha = 2.D-4 |
tAlpha = 2.D-4 |
411 |
C tAlpha = 0.D-4 |
C tAlpha = 0.D-4 |
414 |
|
|
415 |
C Numerics parameters |
C Numerics parameters |
416 |
abEps = 0.1 |
abEps = 0.1 |
417 |
rAja = 9.930555555555556D-1 |
rAja = 9.887999800888535E-01 |
418 |
C rAja = 0. |
C rAja = 0. |
419 |
dUnit = 10 |
dUnit = 10 |
420 |
scrUnit1 = 11 |
scrUnit1 = 11 |
464 |
READ (11) tmpXY |
READ (11) tmpXY |
465 |
C Fu = tmpXY/DELPS(1)*G*0.1*uMask(_I3(1,:,:)) ! Dynes/cm**2 -> m/s**2 |
C Fu = tmpXY/DELPS(1)*G*0.1*uMask(_I3(1,:,:)) ! Dynes/cm**2 -> m/s**2 |
466 |
Fu = tmpXY/DELPS(1)*G*0.1 |
Fu = tmpXY/DELPS(1)*G*0.1 |
467 |
|
Fu = 0.1/DELPS(1)*G |
468 |
CLOSE(11) |
CLOSE(11) |
469 |
|
|
470 |
OPEN ( 11, FILE='windy', STATUS='old' , FORM='unformatted') |
OPEN ( 11, FILE='windy', STATUS='old' , FORM='unformatted') |
471 |
READ (11) tmpXY |
READ (11) tmpXY |
472 |
C Fu = tmpXY/DELPS(1)*G*0.1*uMask(_I3(1,:,:)) ! Dynes/cm**2 -> m/s**2 |
C Fu = tmpXY/DELPS(1)*G*0.1*uMask(_I3(1,:,:)) ! Dynes/cm**2 -> m/s**2 |
473 |
Fv = tmpXY/DELPS(1)*G*0.1 |
Fv = tmpXY/DELPS(1)*G*0. |
474 |
CLOSE(11) |
CLOSE(11) |
475 |
|
|
476 |
C |
C |
484 |
CALL PLOT_FIELD(fv,Nx,Ny) |
CALL PLOT_FIELD(fv,Nx,Ny) |
485 |
C STOP |
C STOP |
486 |
CcnhDebugEnds |
CcnhDebugEnds |
487 |
|
|
488 |
|
C T and S profiles |
489 |
|
OPEN ( 11, FILE='LevCli.ptmp.sun.b', STATUS='old' , FORM='unformatted') |
490 |
|
READ (11) tmpXYZ |
491 |
|
Heat = tmpXYZ(_I3(1,:,:)) |
492 |
|
CLOSE(11) |
493 |
|
C T and S profiles |
494 |
|
OPEN ( 11, FILE='LevCli.salt.sun.b', STATUS='old' , FORM='unformatted') |
495 |
|
READ (11) tmpXYZ |
496 |
|
SSurf = tmpXYZ(_I3(1,:,:)) |
497 |
|
CLOSE(11) |
498 |
C |
C |
499 |
END |
END |