46 character(len=100),
intent(in) :: anfdatname
48 real(dp),
intent(out) :: dxi, deta, z_sl
50 integer(i4b) :: i, j, kc, kt, kr, n
52 character(len=256) :: ch_attr_title, ch_attr_institution, ch_attr_source, &
53 ch_attr_history, ch_attr_references
56 integer(i2b),
dimension(0:IMAX,0:JMAX) :: maske_conv, n_cts_conv
57 real(sp) :: time_conv, dummy_conv, z_sl_conv, &
58 v_tot_conv, a_grounded_conv, a_floating_conv, &
60 xi_conv(0:imax), eta_conv(0:jmax), &
61 sigma_level_c_conv(0:kcmax), sigma_level_t_conv(0:ktmax), &
62 sigma_level_r_conv(0:krmax)
63 real(sp),
dimension(0:IMAX,0:JMAX) :: lambda_conv, phi_conv, &
65 temp_s_conv, as_perp_conv, &
66 zs_conv, zm_conv, zb_conv, zl_conv, h_c_conv, h_t_conv, h_conv, &
67 q_bm_conv, q_tld_conv, &
70 dzs_dtau_conv, dzm_dtau_conv, dzb_dtau_conv, dzl_dtau_conv, &
71 dh_c_dtau_conv, dh_t_dtau_conv, dh_dtau_conv, &
72 vx_b_g_conv, vy_b_g_conv, vz_b_conv, vh_b_conv, &
73 vx_s_g_conv, vy_s_g_conv, vz_s_conv, vh_s_conv, &
74 temp_b_conv, temph_b_conv, &
75 p_b_w_conv, h_w_conv, q_gl_g_conv
76 real(sp),
dimension(0:IMAX,0:JMAX,0:KCMAX) :: vx_c_conv, vy_c_conv, vz_c_conv, &
77 temp_c_conv, age_c_conv
78 real(sp),
dimension(0:IMAX,0:JMAX,0:KTMAX) :: vx_t_conv, vy_t_conv, vz_t_conv, &
79 omega_t_conv, age_t_conv
80 real(sp),
dimension(0:IMAX,0:JMAX,0:KRMAX) :: temp_r_conv
84 open(11, iostat=ios, &
85 file=anfdatpath//
'/'//trim(anfdatname), &
86 status=
'old', form=
'unformatted')
87 if (ios /= 0) stop
' topography3: Error when opening the initial-value file!'
89 read(11) ch_attr_title
90 read(11) ch_attr_institution
91 read(11) ch_attr_source
92 read(11) ch_attr_history
93 read(11) ch_attr_references
99 read(11) sigma_level_c_conv
100 read(11) sigma_level_t_conv
101 read(11) sigma_level_r_conv
107 read(11) as_perp_conv
125 read(11) omega_t_conv
129 read(11) am_perp_conv
134 read(11) dzs_dtau_conv
135 read(11) dzm_dtau_conv
136 read(11) dzb_dtau_conv
137 read(11) dzl_dtau_conv
138 read(11) dh_c_dtau_conv
139 read(11) dh_t_dtau_conv
140 read(11) dh_dtau_conv
141 read(11,end=110) vx_b_g_conv
142 read(11,end=110) vy_b_g_conv
143 read(11,end=110) vz_b_conv
144 read(11,end=110) vh_b_conv
145 read(11,end=110) vx_s_g_conv
146 read(11,end=110) vy_s_g_conv
147 read(11,end=110) vz_s_conv
148 read(11,end=110) vh_s_conv
149 read(11,end=110) temp_b_conv
150 read(11,end=110) temph_b_conv
151 read(11,end=110) p_b_w_conv
152 read(11,end=110) h_w_conv
153 read(11,end=110) q_gl_g_conv
158 write(6, fmt=
'(/1x,a)')
'topography3: End-of-file condition while reading *.erg.'
159 write(6, fmt=
'(1x,a)')
' Some variables will be undefined.'
163 close(11,status=
'keep')
167 z_sl =
real(z_sl_conv,dp)
169 h_r =
real(h_r_conv,dp)
172 xi(i) =
real(xi_conv(i),dp)
176 eta(j) =
real(eta_conv(j),dp)
182 maske(j,i) = maske_conv(i,j)
183 n_cts(j,i) = n_cts_conv(i,j)
184 zs(j,i) =
real(zs_conv(i,j),dp)
185 zm(j,i) =
real(zm_conv(i,j),dp)
186 zb(j,i) =
real(zb_conv(i,j),dp)
187 zl(j,i) =
real(zl_conv(i,j),dp)
188 h_c(j,i) =
real(H_c_conv(i,j),dp)
189 h_t(j,i) =
real(H_t_conv(i,j),dp)
190 q_bm(j,i) =
real(Q_bm_conv(i,j),dp)/year_sec
191 q_tld(j,i) =
real(Q_tld_conv(i,j),dp)/year_sec
192 am_perp(j,i) =
real(am_perp_conv(i,j),dp)/year_sec
193 qx(j,i) =
real(qx_conv(i,j),dp)/year_sec
194 qy(j,i) =
real(qy_conv(i,j),dp)/year_sec
195 dzs_dtau(j,i) =
real(dzs_dtau_conv(i,j),dp)/year_sec
196 dzm_dtau(j,i) =
real(dzm_dtau_conv(i,j),dp)/year_sec
197 dzb_dtau(j,i) =
real(dzb_dtau_conv(i,j),dp)/year_sec
198 dzl_dtau(j,i) =
real(dzl_dtau_conv(i,j),dp)/year_sec
199 dh_c_dtau(j,i) =
real(dH_c_dtau_conv(i,j),dp)/year_sec
200 dh_t_dtau(j,i) =
real(dH_t_dtau_conv(i,j),dp)/year_sec
201 vx_b_g(j,i) =
real(vx_b_g_conv(i,j),dp)/year_sec
202 vy_b_g(j,i) =
real(vy_b_g_conv(i,j),dp)/year_sec
203 vz_b(j,i) =
real(vz_b_conv(i,j),dp)/year_sec
204 vx_s_g(j,i) =
real(vx_s_g_conv(i,j),dp)/year_sec
205 vy_s_g(j,i) =
real(vy_s_g_conv(i,j),dp)/year_sec
206 vz_s(j,i) =
real(vz_s_conv(i,j),dp)/year_sec
207 temp_b(j,i) =
real(temp_b_conv(i,j),dp)
208 temph_b(j,i) =
real(temph_b_conv(i,j),dp)
209 p_b_w(j,i) =
real(p_b_w_conv(i,j),dp)
210 h_w(j,i) =
real(H_w_conv(i,j),dp)
211 q_gl_g(j,i) =
real(q_gl_g_conv(i,j),dp)/year_sec
214 temp_r(kr,j,i) =
real(temp_r_conv(i,j,kr),dp)
218 vx_t(kt,j,i) =
real(vx_t_conv(i,j,kt),dp)/year_sec
219 vy_t(kt,j,i) =
real(vy_t_conv(i,j,kt),dp)/year_sec
220 vz_t(kt,j,i) =
real(vz_t_conv(i,j,kt),dp)/year_sec
221 omega_t(kt,j,i) =
real(omega_t_conv(i,j,kt),dp)
222 age_t(kt,j,i) =
real(age_t_conv(i,j,kt),dp)*year_sec
226 vx_c(kc,j,i) =
real(vx_c_conv(i,j,kc),dp)/year_sec
227 vy_c(kc,j,i) =
real(vy_c_conv(i,j,kc),dp)/year_sec
228 vz_c(kc,j,i) =
real(vz_c_conv(i,j,kc),dp)/year_sec
229 temp_c(kc,j,i) =
real(temp_c_conv(i,j,kc),dp)
230 age_c(kc,j,i) =
real(age_c_conv(i,j,kc),dp)*year_sec
238 open(23, iostat=ios, &
239 file=inpath//
'/grl/'//zl0_file, &
240 recl=16384, status=
'old')
242 if (ios /= 0) stop
' topography3: Error when opening the zl0 file!'
244 do n=1, 6;
read(23, fmt=
'(a)') ch_dummy;
end do
247 read(23, fmt=*) (zl0(j,i), i=0,imax)
250 close(23, status=
'keep')
259 call
geo_coord(phi(j,i), lambda(j,i), xi(i), eta(j))
277 area(j,i) = sq_g11_g(j,i)*sq_g22_g(j,i)*dxi*deta