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

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

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

Revision 1.1 - (hide annotations) (download)
Thu Sep 2 09:11:54 2004 UTC (19 years, 9 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 mlosch 1.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