/[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.3 by cnh, Fri Jun 12 19:33:32 1998 UTC revision 1.4 by cnh, Mon Jun 15 05:13:54 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)        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.         |
# Line 96  C Line 97  C
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
# Line 110  CSTARTIGNORESTARTIGNORESTARTIGNORESTARTI Line 111  CSTARTIGNORESTARTIGNORESTARTIGNORESTARTI
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
# Line 130  CSTARTIGNORESTARTIGNORESTARTIGNORESTARTI Line 138  CSTARTIGNORESTARTIGNORESTARTIGNORESTARTI
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
# Line 397  C end NOBETACODE Line 405  C end NOBETACODE
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
# Line 405  C      tAlpha     = 0.D-4 Line 414  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
# Line 455  C       fu(i,j) = -tauMax Line 464  C       fu(i,j) = -tauMax
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
# Line 474  Cdbg  WRITE(0,*) ' tauMax = ', tauMax Line 484  Cdbg  WRITE(0,*) ' tauMax = ', tauMax
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
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22