/[MITgcm]/MITgcm/pkg/my82/my82_readparms.F
ViewVC logotype

Contents of /MITgcm/pkg/my82/my82_readparms.F

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

Revision 1.1 - (show annotations) (download)
Thu Sep 2 09:11:54 2004 UTC (20 years, 4 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint58l_post, checkpoint57t_post, checkpoint57o_post, checkpoint58e_post, checkpoint57v_post, checkpoint58u_post, checkpoint58w_post, checkpoint57m_post, checkpoint55c_post, checkpoint54e_post, checkpoint57s_post, checkpoint57k_post, checkpoint55d_pre, checkpoint57d_post, checkpoint57g_post, checkpoint60, checkpoint61, checkpoint57b_post, checkpoint57c_pre, checkpoint58r_post, checkpoint55j_post, checkpoint56b_post, checkpoint57i_post, checkpoint57y_post, checkpoint57e_post, checkpoint55h_post, checkpoint58n_post, checkpoint58x_post, checkpoint57g_pre, checkpoint55b_post, checkpoint58t_post, checkpoint58h_post, checkpoint56c_post, checkpoint57y_pre, checkpoint55, checkpoint57f_pre, checkpoint57a_post, checkpoint58q_post, checkpoint54f_post, checkpoint59q, checkpoint59p, checkpoint55g_post, checkpoint59r, checkpoint58j_post, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint55f_post, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint59j, checkpoint57r_post, checkpoint59, checkpoint58, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, eckpoint57e_pre, checkpoint57h_done, checkpoint58f_post, checkpoint57x_post, checkpoint57n_post, checkpoint58d_post, checkpoint58c_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint57f_post, checkpoint58a_post, checkpoint58i_post, checkpoint57q_post, checkpoint58g_post, checkpoint58o_post, checkpoint57z_post, checkpoint57c_post, checkpoint58y_post, checkpoint55e_post, checkpoint58k_post, checkpoint58v_post, checkpoint55a_post, checkpoint58s_post, checkpoint61f, checkpoint61g, checkpoint61d, checkpoint61e, checkpoint61b, checkpoint61c, checkpoint58p_post, checkpoint61a, checkpoint57j_post, checkpoint58b_post, checkpoint57h_pre, checkpoint58m_post, checkpoint57l_post, checkpoint57h_post, checkpoint56a_post, checkpoint55d_post 
o add two new packages
  - pp81 (Packanowski and Philander, 1981), Richardson number and
    stratification dependent mixing
  - my82 (Mellor and Yamada, 1982) level 2 turbulence closure scheme
1 C $Header: $
2 C $Name: $
3
4 #include "MY82_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: MY82_READPARMS
8
9 C !INTERFACE: ==========================================================
10 SUBROUTINE MY82_READPARMS( myThid )
11
12 C !DESCRIPTION:
13 C Initialize MY82 parameters, read in data.pp
14
15 C !USES: ===============================================================
16 IMPLICIT NONE
17 #include "SIZE.h"
18 #include "EEPARAMS.h"
19 #include "MY82.h"
20 #include "PARAMS.h"
21
22 C !INPUT PARAMETERS: ===================================================
23 C myThid :: thread number
24 INTEGER myThid
25
26 C !OUTPUT PARAMETERS: ==================================================
27 C none
28
29 #ifdef ALLOW_MY82
30
31 C !LOCAL VARIABLES: ====================================================
32 C iUnit :: unit number for I/O
33 C msgBuf :: message buffer
34 INTEGER iUnit
35 CHARACTER*(MAX_LEN_MBUF) msgBuf
36 CEOP
37
38 NAMELIST /MY_PARM01/
39 & MYviscMax,
40 & MYdiffMax,
41 & MYhblScale,
42 & MYdumpFreq,
43 & MYMixingMaps,
44 & MYwriteState,
45 & MYtaveFreq
46
47 C This routine has been called by the main model so we set our
48 C internal flag to indicate we are in business
49 MYisON=.TRUE.
50
51 C Set defaults values for parameters in MY82.h
52 MYviscMax = 1. _d 0
53 MYdiffMax = 1. _d 0
54 MYhblScale = 0.1
55 RiMax = 0.1950
56 MYdumpFreq = dumpFreq
57 MYtaveFreq = taveFreq
58 MYMixingMaps = .FALSE.
59 MYwriteState = .FALSE.
60
61 C Open and read the data.pp file
62 _BEGIN_MASTER(myThid)
63 WRITE(msgBuf,'(A)') ' MY82_READPARMS: opening data.my82'
64 CALL PRINT_MESSAGE(msgBuf, standardMessageUnit,
65 & SQUEEZE_RIGHT , 1)
66 CALL OPEN_COPY_DATA_FILE(
67 I 'data.my82', 'MY82_READPARMS',
68 O iUnit,
69 I myThid )
70 READ(UNIT=iUnit,NML=MY_PARM01)
71 WRITE(msgBuf,'(A)')
72 & ' MY82_READPARMS: finished reading data.my82'
73 CALL PRINT_MESSAGE(msgBuf, standardMessageUnit,
74 & SQUEEZE_RIGHT , 1)
75
76 C Close the open data file
77 CLOSE(iUnit)
78 _END_MASTER(myThid)
79
80 C Everyone else must wait for the parameters to be loaded
81 _BARRIER
82
83 C Now set-up any remaining parameters that result from the input parameters
84 IF ( MYviscMax .LE. 0. ) THEN
85 WRITE(msgBuf,'(A)') 'MYviscMax must be greater than zero'
86 CALL PRINT_ERROR( msgBuf , 1)
87 STOP 'ABNORMAL END: S/R MY82_READPARMS'
88 ENDIF
89 IF ( MYdiffMax .LE. 0. ) THEN
90 WRITE(msgBuf,'(A)') 'MYdiffMax must be greater than zero'
91 CALL PRINT_ERROR( msgBuf , 1)
92 STOP 'ABNORMAL END: S/R MY82_READPARMS'
93 ENDIF
94 #endif /* ALLOW_MY82 */
95
96 RETURN
97 END
  ViewVC Help
Powered by ViewVC 1.1.22