-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathhycom.F
executable file
·172 lines (171 loc) · 5.84 KB
/
hycom.F
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
#if defined(USE_ESMF)
program hycom
c
c --- ESMF driver for stand-alone HYCOM ocean model
c
use ESMF_Mod
* use mod_hycom, only : OCEAN_SetServices => HYCOM_SetServices
use mod_hycom, only : end_of_run,
& OCEAN_SetServices => HYCOM_SetServices
c
implicit none
c
c --- Local variables
c
c --- Gridded Components
type(ESMF_GridComp) :: oceanGridComp
c
c --- States, Virtual Machines, and Layouts
type(ESMF_VM) :: worldVM
type(ESMF_State) :: oceanImpState, oceanExpState
integer :: petCount, localPet, split
c
c --- Calendars and clocks
type(ESMF_Clock) :: worldClock
type(ESMF_Clock) :: oceanClock
c
c --- Return codes for error checks
integer :: rc
c
c --- Miscellaneous
integer :: i
c
c-------------------------------------------------------------------------------
c Initialize the ESMF Framework
c-------------------------------------------------------------------------------
c
c --- Set default calendar and log type; get world VM
call ESMF_Initialize(defaultCalendar=ESMF_CAL_GREGORIAN,
& defaultLogType=ESMF_LOG_SINGLE,
& vm=worldVM, rc=rc)
if (rc .ne. ESMF_SUCCESS) stop 99
c
c --- Get VM info
call ESMF_VMGet(worldVM, petCount=petCount, localPET=localPet,
& rc=rc)
if (ESMF_LogMsgFoundError(rc, "ESMF_VMGet failed", rc))
& goto 10
c
c-------------------------------------------------------------------------------
c --- Create section
c-------------------------------------------------------------------------------
c
c --- Create the OCEAN gridded component
oceanGridComp = ESMF_GridCompCreate(vm=worldVM,
& name="OCEAN Gridded Component",
& gridCompType=ESMF_OCEAN,
& rc=rc)
if (ESMF_LogMsgFoundError(rc, "OCEAN GridCompCreate failed", rc))
& goto 10
c
c --- Create empty OCEAN import/export states
oceanImpState = ESMF_StateCreate(stateName="OCEAN Import State",
& stateType=ESMF_STATE_IMPORT,
& rc=rc)
if (ESMF_LogMsgFoundError(rc, "OCEAN ImpState Create failed", rc))
& goto 10
oceanExpState = ESMF_StateCreate(stateName="OCEAN Export State",
& stateType=ESMF_STATE_EXPORT,
& rc=rc)
if (ESMF_LogMsgFoundError(rc, "OCEAN ExpState Create failed", rc))
& goto 10
c
c-------------------------------------------------------------------------------
c --- Register section
c-------------------------------------------------------------------------------
c
c --- Register the OCEAN gridded component
call ESMF_GridCompSetServices(oceanGridComp,
& OCEAN_SetServices, rc)
if (ESMF_LogMsgFoundError(rc, "OCEAN Registration failed", rc))
& goto 10
c
c-------------------------------------------------------------------------------
c --- Initalize Section
c-------------------------------------------------------------------------------
c
c --- Initialize OCEAN gridded component
call ESMF_GridCompInitialize(gridComp=oceanGridComp,
& importState=oceanImpState,
& exportState=oceanExpState,
& clock=worldClock,
& phase=ESMF_SINGLEPHASE,
& blockingflag=ESMF_NONBLOCKING,
& rc=rc)
if (ESMF_LogMsgFoundError(rc, "OCEAN Initialize failed", rc))
& goto 10
*c
*c --- Get copy of OCEAN clock
* call ESMF_GridCompGet(oceanGridComp, clock=oceanClock, rc=rc)
*c
*c --- Initialize WORLD clock using OCEAN clock
* worldClock = ESMF_ClockCreate(clock=oceanClock, rc=rc)
c
c-------------------------------------------------------------------------------
c --- Run Section
c-------------------------------------------------------------------------------
c
do !until end of run
call ESMF_GridCompRun(gridComp=oceanGridComp,
& importState=oceanImpState,
& exportState=oceanExpState,
& clock=worldClock,
& phase=ESMF_SINGLEPHASE,
& blockingflag=ESMF_NONBLOCKING,
& rc=rc)
c
c --- use end_of_run, rather than a ESMF Clock
if (end_of_run) then
exit
endif
enddo
c
c-------------------------------------------------------------------------------
c Finalize Section
c-------------------------------------------------------------------------------
c
c --- Finalize OCEAN gridded component
call ESMF_GridCompFinalize(gridComp=oceanGridComp,
& importState=oceanImpState,
& exportState=oceanExpState,
& clock=worldClock,
& phase=ESMF_SINGLEPHASE,
& blockingflag=ESMF_NONBLOCKING,
& rc=rc)
if (ESMF_LogMsgFoundError(rc, "OCEAN Finalize failed", rc))
& goto 10
c
10 continue
call ESMF_VMBarrier(worldVM)
call ESMF_Finalize(rc=rc)
c
stop
end program hycom
#else
program hycom
c
c --- Non-ESMF driver for stand-alone HYCOM ocean model
c
use mod_hycom, only : end_of_run,
& HYCOM_Init,
& HYCOM_Run,
& HYCOM_Final
c
implicit none
c
c --- Initialize HYCOM.
call HYCOM_Init
c --- Run HYCOM.
do !until end of run
call HYCOM_Run
if (end_of_run) then
exit
endif
enddo
c
c --- Finalize HYCOM.
call HYCOM_Final
c
stop
end program hycom
#endif