Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

geticecover.F90 @ 1009

Last change on this file since 1009 was 332, checked in by aslmd, 13 years ago

LMDZ.MARS and MESOSCALE: added evolutive CO2 seasonal caps according to Titus crocus lines

compiles fine. to be tested in mesoscale polar simulations.

27/10/11 == AS

--> Added geticecover.F90 which computes seasonal ice cover given ls, lati(ngrid), long(ngrid)

as proposed by T. Titus from TES observations [fitting functions for crocus line]
... output is icecover(ngrid) which value is 0 [no ice cover] or 1 [ice cover]
... no calculations are done for latitudes between -40 and +40 [ice cover is directly set to 0]

--> In physiq.F, co2ice is set to a dummy high value to simulate a CO2 cap

wherever icecover(ngrid) is 1. This is done at each timestep before newcondens is called.

--> For the moment this is MESOSCALE only, but potentially useful to everyone.

File size: 14.1 KB

Line
1	SUBROUTINE geticecover( ngrid, ls, lontab, lattab, icecover )
2	!!**********************************************!!
3	! A. Spiga (aymeric.spiga@upmc.fr) October 2011 !
4	!!**********************************************!!
5	IMPLICIT NONE
6	!! INPUTS
7	REAL, DIMENSION(ngrid), INTENT(IN) :: lontab,lattab
8	REAL, INTENT(IN) :: ls
9	INTEGER, INTENT(IN) :: ngrid
10	!! OUTPUTS
11	REAL, DIMENSION(ngrid), INTENT(OUT) :: icecover
12	!! LOCAL
13	REAL :: isitice
14	INTEGER :: ig
15	icecover(:) = 0.
16	DO ig=1,ngrid
17	icecover(ig) = isitice( ls, lontab(ig), lattab(ig) )
18	ENDDO
19	END SUBROUTINE geticecover
20
21	REAL FUNCTION isitice( ls, lon, lat )
22	!!**********************************************!!
23	! A. Spiga (aymeric.spiga@upmc.fr) October 2011 !
24	!!**********************************************!!
25	IMPLICIT NONE
26	REAL, INTENT(IN) :: ls,lon,lat
27	REAL :: nplatcrocus, splatcrocus
28	isitice = 0.
29	!! for speedup purposes:
30	!! useless to call functions for lats between -40. and 40.
31	IF (lat .le. -40.) THEN
32	IF ( lat .le. splatcrocus(ls,lon) ) isitice = 1.
33	ENDIF
34	IF (lat .ge. 40.) THEN
35	IF ( lat .ge. nplatcrocus(ls,lon) ) isitice = 1.
36	ENDIF
37	END FUNCTION isitice
38
39	REAL FUNCTION nplatcrocus( ls, lon )
40	!;Purpose: To return the areocentric latitude of the IR cap edge for TES year "1".
41	!;
42	!;Inputs:
43	!; ls: Season in degrees of L_s.
44	!; lo: East Longitude of interest.
45	!;
46	!;Dates of TES data used: 1999 FEB 28 to 2001 JAN 15
47	!;Dates of TES data used: 2001 JAN 15 to 2002 DEC 2
48	!;Dates of TES data used: 2002 DEC 2 to 1980 JAN 1
49	!;
50	!;Written by T.N. Titus Thu Dec 15 11:54:26 2005
51	!;U.S. Geological Survey Astrogeology Team.
52	!;2255 North Gemini Drive
53	!;Flagstaff, AZ 86001 USA
54	!;http://www.mars-ice.org
55	!;
56	!;If using this function for research or publication, please cite:
57	!; T. N. Titus (2005), Mars Polar Cap Edges Tracked over 3 Full Mars Years,
58	!; 36th Annual Lunar and Planetary Science Conference, March 14-18, 2005,
59	!; in League City, Texas, abstract no.1993
60	!!**********************************************************!!
61	! ADAPTATION: A. Spiga (aymeric.spiga@upmc.fr) October 2011 !
62	!!**********************************************************!!
63	IMPLICIT NONE
64	!! INPUT
65	REAL, INTENT(IN) :: ls,lon
66	!! LOCAL
67	REAL, DIMENSION(40) :: gen1, gen2, quan, c, s
68	REAL, DIMENSION(40,2,3) :: tes
69	INTEGER :: year
70	REAL, DIMENSION(3) :: const, line
71	!!*****************************************************************************************
72	!! CONSTANT
73	const(:) = (/ 71.4648,70.4607,71.5339 /)
74	!! GENERIC TABS
75	gen1(:) = (/ &
76	& 0.00000,1.00000,2.00000,3.00000,4.00000, &
77	& 0.00000,1.00000,2.00000,3.00000,4.00000, &
78	& 0.00000,1.00000,2.00000,3.00000,4.00000, &
79	& 0.00000,1.00000,2.00000,3.00000,4.00000, &
80	& 1.00000,2.00000,3.00000,4.00000, &
81	& 1.00000,2.00000,3.00000,4.00000, &
82	& 1.00000,2.00000,3.00000,4.00000, &
83	& 1.00000,2.00000,3.00000,4.00000, &
84	& 1.00000,2.00000,3.00000,4.00000 &
85	/)
86	gen2(:) = (/ &
87	& 1.00000,1.00000,1.00000,1.00000,1.00000, &
88	& 2.00000,2.00000,2.00000,2.00000,2.00000, &
89	& 3.00000,3.00000,3.00000,3.00000,3.00000, &
90	& 4.00000,4.00000,4.00000,4.00000,4.00000, &
91	& 0.00000, 0.00000, 0.00000, 0.00000, &
92	& -1.00000,-1.00000,-1.00000,-1.00000, &
93	& -2.00000,-2.00000,-2.00000,-2.00000, &
94	& -3.00000,-3.00000,-3.00000,-3.00000, &
95	& -4.00000,-4.00000,-4.00000,-4.00000 &
96	/)
97	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! NORTHERN CAP
98	!!! FIRST YEAR TES
99	tes(:,1,1) = (/ &
100	& -1.04928, -0.217329, 0.101174, -0.0374879, 0.00576928, &
101	& -0.467708, 0.213426, -0.322449, 0.248625, -0.0879115, &
102	& -0.0852889,-0.0280091,-0.185435, 0.0160633, 0.0485029, &
103	& -0.0402631, 0.163700, -0.115792, -0.0468377,-0.00756858, &
104	& -11.1382, 0.918513, 0.754372, -0.288945, 0.0435517, &
105	& 0.143201, -0.195657, 0.0378341,-0.640119, 0.429567, &
106	& -0.114784, -0.0829540,-0.179708, 0.0328001, 0.0342086, &
107	& -0.0193232, 0.0279848, 0.180155, -0.154807, 0.0971645 &
108	/)
109	tes(:,2,1) = (/ &
110	& 0.673412, 0.194463, 0.0213440, -0.219141, 0.124593, &
111	& 0.0526366, -0.604908, 0.113683, 0.0372814, -0.0714123, &
112	& -0.00710610,-0.0933308, 0.0241633, 0.0559299, -0.0924497, &
113	& 0.0747056, 0.0400294,-0.127268, -0.0222301, 0.0272872, &
114	& -19.2296, 4.86860, -0.454119, -0.375147, -0.119044, &
115	& 0.209839, -0.271209, 0.367698, 0.171835, 0.350654, &
116	& -0.288149, -0.0152939, 0.0195934, 0.184573, -0.0272326, &
117	& -0.0466995, 0.299228,- 0.0593913, -0.0660705, 0.0282595 &
118	/)
119	!!! SECOND YEAR TES
120	tes(:,1,2) = (/ &
121	& -0.728762,-0.671618,0.154561,0.112470,-0.152911,-0.633290, &
122	& 0.372959,-0.241842,0.100173,0.0558815,-0.0258014,0.0138332, &
123	& -0.113160,0.0965581,-0.0478080,0.0199076,0.0360519,-0.119565, &
124	& -0.0564988,0.00166589,-10.8664,1.34367,0.334286,-0.0889369, &
125	& 0.170191,-0.241918,0.0406942,-0.0312268,-0.542561,0.440734,&
126	& -0.181151,-0.178009,-0.216867,-0.00154585,0.0381691,0.0333462,&
127	& 0.0558955,0.0512524,-0.0389391,0.0705681 &
128	/)
129	tes(:,2,2) = (/ &
130	& 0.133511,0.301999,0.354979,-0.471594,0.151863,0.0575986, &
131	& -0.478402,-0.0168485,0.148232,-0.128787,-0.0914331,-0.0603194, &
132	& 0.0611823,0.0668322,-0.0833629,0.0176089,0.0570124,-0.0641502, &
133	& 0.0196993,-0.0272954,-18.8258,4.08857,-0.252963,-0.0642321, &
134	& -0.463212,0.380846,-0.00712758,0.0374067,0.230256,0.154810, &
135	& -0.235321,0.0838983,0.0737168,0.115381,-0.00787243,-0.0289239, &
136	& 0.169424,-0.00376776,-0.0134735,0.00663328 &
137	/)
138	!!! THIRD YEAR TES
139	tes(:,1,3) = (/ &
140	& -0.170075,-0.706807,-0.0725778,0.266425,-0.166670,-0.408371, &
141	& 0.106833,-0.318539,0.319422,-0.0361432,-0.0366784,-0.0231430, &
142	& -0.129665,0.0478013,0.0283050,0.399130,-0.177498,-0.113099, &
143	& 0.0885522,-0.0503254,-12.2470,1.40947,0.401761,-0.149156, &
144	& -0.256222,-0.0307999,0.173418,-0.112992,-0.628345,0.286791, &
145	& -0.0395388,-0.0634619,-0.263644,0.0727097,0.0679133,0.0350731, &
146	& -0.140744,0.0673651,0.0457041,-0.0454068 &
147	/)
148	tes(:,2,3) = (/ &
149	& 0.447200,-0.251112,0.597279,-0.374734,-0.0290520,-0.0208475, &
150	& -0.682811,0.146231,0.0122089,-0.0917365,0.00708168,-0.143992, &
151	& 0.0724996,0.102962,-0.0703368,-0.0429708,-0.173888,0.189929, &
152	& -0.107449,-0.0505178,-19.3010,4.83420,-1.09834,0.430944, &
153	& -0.643534,0.710448,-0.299251,0.135935,0.134539,0.320505, &
154	& -0.200302,-0.0106718,0.0166824,0.260126,-0.0992902,-0.0182585, &
155	& -0.0627401,0.232410,-0.126134,-0.0100335 &
156	/)
157	!!! MAIN CALCULATIONS
158	quan = (acos(-1.)/180.)lsgen1 + (acos(-1.)/180.)longen2
159	c = cos(quan)
160	s = sin(quan)
161	line(:) = 0.
162	DO year=1,3
163	line(year) = DOT_PRODUCT(c,tes(:,1,year)) - DOT_PRODUCT(s,tes(:,2,year)) + const(year)
164	!!! this is taken into account in isitice
165	!if (line(year) < -90.) line(year) = -90.
166	!if (line(year) > 90.) line(year) = 90.
167	ENDDO
168	!!! THIS IS TEMPORARY, BUT WHY NOT FOR MODELING
169	nplatcrocus = (line(1)+line(2)+line(3))/3.
170	END FUNCTION nplatcrocus
171
172	REAL FUNCTION splatcrocus( ls, lon )
173	!;Purpose: To return the areocentric latitude of the IR cap edge for TES year "1".
174	!;
175	!;Inputs:
176	!; ls: Season in degrees of L_s.
177	!; lo: East Longitude of interest.
178	!;
179	!;Dates of TES data used: 1999 FEB 28 to 2001 JAN 15
180	!;Dates of TES data used: 2001 JAN 15 to 2002 DEC 2
181	!;Dates of TES data used: 2002 DEC 2 to 1980 JAN 1
182	!;
183	!;Written by T.N. Titus Thu Dec 15 11:54:26 2005
184	!;U.S. Geological Survey Astrogeology Team.
185	!;2255 North Gemini Drive
186	!;Flagstaff, AZ 86001 USA
187	!;http://www.mars-ice.org
188	!;
189	!;If using this function for research or publication, please cite:
190	!; T. N. Titus (2005), Mars Polar Cap Edges Tracked over 3 Full Mars Years,
191	!; 36th Annual Lunar and Planetary Science Conference, March 14-18, 2005,
192	!; in League City, Texas, abstract no.1993
193	!!**********************************************************!!
194	! ADAPTATION: A. Spiga (aymeric.spiga@upmc.fr) October 2011 !
195	!!**********************************************************!!
196	IMPLICIT NONE
197	!! INPUT
198	REAL, INTENT(IN) :: ls,lon
199	!! LOCAL
200	REAL, DIMENSION(60) :: gen1, gen2, quan, c, s
201	REAL, DIMENSION(60,2,2) :: tes
202	INTEGER :: year
203	REAL, DIMENSION(2) :: const, line
204	!!*****************************************************************************************
205	!! CONSTANT
206	const(:) = (/ -69.6039, -68.8420 /)
207	!! GENERIC TABS
208	gen1(:) = (/ &
209	& 0.00000,1.00000,2.00000,3.00000,4.00000, &
210	& 5.00000,0.00000,1.00000,2.00000,3.00000, &
211	& 4.00000,5.00000,0.00000,1.00000,2.00000, &
212	& 3.00000,4.00000,5.00000,0.00000,1.00000, &
213	& 2.00000,3.00000,4.00000,5.00000,0.00000, &
214	& 1.00000,2.00000,3.00000,4.00000,5.00000, &
215	& 1.00000,2.00000,3.00000,4.00000,5.00000, &
216	& 1.00000,2.00000,3.00000,4.00000,5.00000, &
217	& 1.00000,2.00000,3.00000,4.00000,5.00000, &
218	& 1.00000,2.00000,3.00000,4.00000,5.00000, &
219	& 1.00000,2.00000,3.00000,4.00000,5.00000, &
220	& 1.00000,2.00000,3.00000,4.00000,5.00000 &
221	/)
222	gen2(:) = (/ &
223	& 1.00000, 1.00000, 1.00000, 1.00000, 1.00000, &
224	& 1.00000, 2.00000, 2.00000, 2.00000, 2.00000, &
225	& 2.00000, 2.00000, 3.00000, 3.00000, 3.00000, &
226	& 3.00000, 3.00000, 3.00000, 4.00000, 4.00000, &
227	& 4.00000, 4.00000, 4.00000, 4.00000, 5.00000, &
228	& 5.00000, 5.00000, 5.00000, 5.00000, 5.00000, &
229	& 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, &
230	& -1.00000,-1.00000,-1.00000,-1.00000,-1.00000, &
231	& -2.00000,-2.00000,-2.00000,-2.00000,-2.00000, &
232	& -3.00000,-3.00000,-3.00000,-3.00000,-3.00000, &
233	& -4.00000,-4.00000,-4.00000,-4.00000,-4.00000, &
234	& -5.00000,-5.00000,-5.00000,-5.00000,-5.00000 &
235	/)
236	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! SOUTHERN CAP
237	!!! FIRST YEAR TES
238	tes(:,1,1) = (/ &
239	& 3.53762, -1.57584, -2.97332, -1.07670, 0.272623, &
240	& 0.311819, 0.842387, 1.11375, 0.388617, -0.575892, &
241	& -0.502105, -0.0357803,-0.689386, 0.743998, 0.796186, &
242	& 0.304092, -0.172464, -0.290923, -0.542742, -0.577187, &
243	& -0.439117, 0.266552, 0.393525, 0.0654230, 0.296123, &
244	& 0.0949434,-0.162267, -0.175361, 0.000162221,0.105111, &
245	& -12.4148, 0.177898, 1.45822, 0.309320, -0.148393, 3.51232, &
246	& 0.355009, -0.782760, -0.427043, -0.178686, -0.429962, &
247	& -0.649747, -0.487551, 0.259459, 0.139186, -0.964653, &
248	& -0.409536, 0.137911, 0.293971, 0.0248885, 0.238855, &
249	& 0.612519, 0.172469, -0.0462983,-0.0129877, 0.00532444,&
250	& -0.189464, -0.139035, -0.0231797, 0.0903127 &
251	/)
252	tes(:,2,1) = (/ &
253	& 3.21178, 3.81330, 0.357977, -1.85466, -1.01477, &
254	& -0.209586, -0.797844, 0.334659, 0.828534, 0.519493, &
255	& -0.0512511,-0.204263, -0.924137, -0.694012, -0.0666508, &
256	& 0.632104, 0.397041, -0.00521332, 0.933057, 0.0136105, &
257	& -0.646719, -0.480330, -0.0122939, 0.214389, 0.0880269, &
258	& 0.0896852, 0.158164, -0.00717643,-0.130884, -0.0425004, &
259	& -21.5000, -5.88472, -0.767056, 0.261383, 0.232822, &
260	& 1.18087, 2.13123, 1.00744, -0.0527458,-0.0417051, &
261	& 1.05940, 0.00149598,-0.501423, -0.290405, -0.0562713, &
262	& 0.335254, -0.480024, -0.604500, 0.0411040, 0.0403603, &
263	& -0.489831, 0.109007, 0.217500, 0.109769, -0.000845153, &
264	& 0.152139, 0.206418, -0.0980599, -0.180992, -0.00685824 &
265	/)
266	!!! SECOND YEAR TES
267	tes(:,1,2) = (/ &
268	& 3.10733, -0.0926926, -1.38769, -0.638410, 0.0598460, &
269	& 0.0212056,-0.574512, -0.687073, -0.193427, 0.161192, &
270	& 0.242137, 0.0178444, -0.427617, 0.892971, 0.650940, &
271	& -0.0372970,-0.311147, -0.123042, -0.221688, -0.467191, &
272	& -0.659134, 0.0246518, 0.257168, 0.124850, 0.174360, &
273	& 0.320583, 0.136864, -0.0606024, -0.0941828, 0.00372431,&
274	& -11.1678, -0.376356, 0.316746, -0.100600, -0.119510, &
275	& 1.35807, -1.53000, -0.867922, 0.301309, 0.258290, &
276	& 0.0813258, 0.601762, 0.268272, 0.0989900, -0.0737015,&
277	& -0.968944, -0.599203, 0.217550, 0.466813, 0.0202996,&
278	& 0.693900, 0.653567, -0.0888332, -0.219429, -0.0298150,&
279	& -0.186738, -0.329258, -0.155475, 0.0241376, 0.0613455 &
280	/)
281	tes(:,2,2) = (/ &
282	& 0.745574, 2.34362, 0.617967, -0.847462, -0.440483, &
283	& -0.135479, 0.546528, -0.0659467, -0.627697, -0.358895, &
284	& 0.0292012, 0.0797757, -0.953241, -0.386098, 0.325419, &
285	& 0.543846, 0.108148, -0.0596165, 1.15557, 0.436356, &
286	& -0.392191, -0.500867, -0.205739, 0.194082, -0.239364, &
287	& -0.136630, 0.135228, 0.213587, 0.0299697, -0.0101247,&
288	& -19.9097, -5.34634, 0.103437, -0.343399, -0.348643, &
289	& 2.20957, 1.02330, -0.722428, -0.574345, 0.149871, &
290	& -0.635328, -0.317596, 0.154507, 0.110691, 0.0265673,&
291	& 0.593630, -0.516898 ,-0.689954, 0.00465230, 0.120493, &
292	& -0.323229, 0.491035, 0.508041, 0.0373988, -0.106504, &
293	& 0.241758, 0.0470212, -0.236231, -0.186565, 0.0105776 &
294	/)
295	!!! MAIN CALCULATIONS
296	quan = (acos(-1.)/180.)lsgen1 + (acos(-1.)/180.)longen2
297	c = cos(quan)
298	s = sin(quan)
299	line(:) = 0.
300	DO year=1,2
301	line(year) = DOT_PRODUCT(c,tes(:,1,year)) - DOT_PRODUCT(s,tes(:,2,year)) + const(year)
302	!!! this is taken into account in isitice
303	!if (line(year) < -90.) line(year) = -90.
304	!if (line(year) > 90.) line(year) = 90.
305	ENDDO
306	!!! THIS IS TEMPORARY, BUT WHY NOT FOR MODELING
307	splatcrocus = (line(1)+line(2))/2.
308	END FUNCTION splatcrocus
309
310	!PROGRAM main
311	!
312	!implicit none
313	!INTEGER, PARAMETER :: ngrid = 12
314	!REAL,DIMENSION(ngrid) :: lontab,lattab,icecover
315	!REAL :: ls,lo,outputs,isitice
316	!REAL :: nplatcrocus,splatcrocus
317	!INTEGER :: i,j,zels
318	!
319	!DO i=0,360,30
320	!print *, 'ls,np,sp ', float(i), nplatcrocus(float(i),0.), splatcrocus(float(i),0.)
321	!ENDDO
322	!
323	!print *, 'isitice(90.,0.,50.)',isitice(90.,0.,50.)
324	!print *, 'isitice(90.,0.,87.)',isitice(90.,0.,87.)
325	!print *, 'isitice(90.,0.,-50.)',isitice(90.,0.,-50.)
326	!print *, 'isitice(90.,0.,-87.)',isitice(90.,0.,-87.)
327	!print *, 'isitice(290.,0.,50.)',isitice(290.,0.,50.)
328	!print *, 'isitice(290.,0.,87.)',isitice(290.,0.,87.)
329	!print *, 'isitice(290.,0.,-50.)',isitice(290.,0.,-50.)
330	!print *, 'isitice(290.,0.,-87.)',isitice(290.,0.,-87.)
331	!print *, 'isitice(120.,0.,90.)',isitice(120.,0.,90.)
332	!
333	!lontab(:) = 0.
334	!lattab = (/ -90., -75., -60., -45., -30., -15., 15., 30., 45., 60., 75., 90. /)
335	!print *,'ls. lat:', lattab
336	!DO zels=0,360,30
337	! CALL geticecover( ngrid, float(zels), lontab, lattab, icecover )
338	! print *, zels, icecover
339	!ENDDO
340	!
341	!END PROGRAM main

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format