/[MITgcm]/MITgcm/pkg/streamice/streamice_readparms.F
ViewVC logotype

Diff of /MITgcm/pkg/streamice/streamice_readparms.F

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

revision 1.7 by dgoldberg, Fri Sep 5 14:25:11 2014 UTC revision 1.8 by jmc, Fri Sep 12 22:46:06 2014 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
   
4  C this needs changes  C this needs changes
5    
6  #include "STREAMICE_OPTIONS.h"  #include "STREAMICE_OPTIONS.h"
# Line 21  C     !USES: Line 20  C     !USES:
20  #include "STREAMICE.h"  #include "STREAMICE.h"
21  #include "STREAMICE_BDRY.h"  #include "STREAMICE_BDRY.h"
22  #ifdef ALLOW_STREAMICE_FLUX_CONTROL  #ifdef ALLOW_STREAMICE_FLUX_CONTROL
23  #include "STREAMICE_CTRL_FLUX.h"        #include "STREAMICE_CTRL_FLUX.h"
24  #endif  #endif
25    
26  C     !INPUT PARAMETERS:  C     !INPUT PARAMETERS:
# Line 39  C     iUnit      :: Work variable for IO Line 38  C     iUnit      :: Work variable for IO
38        NAMELIST /STREAMICE_PARM01/        NAMELIST /STREAMICE_PARM01/
39       &     streamice_density, streamice_density_ocean_avg,       &     streamice_density, streamice_density_ocean_avg,
40       &     B_glen_isothermal, n_glen, eps_glen_min, eps_u_min,       &     B_glen_isothermal, n_glen, eps_glen_min, eps_u_min,
41       &     C_basal_fric_const, n_basal_friction,       &     C_basal_fric_const, n_basal_friction,
42       &     streamice_vel_update,streamice_cg_tol,streamice_nonlin_tol,       &     streamice_vel_update,streamice_cg_tol,streamice_nonlin_tol,
43       &     streamice_nonlin_tol_fp,       &     streamice_nonlin_tol_fp,
44       &     streamice_max_cg_iter, streamice_max_nl_iter,       &     streamice_max_cg_iter, streamice_max_nl_iter,
# Line 50  C     iUnit      :: Work variable for IO Line 49  C     iUnit      :: Work variable for IO
49       &     STREAMICEthickFile,       &     STREAMICEthickFile,
50       &     STREAMICEcalveMaskFile,       &     STREAMICEcalveMaskFile,
51       &     STREAMICEcostMaskFile,       &     STREAMICEcostMaskFile,
      &     STREAMICEison,  
52       &     STREAMICE_dump_mdsio, STREAMICE_tave_mdsio,       &     STREAMICE_dump_mdsio, STREAMICE_tave_mdsio,
53       &     STREAMICE_dump_mnc, STREAMICE_tave_mnc,       &     STREAMICE_dump_mnc, STREAMICE_tave_mnc,
54       &     STREAMICE_move_front,       &     STREAMICE_move_front,
# Line 74  C     iUnit      :: Work variable for IO Line 72  C     iUnit      :: Work variable for IO
72       &     STREAMICEuFaceBdryFile,       &     STREAMICEuFaceBdryFile,
73       &     STREAMICEvFaceBdryFile,       &     STREAMICEvFaceBdryFile,
74       &     STREAMICEuDirichValsFile,       &     STREAMICEuDirichValsFile,
75       &     STREAMICEvDirichValsFile,       &     STREAMICEvDirichValsFile,
76       &     STREAMICEuMassFluxFile,       &     STREAMICEuMassFluxFile,
77       &     STREAMICEvMassFluxFile,       &     STREAMICEvMassFluxFile,
78       &     STREAMICEuNormalStressFile,       &     STREAMICEuNormalStressFile,
79       &     STREAMICEvNormalStressFile,       &     STREAMICEvNormalStressFile,
80       &     STREAMICEuShearStressFile,       &     STREAMICEuShearStressFile,
# Line 89  C     iUnit      :: Work variable for IO Line 87  C     iUnit      :: Work variable for IO
87       &     STREAMICE_ppm_driving_stress,       &     STREAMICE_ppm_driving_stress,
88       &     STREAMICE_h_ctrl_const_surf,       &     STREAMICE_h_ctrl_const_surf,
89       &     streamice_wgt_drift,streamice_wgt_surf,streamice_wgt_vel,       &     streamice_wgt_drift,streamice_wgt_surf,streamice_wgt_vel,
90       &     streamice_wgt_avthick, streamice_wgt_tikh,       &     streamice_wgt_avthick, streamice_wgt_tikh,
91       &     streamice_addl_backstress,       &     streamice_addl_backstress,
92       &     streamice_smooth_gl_width,       &     streamice_smooth_gl_width,
93       &     streamice_adot_uniform,       &     streamice_adot_uniform,
# Line 117  C     iUnit      :: Work variable for IO Line 115  C     iUnit      :: Work variable for IO
115       &     streamice_ctrl_flux_scaleVel       &     streamice_ctrl_flux_scaleVel
116  #endif  #endif
117    
   
118        NAMELIST /STREAMICE_PARM02/        NAMELIST /STREAMICE_PARM02/
119       &     shelf_max_draft,       &     shelf_max_draft,
120       &     shelf_min_draft,       &     shelf_min_draft,
# Line 151  C     iUnit      :: Work variable for IO Line 148  C     iUnit      :: Work variable for IO
148       &     flux_bdry_val_WEST, flux_bdry_val_EAST,       &     flux_bdry_val_WEST, flux_bdry_val_EAST,
149       &     STREAMICE_NS_periodic, STREAMICE_EW_periodic       &     STREAMICE_NS_periodic, STREAMICE_EW_periodic
150    
151    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
152    
153          IF ( .NOT.useStreamIce ) THEN
154    C-    pkg STREAMICE is not used
155            _BEGIN_MASTER(myThid)
156    C-    Track pkg activation status:
157             STREAMICEisOn = .FALSE.
158    C     print a (weak) warning if data.streamice is found
159             CALL PACKAGES_UNUSED_MSG( 'useStreamIce', ' ', ' ' )
160            _END_MASTER(myThid)
161            RETURN
162          ENDIF
163    
164        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
165    
166    C This routine has been called by the main model so we set our
167    C internal flag to indicate we are in business
168          STREAMICEisOn = .TRUE.
169    
170  C--   Default values for STREAMICE  C--   Default values for STREAMICE
171    
172        streamice_density = 917.        streamice_density = 917.
173        streamice_density_ocean_avg = 1024.        streamice_density_ocean_avg = 1024.
174        B_glen_isothermal = 9.461e-18   !  Pa (-1/3) a        B_glen_isothermal = 9.461e-18   !  Pa (-1/3) a
175        n_glen = 3.        n_glen = 3.
# Line 180  C--   Default values for STREAMICE Line 194  C--   Default values for STREAMICE
194        streamice_ky_b_init = 1.        streamice_ky_b_init = 1.
195        streamice_wgt_drift = 0.        streamice_wgt_drift = 0.
196        streamice_wgt_tikh = 0.        streamice_wgt_tikh = 0.
197        streamice_wgt_surf = 0.            streamice_wgt_surf = 0.
198        streamice_wgt_vel = 0.            streamice_wgt_vel = 0.
199        streamice_wgt_avthick = 0.        streamice_wgt_avthick = 0.
200        streamice_addl_backstress = 0.0        streamice_addl_backstress = 0.0
201        streamice_smooth_gl_width = 0.0        streamice_smooth_gl_width = 0.0
# Line 243  C--   Default values for STREAMICE Line 257  C--   Default values for STREAMICE
257        ENDDO        ENDDO
258  #endif  #endif
259    
       STREAMICEison = .TRUE.  
260        STREAMICE_tave_mdsio = .TRUE.        STREAMICE_tave_mdsio = .TRUE.
261        STREAMICE_dump_mdsio = .TRUE.        STREAMICE_dump_mdsio = .TRUE.
262        STREAMICE_dump_mnc = .FALSE.        STREAMICE_dump_mnc = .FALSE.
263        STREAMICE_tave_mnc = .FALSE.        STREAMICE_tave_mnc = .FALSE.
264  !      STREAMICE_GL_regularize = .FALSE.  !      STREAMICE_GL_regularize = .FALSE.
265        STREAMICE_move_front = .FALSE.        STREAMICE_move_front = .FALSE.
# Line 304  C--   Default values for STREAMICE Line 317  C--   Default values for STREAMICE
317        flux_bdry_val_EAST = 0.        flux_bdry_val_EAST = 0.
318    
319        STREAMICE_NS_periodic = .FALSE.        STREAMICE_NS_periodic = .FALSE.
320        STREAMICE_EW_periodic = .FALSE.        STREAMICE_EW_periodic = .FALSE.
321    
322        WRITE(msgBuf,'(A)') 'STREAMICE_READPARMS: opening data.streamice'        WRITE(msgBuf,'(A)') 'STREAMICE_READPARMS: opening data.streamice'
323        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
# Line 345  C     Read parameters from open data fil Line 358  C     Read parameters from open data fil
358       &    SQUEEZE_RIGHT , 1)       &    SQUEEZE_RIGHT , 1)
359  #endif  #endif
360    
   
361        READ(UNIT=iUnit,NML=STREAMICE_PARM03)        READ(UNIT=iUnit,NML=STREAMICE_PARM03)
362        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)')
363       &    'STREAMICE_READPARMS: read third param block'       &    'STREAMICE_READPARMS: read third param block'
# Line 353  C     Read parameters from open data fil Line 365  C     Read parameters from open data fil
365       &    SQUEEZE_RIGHT , 1)       &    SQUEEZE_RIGHT , 1)
366        CLOSE(iUnit)        CLOSE(iUnit)
367    
   
368  #ifdef ALLOW_STREAMICE_FLUX_CONTROL  #ifdef ALLOW_STREAMICE_FLUX_CONTROL
369    
370        CALL OPEN_COPY_DATA_FILE(        CALL OPEN_COPY_DATA_FILE(
# Line 361  C     Read parameters from open data fil Line 372  C     Read parameters from open data fil
372       O                     iUnit,       O                     iUnit,
373       I                     myThid )       I                     myThid )
374    
   
375        READ(UNIT=iUnit,NML=STREAMICE_PARMFLUXCTRL)        READ(UNIT=iUnit,NML=STREAMICE_PARMFLUXCTRL)
376        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)')
377       &    'STREAMICE_READPARMS: read flux_ctrl param block'       &    'STREAMICE_READPARMS: read flux_ctrl param block'
# Line 370  C     Read parameters from open data fil Line 380  C     Read parameters from open data fil
380        CLOSE(iUnit)        CLOSE(iUnit)
381  #endif  #endif
382    
   
   
383        streamice_nstep_velocity = NINT (streamice_vel_update / deltaT)        streamice_nstep_velocity = NINT (streamice_vel_update / deltaT)
384    
385  C-    Set Output type flags :  C-    Set Output type flags :
386          
387  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
388        IF (useMNC) THEN        IF (useMNC) THEN
389          IF ( .NOT.outputTypesInclusive          IF ( .NOT.outputTypesInclusive
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22