31 |
C iUnit :: Work variable for IO unit number |
C iUnit :: Work variable for IO unit number |
32 |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
33 |
INTEGER iUnit |
INTEGER iUnit |
34 |
|
INTEGER k |
35 |
|
|
36 |
NAMELIST /ATM_PHYS_PARM01/ |
NAMELIST /ATM_PHYS_PARM01/ |
|
c & myPa_MNC, |
|
|
c & myPa_StaV_Cgrid, myPa_Tend_Cgrid, |
|
37 |
& atmPhys_addTendT, atmPhys_addTendS, |
& atmPhys_addTendT, atmPhys_addTendS, |
38 |
& atmPhys_addTendU, atmPhys_addTendV, |
& atmPhys_addTendU, atmPhys_addTendV, |
39 |
|
& atmPhys_tauDampUV, atmPhys_dampUVfac, |
40 |
& atmPhys_stepSST, |
& atmPhys_stepSST, |
|
c & myPa_index1, myPa_index2, myPa_doSwitch2, |
|
|
c & myPa_param1, myPa_param2, |
|
|
c & myPa_string1, myPa_string2, |
|
|
c & myPa_Scal1File, myPa_Scal2File, |
|
|
c & myPa_VelUFile, myPa_VelVFile, |
|
41 |
& atmPhys_SSTFile, atmPhys_QflxFile |
& atmPhys_SSTFile, atmPhys_QflxFile |
42 |
|
|
43 |
_BEGIN_MASTER(myThid) |
_BEGIN_MASTER(myThid) |
44 |
|
|
45 |
C-- Default values for ATM_PHYS |
C-- Default values for ATM_PHYS |
|
c myPa_MNC = useMNC |
|
|
c myPa_StaV_Cgrid = .TRUE. |
|
|
c myPa_Tend_Cgrid = .TRUE. |
|
46 |
atmPhys_addTendT = .TRUE. |
atmPhys_addTendT = .TRUE. |
47 |
atmPhys_addTendS = .TRUE. |
atmPhys_addTendS = .TRUE. |
48 |
atmPhys_addTendU = .TRUE. |
atmPhys_addTendU = .TRUE. |
49 |
atmPhys_addTendV = .TRUE. |
atmPhys_addTendV = .TRUE. |
50 |
|
atmPhys_tauDampUV = 0. _d 0 |
51 |
|
DO k=1,Nr |
52 |
|
atmPhys_dampUVfac(k) = 1. _d 0 |
53 |
|
ENDDO |
54 |
C- additional parameters: |
C- additional parameters: |
55 |
atmPhys_stepSST = .FALSE. |
atmPhys_stepSST = .FALSE. |
|
c myPa_doSwitch2 = .FALSE. |
|
|
c myPa_index1 = 0 |
|
|
c myPa_index2 = 0 |
|
|
c myPa_param1 = 0. _d 0 |
|
|
c myPa_param2 = 0. _d 0 |
|
|
c myPa_string1 = ' ' |
|
|
c myPa_string2 = ' ' |
|
|
C- file names for initial conditions: |
|
|
c myPa_Scal1File = ' ' |
|
|
c myPa_Scal2File = ' ' |
|
|
c myPa_VelUFile = ' ' |
|
|
c myPa_VelVFile = ' ' |
|
56 |
atmPhys_SSTFile = ' ' |
atmPhys_SSTFile = ' ' |
57 |
atmPhys_QflxFile = ' ' |
atmPhys_QflxFile = ' ' |
58 |
|
|
73 |
C Close the open data file |
C Close the open data file |
74 |
CLOSE(iUnit) |
CLOSE(iUnit) |
75 |
|
|
76 |
|
IF ( atmPhys_tauDampUV.LE.zeroRL ) THEN |
77 |
|
DO k=1,Nr |
78 |
|
atmPhys_dampUVfac(k) = 0. _d 0 |
79 |
|
ENDDO |
80 |
|
ELSE |
81 |
|
DO k=1,Nr |
82 |
|
atmPhys_dampUVfac(k) = atmPhys_dampUVfac(k)/atmPhys_tauDampUV |
83 |
|
ENDDO |
84 |
|
ENDIF |
85 |
|
|
86 |
C Make sure that we locally honor the global MNC on/off flag |
C Make sure that we locally honor the global MNC on/off flag |
87 |
c myPa_MNC = myPa_MNC .AND. useMNC |
c myPa_MNC = myPa_MNC .AND. useMNC |
88 |
#ifndef ALLOW_MNC |
#ifndef ALLOW_MNC |