dgoldberg/streamice/streamice_readparms.F

C $Header: /u/gcmpack/MITgcm_contrib/dgoldberg/streamice/streamice_readparms.F,v 1.4 2012/07/27 21:07:13 heimbach Exp $
C $Name:  $


C this needs changes

#include "STREAMICE_OPTIONS.h"

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
CBOP 0
      SUBROUTINE STREAMICE_READPARMS( myThid )

C     !DESCRIPTION:
C     Initialize STREAMICE variables and constants.

C     !USES:
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "STREAMICE.h"
#include "STREAMICE_BDRY.h"

C     !INPUT PARAMETERS:
      INTEGER myThid
CEOP

#ifdef ALLOW_STREAMICE

C     !LOCAL VARIABLES:
C     msgBuf     :: Informational/error message buffer
C     iUnit      :: Work variable for IO unit number
      CHARACTER*(MAX_LEN_MBUF) msgBuf
      INTEGER iUnit

      NAMELIST /STREAMICE_PARM01/
     &     streamice_density, streamice_density_ocean_avg,
     &     A_glen_isothermal, n_glen, eps_glen_min,
     &     C_basal_fric_const, n_basal_friction, 
     &     streamice_vel_update,streamice_cg_tol,streamice_nonlin_tol,
     &     streamice_nonlin_tol_fp,
     &     streamice_max_cg_iter, streamice_max_nl_iter,
     &     STREAMICE_GL_regularize,
     &     STREAMICEthickInit,
     &     STREAMICEthickFile,
     &     STREAMICEcalveMaskFile,
     &     STREAMICEison,
     &     STREAMICE_dump_mdsio, STREAMICE_tave_mdsio,
     &     STREAMICE_dump_mnc, STREAMICE_tave_mnc,
     &     STREAMICE_GL_regularize, STREAMICE_move_front,
     &     STREAMICE_calve_to_mask,
!     &     STREAMICE_construct_matrix,
     &     STREAMICE_lower_cg_tol,
     &     streamice_CFL_factor,
     &     streamice_adjDump

      NAMELIST /STREAMICE_PARM02/
     &     shelf_max_draft,
     &     shelf_min_draft,
     &     shelf_edge_pos,
     &     shelf_slope_scale,
     &     shelf_flat_width,
     &     flow_dir

      NAMELIST /STREAMICE_PARM03/
     &     min_x_noflow_NORTH, max_x_noflow_NORTH,
     &     min_x_noflow_SOUTH, max_x_noflow_SOUTH,
     &     min_y_noflow_WEST, max_y_noflow_WEST,
     &     min_y_noflow_EAST, max_y_noflow_EAST,
     &     min_x_noStress_NORTH, max_x_noStress_NORTH,
     &     min_x_noStress_SOUTH, max_x_noStress_SOUTH,
     &     min_y_noStress_WEST, max_y_noStress_WEST,
     &     min_y_noStress_EAST, max_y_noStress_EAST,
     &     min_x_FluxBdry_NORTH, max_x_FluxBdry_NORTH,
     &     min_x_FluxBdry_SOUTH, max_x_FluxBdry_SOUTH,
     &     min_y_FluxBdry_WEST, max_y_FluxBdry_WEST,
     &     min_y_FluxBdry_EAST, max_y_FluxBdry_EAST,
     &     min_x_Dirich_NORTH, max_x_Dirich_NORTH,
     &     min_x_Dirich_SOUTH, max_x_Dirich_SOUTH,
     &     min_y_Dirich_WEST, max_y_Dirich_WEST,
     &     min_y_Dirich_EAST, max_y_Dirich_EAST,
     &     min_x_CFBC_NORTH, max_x_CFBC_NORTH,
     &     min_x_CFBC_SOUTH, max_x_CFBC_SOUTH,
     &     min_y_CFBC_WEST, max_y_CFBC_WEST,
     &     min_y_CFBC_EAST, max_y_CFBC_EAST,
     &     flux_bdry_val_SOUTH, flux_bdry_val_NORTH,
     &     flux_bdry_val_WEST, flux_bdry_val_EAST

      _BEGIN_MASTER(myThid)

C--   Default values for STREAMICE

      streamice_density = 917. 
      streamice_density_ocean_avg = 1024.
      A_glen_isothermal = 9.461e-18   !  Pa (-1/3) a
      n_glen = 3.
      eps_glen_min = 1.0e-12
      C_basal_fric_const = 31.71 ! Pa (m/a)-1n
      n_basal_friction = 1.
      streamice_vel_update = 169200. ! seconds
      streamice_cg_tol = 1e-6
      streamice_nonlin_tol = 1e-6
      streamice_nonlin_tol_fp = 1.e-14
      streamice_max_cg_iter = 2000
      streamice_max_nl_iter = 100
      streamice_n_sub_regularize = 4
      streamice_CFL_factor = .5
      streamice_adjDump = 0.
      
      STREAMICEthickInit = 'FILE'
      STREAMICEthickFile = ' '
      STREAMICEcalveMaskFile = ' '

      STREAMICEison = .TRUE.
      STREAMICE_tave_mdsio = .TRUE.
      STREAMICE_dump_mdsio = .TRUE.
      STREAMICE_dump_mnc = .FALSE. 
      STREAMICE_tave_mnc = .FALSE.
      STREAMICE_GL_regularize = .FALSE.
      STREAMICE_move_front = .FALSE.
      STREAMICE_calve_to_mask = .FALSE.
!      STREAMICE_construct_matrix = .TRUE.
      STREAMICE_lower_cg_tol = .FALSE.

      min_x_noflow_NORTH = 0.
      max_x_noflow_NORTH = 0.
      min_x_noflow_SOUTH = 0.
      max_x_noflow_SOUTH = 0.
      min_y_noflow_WEST = 0.
      max_y_noflow_WEST = 0.
      min_y_noflow_EAST = 0.
      max_y_noflow_EAST = 0.
      min_x_noStress_NORTH = 0.
      max_x_noStress_NORTH = 0.
      min_x_noStress_SOUTH = 0.
      max_x_noStress_SOUTH = 0.
      min_y_noStress_WEST = 0.
      max_y_noStress_WEST = 0.
      min_y_noStress_EAST = 0.
      max_y_noStress_EAST = 0.
      min_x_FluxBdry_NORTH = 0.
      max_x_FluxBdry_NORTH = 0.
      min_x_FluxBdry_SOUTH = 0.
      max_x_FluxBdry_SOUTH = 0.
      min_y_FluxBdry_WEST = 0.
      max_y_FluxBdry_WEST = 0.
      min_y_FluxBdry_EAST = 0.
      max_y_FluxBdry_EAST = 0.
      min_x_Dirich_NORTH = 0.
      max_x_Dirich_NORTH = 0.
      min_x_Dirich_SOUTH = 0.
      max_x_Dirich_SOUTH = 0.
      min_y_Dirich_WEST = 0.
      max_y_Dirich_WEST = 0.
      min_y_Dirich_EAST = 0.
      max_y_Dirich_EAST = 0.
      min_y_CFBC_WEST = 0.
      max_y_CFBC_WEST = 0.
      min_y_CFBC_EAST = 0.
      max_y_CFBC_EAST = 0.
      flux_bdry_val_SOUTH = 0.
      flux_bdry_val_NORTH = 0.
      flux_bdry_val_WEST = 0.
      flux_bdry_val_EAST = 0.


      WRITE(msgBuf,'(A)') 'STREAMICE_READPARMS: opening data.streamice'
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &                    SQUEEZE_RIGHT , 1)
      CALL OPEN_COPY_DATA_FILE(
     I                     'data.streamice', 'STREAMICE_READPARMS',
     O                     iUnit,
     I                     myThid )

C     Read parameters from open data file
      READ(UNIT=iUnit,NML=STREAMICE_PARM01)
      WRITE(msgBuf,'(A)')
     &    'STREAMICE_READPARMS: read first param block'
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &    SQUEEZE_RIGHT , 1)

      IF (TRIM(STREAMICEthickInit) .eq. "PARAM") THEN
       READ(UNIT=iUnit,NML=STREAMICE_PARM02)
        WRITE(msgBuf,'(A)')
     &    'STREAMICE_READPARMS: read second param block'
       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &    SQUEEZE_RIGHT , 1)
      ENDIF

      READ(UNIT=iUnit,NML=STREAMICE_PARM03)
      WRITE(msgBuf,'(A)')
     &    'STREAMICE_READPARMS: read third param block'
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &    SQUEEZE_RIGHT , 1)
C     Close the open data file
      CLOSE(iUnit)

      streamice_nstep_velocity = NINT (streamice_vel_update / deltaT)

C-    Set Output type flags :
      
#ifdef ALLOW_MNC
      IF (useMNC) THEN
        IF ( .NOT.outputTypesInclusive
     &       .AND. STREAMICE_tave_mnc ) STREAMICE_tave_mdsio = .FALSE.
        IF ( .NOT.outputTypesInclusive
     &       .AND. STREAMICE_dump_mnc ) STREAMICE_dump_mdsio = .FALSE.
      ENDIF
#endif

      _END_MASTER(myThid)

C--   Everyone else must wait for the parameters to be loaded
      _BARRIER

#endif /* ALLOW_STREAMICE */

      RETURN
      END
1	dgoldberg	1.5	C $Header: /u/gcmpack/MITgcm_contrib/dgoldberg/streamice/streamice_readparms.F,v 1.4 2012/07/27 21:07:13 heimbach Exp $
2	heimbach	1.1	C $Name: $
3
4
5			C this needs changes
6
7			#include "STREAMICE_OPTIONS.h"
8
9			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
10			CBOP 0
11			SUBROUTINE STREAMICE_READPARMS( myThid )
12
13			C !DESCRIPTION:
14			C Initialize STREAMICE variables and constants.
15
16			C !USES:
17			IMPLICIT NONE
18			#include "SIZE.h"
19			#include "EEPARAMS.h"
20			#include "PARAMS.h"
21			#include "STREAMICE.h"
22			#include "STREAMICE_BDRY.h"
23
24			C !INPUT PARAMETERS:
25			INTEGER myThid
26			CEOP
27
28			#ifdef ALLOW_STREAMICE
29
30			C !LOCAL VARIABLES:
31			C msgBuf :: Informational/error message buffer
32			C iUnit :: Work variable for IO unit number
33			CHARACTER*(MAX_LEN_MBUF) msgBuf
34			INTEGER iUnit
35
36			NAMELIST /STREAMICE_PARM01/
37			& streamice_density, streamice_density_ocean_avg,
38			& A_glen_isothermal, n_glen, eps_glen_min,
39			& C_basal_fric_const, n_basal_friction,
40			& streamice_vel_update,streamice_cg_tol,streamice_nonlin_tol,
41	dgoldberg	1.5	& streamice_nonlin_tol_fp,
42	heimbach	1.1	& streamice_max_cg_iter, streamice_max_nl_iter,
43			& STREAMICE_GL_regularize,
44			& STREAMICEthickInit,
45			& STREAMICEthickFile,
46			& STREAMICEcalveMaskFile,
47			& STREAMICEison,
48			& STREAMICE_dump_mdsio, STREAMICE_tave_mdsio,
49			& STREAMICE_dump_mnc, STREAMICE_tave_mnc,
50			& STREAMICE_GL_regularize, STREAMICE_move_front,
51			& STREAMICE_calve_to_mask,
52	dgoldberg	1.2	! & STREAMICE_construct_matrix,
53	heimbach	1.1	& STREAMICE_lower_cg_tol,
54	heimbach	1.4	& streamice_CFL_factor,
55			& streamice_adjDump
56	heimbach	1.1
57			NAMELIST /STREAMICE_PARM02/
58			& shelf_max_draft,
59			& shelf_min_draft,
60			& shelf_edge_pos,
61			& shelf_slope_scale,
62			& shelf_flat_width,
63			& flow_dir
64
65			NAMELIST /STREAMICE_PARM03/
66			& min_x_noflow_NORTH, max_x_noflow_NORTH,
67			& min_x_noflow_SOUTH, max_x_noflow_SOUTH,
68			& min_y_noflow_WEST, max_y_noflow_WEST,
69			& min_y_noflow_EAST, max_y_noflow_EAST,
70			& min_x_noStress_NORTH, max_x_noStress_NORTH,
71			& min_x_noStress_SOUTH, max_x_noStress_SOUTH,
72			& min_y_noStress_WEST, max_y_noStress_WEST,
73			& min_y_noStress_EAST, max_y_noStress_EAST,
74			& min_x_FluxBdry_NORTH, max_x_FluxBdry_NORTH,
75			& min_x_FluxBdry_SOUTH, max_x_FluxBdry_SOUTH,
76			& min_y_FluxBdry_WEST, max_y_FluxBdry_WEST,
77			& min_y_FluxBdry_EAST, max_y_FluxBdry_EAST,
78			& min_x_Dirich_NORTH, max_x_Dirich_NORTH,
79			& min_x_Dirich_SOUTH, max_x_Dirich_SOUTH,
80			& min_y_Dirich_WEST, max_y_Dirich_WEST,
81			& min_y_Dirich_EAST, max_y_Dirich_EAST,
82			& min_x_CFBC_NORTH, max_x_CFBC_NORTH,
83			& min_x_CFBC_SOUTH, max_x_CFBC_SOUTH,
84			& min_y_CFBC_WEST, max_y_CFBC_WEST,
85			& min_y_CFBC_EAST, max_y_CFBC_EAST,
86			& flux_bdry_val_SOUTH, flux_bdry_val_NORTH,
87			& flux_bdry_val_WEST, flux_bdry_val_EAST
88
89			_BEGIN_MASTER(myThid)
90
91			C-- Default values for STREAMICE
92
93			streamice_density = 917.
94			streamice_density_ocean_avg = 1024.
95			A_glen_isothermal = 9.461e-18 ! Pa (-1/3) a
96			n_glen = 3.
97			eps_glen_min = 1.0e-12
98			C_basal_fric_const = 31.71 ! Pa (m/a)-1n
99			n_basal_friction = 1.
100			streamice_vel_update = 169200. ! seconds
101			streamice_cg_tol = 1e-6
102			streamice_nonlin_tol = 1e-6
103	dgoldberg	1.5	streamice_nonlin_tol_fp = 1.e-14
104	heimbach	1.1	streamice_max_cg_iter = 2000
105			streamice_max_nl_iter = 100
106			streamice_n_sub_regularize = 4
107			streamice_CFL_factor = .5
108	heimbach	1.4	streamice_adjDump = 0.
109	heimbach	1.1
110			STREAMICEthickInit = 'FILE'
111			STREAMICEthickFile = ' '
112			STREAMICEcalveMaskFile = ' '
113
114			STREAMICEison = .TRUE.
115			STREAMICE_tave_mdsio = .TRUE.
116			STREAMICE_dump_mdsio = .TRUE.
117			STREAMICE_dump_mnc = .FALSE.
118			STREAMICE_tave_mnc = .FALSE.
119			STREAMICE_GL_regularize = .FALSE.
120			STREAMICE_move_front = .FALSE.
121			STREAMICE_calve_to_mask = .FALSE.
122	dgoldberg	1.2	! STREAMICE_construct_matrix = .TRUE.
123	heimbach	1.1	STREAMICE_lower_cg_tol = .FALSE.
124
125			min_x_noflow_NORTH = 0.
126			max_x_noflow_NORTH = 0.
127			min_x_noflow_SOUTH = 0.
128			max_x_noflow_SOUTH = 0.
129			min_y_noflow_WEST = 0.
130			max_y_noflow_WEST = 0.
131			min_y_noflow_EAST = 0.
132			max_y_noflow_EAST = 0.
133			min_x_noStress_NORTH = 0.
134			max_x_noStress_NORTH = 0.
135			min_x_noStress_SOUTH = 0.
136			max_x_noStress_SOUTH = 0.
137			min_y_noStress_WEST = 0.
138			max_y_noStress_WEST = 0.
139			min_y_noStress_EAST = 0.
140			max_y_noStress_EAST = 0.
141			min_x_FluxBdry_NORTH = 0.
142			max_x_FluxBdry_NORTH = 0.
143			min_x_FluxBdry_SOUTH = 0.
144			max_x_FluxBdry_SOUTH = 0.
145			min_y_FluxBdry_WEST = 0.
146			max_y_FluxBdry_WEST = 0.
147			min_y_FluxBdry_EAST = 0.
148			max_y_FluxBdry_EAST = 0.
149			min_x_Dirich_NORTH = 0.
150			max_x_Dirich_NORTH = 0.
151			min_x_Dirich_SOUTH = 0.
152			max_x_Dirich_SOUTH = 0.
153			min_y_Dirich_WEST = 0.
154			max_y_Dirich_WEST = 0.
155			min_y_Dirich_EAST = 0.
156			max_y_Dirich_EAST = 0.
157			min_y_CFBC_WEST = 0.
158			max_y_CFBC_WEST = 0.
159			min_y_CFBC_EAST = 0.
160			max_y_CFBC_EAST = 0.
161			flux_bdry_val_SOUTH = 0.
162			flux_bdry_val_NORTH = 0.
163			flux_bdry_val_WEST = 0.
164			flux_bdry_val_EAST = 0.
165
166
167			WRITE(msgBuf,'(A)') 'STREAMICE_READPARMS: opening data.streamice'
168			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
169			& SQUEEZE_RIGHT , 1)
170			CALL OPEN_COPY_DATA_FILE(
171			I 'data.streamice', 'STREAMICE_READPARMS',
172			O iUnit,
173			I myThid )
174
175			C Read parameters from open data file
176			READ(UNIT=iUnit,NML=STREAMICE_PARM01)
177			WRITE(msgBuf,'(A)')
178			& 'STREAMICE_READPARMS: read first param block'
179			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
180			& SQUEEZE_RIGHT , 1)
181
182			IF (TRIM(STREAMICEthickInit) .eq. "PARAM") THEN
183	dgoldberg	1.3	READ(UNIT=iUnit,NML=STREAMICE_PARM02)
184	heimbach	1.1	WRITE(msgBuf,'(A)')
185			& 'STREAMICE_READPARMS: read second param block'
186			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
187			& SQUEEZE_RIGHT , 1)
188			ENDIF
189
190	dgoldberg	1.3	READ(UNIT=iUnit,NML=STREAMICE_PARM03)
191	heimbach	1.1	WRITE(msgBuf,'(A)')
192			& 'STREAMICE_READPARMS: read third param block'
193			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
194			& SQUEEZE_RIGHT , 1)
195			C Close the open data file
196			CLOSE(iUnit)
197
198			streamice_nstep_velocity = NINT (streamice_vel_update / deltaT)
199
200			C- Set Output type flags :
201
202			#ifdef ALLOW_MNC
203			IF (useMNC) THEN
204			IF ( .NOT.outputTypesInclusive
205			& .AND. STREAMICE_tave_mnc ) STREAMICE_tave_mdsio = .FALSE.
206			IF ( .NOT.outputTypesInclusive
207			& .AND. STREAMICE_dump_mnc ) STREAMICE_dump_mdsio = .FALSE.
208			ENDIF
209			#endif
210
211			_END_MASTER(myThid)
212
213			C-- Everyone else must wait for the parameters to be loaded
214			_BARRIER
215
216			#endif /* ALLOW_STREAMICE */
217
218			RETURN
219			END