1 | !> \file lect-clim-act-anteis_mod.f90 |
---|
2 | !!Module pour la lecture du climat |
---|
3 | !< |
---|
4 | |
---|
5 | !> \namespace lect_clim_act_anteis |
---|
6 | !! Module pour la lecture du climat |
---|
7 | !! \author ... |
---|
8 | !! \date ... |
---|
9 | !! @note Used module |
---|
10 | !! @note - use module3D_phy |
---|
11 | !! @note - use ablation_ann |
---|
12 | !< |
---|
13 | module lect_clim_act_anteis |
---|
14 | |
---|
15 | use module3D_phy |
---|
16 | use ablation_ann ! le module pdd base sur Tann et Tjuly |
---|
17 | |
---|
18 | character(len=80) :: filin |
---|
19 | |
---|
20 | contains |
---|
21 | |
---|
22 | subroutine input_climat_ref() |
---|
23 | |
---|
24 | ! accumulation de Philippe |
---|
25 | filin='accumHUY40km.dat' |
---|
26 | call lect_eis(nx,ny,precip,filin,DIRNAMEINP) |
---|
27 | !====================================== La reponse est 42 =========== |
---|
28 | write(num_rep_42,*) 'fichier accum : ', filin |
---|
29 | |
---|
30 | ! cas particulier de Vostok |
---|
31 | ivo=101 |
---|
32 | jvo=62 |
---|
33 | do j=jvo-1,jvo+1 |
---|
34 | do i=ivo-1,ivo+1 |
---|
35 | precip(i,j)=0.02 ! valeur plus faible a Vostok. |
---|
36 | end do |
---|
37 | end do |
---|
38 | acc(:,:)=precip(:,:) |
---|
39 | |
---|
40 | |
---|
41 | ! temperature en surface : |
---|
42 | ! parametrisation de Fortuin pour la temperature annuelle. |
---|
43 | do j=1,ny |
---|
44 | do i=1,nx |
---|
45 | |
---|
46 | if (s0(i,j).le.200.) then ! shelfs |
---|
47 | tann(i,j)=49.642-0.943*abs(ylat(i,j)) |
---|
48 | else if ((s0(i,j).gt.200.).and.(s0(i,j).lt.1500.)) then ! pente |
---|
49 | tann(i,j)=36.689-0.005102*s0(i,j)-0.725*abs(ylat(i,j)) |
---|
50 | else if (s0(i,j).ge.1500.) then ! plateau |
---|
51 | tann(i,j)=7.405-0.014285*s0(i,j)-0.180*abs(ylat(i,j)) |
---|
52 | endif |
---|
53 | |
---|
54 | ta0(i,j)=tann(i,j) |
---|
55 | ! pour la temperature d'ete, idem parametrisation huybrechts |
---|
56 | tjuly(i,j)=tann(i,j)-17.65+0.00222*s0(i,j)& |
---|
57 | +0.40802*abs(ylat(i,j)) |
---|
58 | end do |
---|
59 | end do |
---|
60 | |
---|
61 | !====================================== La reponse est 42 =========== |
---|
62 | write(num_rep_42,*) '! temperatures annuelles : parametrisation Fortuin' |
---|
63 | write(num_rep_42,*) '! temperature ete : idem parametrisation Huybrechts :' |
---|
64 | write(num_rep_42,*) '! TJULY(I,J)=TANN(I,J)-17.65+0.00222*S(I,J)+0.40802*abs(YLAT(I,J))' |
---|
65 | write(num_rep_42,*) |
---|
66 | |
---|
67 | |
---|
68 | end subroutine input_climat_ref |
---|
69 | |
---|
70 | end module lect_clim_act_anteis |
---|