forked from NOAA-GFDL/SIS2
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathice_spec.F90
166 lines (136 loc) · 6.04 KB
/
ice_spec.F90
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
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~!
! ice_spec_mod - sea ice and SST specified from data as per GFDL climate group !
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~!
module ice_spec_mod
use fms_mod, only: open_namelist_file, check_nml_error, close_file, &
stdlog, stdout, mpp_pe, mpp_root_pe, write_version_number
use mpp_mod, only: input_nml_file
use mpp_domains_mod, only : domain2d
use time_manager_mod, only: time_type, get_date, set_date
use data_override_mod, only : data_override, data_override_init, data_override_unset_domains
implicit none
include 'netcdf.inc'
private
public :: get_sea_surface
character(len=128), parameter :: version = '$Id: ice_spec.F90,v 1.1.2.1.6.1.2.1.2.1 2013/06/18 22:24:14 nnz Exp $'
character(len=128), parameter :: tagname = '$Name: siena_201305_ice_sis2_5layer_dEdd_nnz $'
logical :: module_is_initialized = .false.
logical :: mcm_ice = .false. ! When mcm_ice=.true., ice is handled as in supersource
real :: sst_pert = 0. ! global temperature perturbation used for sensitivity experiments
real :: minimum_ice_concentration = 0.2
real :: minimum_ice_thickness = 1.0
logical :: do_leads = .true. ! when do_leads=false there is no fractional ice concentration
! also you should set the minimum_ice_concentration = 0.5
logical :: sst_degk = .false. ! when sst_degk=true the input sst data is in degrees Kelvin
! otherwise it is assumed to be in degrees Celsius
!amip date for repeating single day (rsd) option
integer :: repeat_date(3)=(/-1,-1,-1/)
namelist / ice_spec_nml / mcm_ice, do_leads, minimum_ice_concentration, &
minimum_ice_thickness, sst_degk, sst_pert, repeat_date
contains
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~!
! get_sea_surface - get SST, ice concentration and thickness from data !
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~!
subroutine get_sea_surface(Time, ts, cn, iceh, ice_domain, ice_domain_end, ts_in_K)
type (time_type), intent(in) :: Time
real, dimension(:, :), intent(out) :: ts
real, dimension(size(ts,1),size(ts,2),2), intent(out) :: cn
real, dimension(size(ts,1),size(ts,2)), intent(out) :: iceh
type(domain2d), optional, intent(in) :: ice_domain
type(domain2d), optional, intent(in) :: ice_domain_end
logical, optional, intent(in) :: ts_in_K
real, dimension(size(ts,1),size(ts,2)) :: sst, icec
real :: t_sw_freeze0 = -1.8
real :: t_sw_freeze
real, parameter :: T_0degC = 273.15 ! 0 degrees C in Kelvin
logical :: ts_in_degK
integer :: ierr, io, unit
type(time_type) :: Spec_Time
integer :: tod(3), dum1, dum2, dum3
integer :: stdoutunit,stdlogunit
stdoutunit=stdout()
stdlogunit=stdlog()
if (.not.module_is_initialized) then
#ifdef INTERNAL_FILE_NML
read (input_nml_file, nml=ice_spec_nml, iostat=io)
#else
unit = open_namelist_file()
read (unit, ice_spec_nml,iostat=io)
call close_file (unit)
#endif
ierr = check_nml_error(io,'ice_spec_nml')
write (stdoutunit,'(/)')
write (stdoutunit, ice_spec_nml)
write (stdlogunit, ice_spec_nml)
call write_version_number(version, tagname)
module_is_initialized = .true.
endif
if (present(ice_domain)) then
call data_override_unset_domains(unset_Ice=.true., must_be_set=.false.)
call data_override_init(Ice_domain_in = Ice_domain)
endif
ts_in_degK = .true. ; if (present(ts_in_K)) ts_in_degK = ts_in_K
! modify time repeating single day option
if (all(repeat_date>0)) then
call get_date(Time,dum1,dum2,dum3,tod(1),tod(2),tod(3))
Spec_Time = set_date(repeat_date(1),repeat_date(2),repeat_date(3),tod(1),tod(2),tod(3))
else
Spec_Time = Time
endif
t_sw_freeze = t_sw_freeze0
if (sst_degk) then
t_sw_freeze = t_sw_freeze0 + T_0degC ! convert sea water freeze point to degK
endif
icec = 0.0; iceh(:,:) = 0.0; sst(:,:) = t_sw_freeze;
call data_override('ICE', 'sic_obs', icec, Spec_Time)
call data_override('ICE', 'sit_obs', iceh, Spec_Time)
call data_override('ICE', 'sst_obs', sst, Spec_Time)
if (present(ice_domain)) then
call data_override_unset_domains(unset_Ice=.true.)
! Reset the data_override ice domain back to the one expected by the coupler.
if (present(ice_domain_end)) call data_override_init(Ice_domain_in = Ice_domain_end)
endif
if (mcm_ice) then
icec = 0.0
! TK Mod: Limit minimum non-zero sea ice thickness to 0.01m.
! This is to eliminate some very thin but non-zero
! sea ice thickness values, where they really should be zero
! but have become nonzero due to spatial interpolation
! where the input grid and model grid are not
! EXACTLY the same. 0.01 was obtained by trial and
! error to roughly match supersource behavior.
! 5/22/01; 8/23/01
where (iceh < 0.01) iceh=0.0
where (iceh>0.0)
icec = 1.0
sst = t_sw_freeze
end where
else
where (icec >= minimum_ice_concentration)
iceh = max(iceh, minimum_ice_thickness)
sst = t_sw_freeze
elsewhere
icec = 0.0
iceh = 0.0
end where
if (.not.do_leads) then
where (icec >= minimum_ice_concentration) icec = 1.0
endif
endif
where (icec==0.0 .and. sst<=t_sw_freeze) sst = t_sw_freeze+1e-10
! add on non-zero sea surface temperature perturbation (namelist option)
! this perturbation may be useful in accessing model sensitivities
if ( abs(sst_pert) > 0.0001 ) then
sst(:,:) = sst(:,:) + sst_pert
endif
cn(:,:,2) = icec(:,:)
cn(:,:,1) = 1.0 - cn(:,:,2)
if (sst_degk .eqv. ts_in_degK) then ! ts and sst have the same units.
ts(:,:) = sst(:,:)
elseif (ts_in_degK) then ! ts is in K and sst in C.
ts(:,:) = sst(:,:) + T_0degC
else ! ts is in C and sst in K.
ts(:,:) = sst(:,:) - T_0degC
endif
end subroutine get_sea_surface
end module ice_spec_mod