7 dtco, dts, o, or, sg, s, u, &
8 hprogram,hatmfile,hatmfiletype,hpgdfile,hpgdfiletype)
50 USE modi_get_type_dim_n
51 USE modi_gltools_readnam
56 USE modi_prep_hor_seaflux_field
59 USE modi_gltools_alloc
61 USE yomhook
,ONLY : lhook, dr_hook
62 USE parkind1
,ONLY : jprb
74 TYPE(ocean_t),
INTENT(INOUT) :: o
80 CHARACTER(LEN=6),
INTENT(IN) :: hprogram
81 CHARACTER(LEN=28),
INTENT(IN) :: hatmfile
82 CHARACTER(LEN=6),
INTENT(IN) :: hatmfiletype
83 CHARACTER(LEN=28),
INTENT(IN) :: hpgdfile
84 CHARACTER(LEN=6),
INTENT(IN) :: hpgdfiletype
93 REAL(KIND=JPRB) :: zhook_handle
97 IF (lhook) CALL dr_hook(
'PREP_SEAICE',0,zhook_handle)
108 s%CSEAICE_SCHEME=cprep_seaice_scheme
109 IF ( s%CSEAICE_SCHEME ==
'GELATO' )
THEN
113 s%LHANDLE_SIC = .false.
114 IF(trim(s%CSEAICE_SCHEME)/=
'NONE' .OR. trim(s%CINTERPOL_SIC)/=
'NONE' )
THEN
122 IF (s%LHANDLE_SIC)
THEN
125 hprogram,
'SIC ',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype)
133 s%LINTERPOL_SIC=.false.
134 IF(trim(s%CINTERPOL_SIC)/=
'NONE')
THEN
135 s%LINTERPOL_SIC=.true.
138 IF(trim(s%CINTERPOL_SIT)/=
'NONE')
THEN
139 s%LINTERPOL_SIT=.true.
142 IF(s%LINTERPOL_SIC)
THEN
147 ALLOCATE(s%XSIC_MTH(
SIZE(s%XSIC),inmth))
149 s%XSIC_MTH(:,jmth)=s%XSIC(:)
154 IF(s%LINTERPOL_SIT)
THEN
159 ALLOCATE(s%XSIT_MTH(
SIZE(s%XSIC),inmth))
161 s%XSIT_MTH(:,jmth)=xundef
184 s%TGLT%sit(:,:,1)%age=0.
186 s%TGLT%sit(:,:,1)%vmp=0.
188 s%TGLT%sit(:,:,1)%asn=0.
190 s%TGLT%sit(:,:,1)%fsi=0.
192 s%TGLT%sit(:,:,1)%hsi=1.*s%TGLT%sit(:,:,1)%fsi
194 s%TGLT%sit(:,:,1)%ssi=0.
196 s%TGLT%sit(:,:,1)%tsf=260.
198 s%TGLT%sit(:,:,1)%hsn=0.
200 s%TGLT%sit(:,:,1)%rsn=100.
205 s%TGLT%sil(:,:,:,1)%ent=-1000.
207 IF (lhook) CALL dr_hook(
'PREP_SEAICE',1,zhook_handle)
subroutine get_type_dim_n(DTCO, U, HTYPE, KDIM)
subroutine prep_hor_seaflux_field(DTCO, UG, U, DTS, O, OR, SG, S, HPROGRAM, HSURF, HATMFILE, HATMFILETYPE, HPGDFILE, HPGDFILETYPE)
subroutine prep_seaice(UG, DTCO, DTS, O, OR, SG, S, U, HPROGRAM, HATMFILE, HATMFILETYPE, HPGDFILE, HPGDFILETYPE)
subroutine get_luout(HPROGRAM, KLUOUT)