1 | ! |
---|
2 | !---5----|----5----|----5----|----5----|----5----|----5----|----5----|72--5----| |
---|
3 | ! Thermo.com is incorparated by an instruction include in ice_th.f , |
---|
4 | ! fontbc.f and acrlbq.f. It comprises the commons associated to |
---|
5 | ! thermodynamic ice computation |
---|
6 | ! CLEANING REQUIRED |
---|
7 | ! |
---|
8 | ! Correspondance between the variables |
---|
9 | ! qlbqb qlbq |
---|
10 | ! qcmbqb qcmbq |
---|
11 | ! thcmb thcm |
---|
12 | ! fstbqb fstrbq |
---|
13 | ! fltbqb ffltbq |
---|
14 | ! fscbqb fscmbq |
---|
15 | ! fsolgb fsolg |
---|
16 | ! ratbqb ratbqg |
---|
17 | ! psbqb psbq |
---|
18 | ! tabqb tabq |
---|
19 | ! qabqb qabq |
---|
20 | ! vabqb vabq |
---|
21 | ! qfvbqb qfvbq |
---|
22 | ! tsb ts |
---|
23 | ! tfub tfu |
---|
24 | ! hnpbqb zhnpbq |
---|
25 | ! hnbqb hnbq |
---|
26 | ! hgbqb hgbq |
---|
27 | ! albqb albq |
---|
28 | ! qstbqb qstobq |
---|
29 | ! fbbqb fbbq |
---|
30 | ! tbqb tbq |
---|
31 | ! dmgbqb dmgbq |
---|
32 | ! dmnbqb dmnbq |
---|
33 | ! qlbbqb qlbsbq |
---|
34 | ! cldqb cloud |
---|
35 | ! dmgwib dmgwi |
---|
36 | ! npb number of points where computations has to be done |
---|
37 | ! npac correspondance between the points |
---|
38 | ! fratsb firg |
---|
39 | ! fcsb fcsg |
---|
40 | ! fleb fleg |
---|
41 | ! albgb albg 06/08/2001 |
---|
42 | ! sal salic 08/03/2002 |
---|
43 | |
---|
44 | COMMON/combq/qlbqb(nbpt),qcmbqb(nbpt),thcmb(nbpt),fstbqb(nbpt), |
---|
45 | & fltbqb(nbpt),fscbqb(nbpt),fsolgb(nbpt),fsolgb2(nbpt), |
---|
46 | & ratbqb(nbpt),psbqb(nbpt),tabqb(nbpt), |
---|
47 | & qabqb(nbpt),vabqb(nbpt),qfvbqb(nbpt),tsb(nbpt), |
---|
48 | & tfub(nbpt),hnpbqb(nbpt),hnbqb(nbpt),hgbqb(nbpt), |
---|
49 | & albqb(nbpt),qstbqb(nbpt),fbbqb(nbpt), |
---|
50 | & dmgbqb(nbpt),dmnbqb(nbpt),qlbbqb(nbpt),cldqb(nbpt), |
---|
51 | & dmgwib(nbpt),npb(nbpt),npac(nbpt),tfs(nbpt), |
---|
52 | & albgb(nbpt),focea(nbpt),fsup(nbpt), qsfcb(nbpt), |
---|
53 | & tdewb(nbpt) |
---|
54 | |
---|
55 | ! heat diffusion equation |
---|
56 | common/diff/ht_s_b(nbpt),ht_i_b(nbpt),ht_w_b(nbpt),t_su_b(nbpt), |
---|
57 | & t_bo_b(nbpt),t_i_b(nbpt,maxnlay),t_s_b(nbpt,maxnlay), |
---|
58 | & s_i_b(nbpt,maxnlay),h_i(maxnlay),h_s(maxnlay), |
---|
59 | & sn_i_b(maxnlay), |
---|
60 | & ab(nbpt),s_i_mean, q_i_b(nbpt,0:maxnlay+2), |
---|
61 | & q_sal_i(nbpt,0:maxnlay+2), |
---|
62 | & q_s_b(nbpt,0:maxnlay+2), |
---|
63 | & old_ht_i_b(nbpt), e_i_b(maxnlay) |
---|
64 | |
---|
65 | common/snow/f_s_im(nbpt) |
---|
66 | |
---|
67 | ! vertical grid |
---|
68 | common/vgrid/z_i(0:maxnlay),z_s(0:maxnlay), |
---|
69 | & dh_s_subl(nbpt), dh_s_prec(nbpt), dh_s_font(nbpt), |
---|
70 | & dh_s_snic(nbpt), |
---|
71 | & dh_i_subl(nbpt), dh_i_font(nbpt), dh_i_bott(nbpt), |
---|
72 | & dh_i_snic(nbpt),dh_snowice(nbpt), |
---|
73 | & dh_s_melt(nbpt),dh_i_melt(nbpt), |
---|
74 | & dh_i_surf(nbpt),dh_s_tot(nbpt) |
---|
75 | ! |
---|
76 | common/comdbq/fratsb(nbpt),fcsb(nbpt), |
---|
77 | & fleb(nbpt),dvsbqb(nbpt),dvbbqb(nbpt),dvlbqb(nbpt), |
---|
78 | & dvnbqb(nbpt), fc_su(nbpt), fc_bo_i(nbpt), |
---|
79 | & fatm(nbpt,jpl), fc_i(nbpt,0:maxnlay), |
---|
80 | & fc_s(nbpt,0:maxnlay) |
---|
81 | |
---|
82 | ! for computation of temperatures after accr/abl |
---|
83 | common/ comtem /enthal(0:maxnlay), |
---|
84 | & aaa,bbb,ccc,discrim, |
---|
85 | & enthalsi |
---|
86 | common/ brintrsp /q_br(maxnlay), dq_i_brf(maxnlay) |
---|
87 | common/ bitzlip /ftrice,tmelts,ftroce |
---|
88 | common/ vgrid /himax,hsmax |
---|
89 | common/ conserv /e_i_o,e_s_o,de_i,de_s,isnow,imelt, |
---|
90 | & z_e_i(maxnlay),z_f_i(0:maxnlay), |
---|
91 | & e_t_o,ti_old(maxnlay),si_old(maxnlay), |
---|
92 | & s_i_mmean(maxnlay),h_i_mmean, |
---|
93 | & dq_i(nbpt,jpl), sum_fluxq(nbpt,jpl), |
---|
94 | & cons_error(nbpt,jpl), surf_error(nbpt,jpl), |
---|
95 | & dq_i_layer(nbpt, maxnlay), |
---|
96 | & qt_i_in(nbpt,jpl), qt_s_in(nbpt,jpl), |
---|
97 | & q_i_layer_in(nbpt,maxnlay), |
---|
98 | & qt_i_fin(nbpt,jpl), qt_s_fin(nbpt,jpl), |
---|
99 | & q_i_layer_fin(nbpt,maxnlay), fprec, fsnic |
---|
100 | |
---|
101 | common/ salt / beta_sal, s_i_new, s_i_snic, e_skel, q_summer, |
---|
102 | & diff_br(maxnlay), rayleigh(maxnlay), fsb, fsbp, |
---|
103 | & w_flood, w_flush |
---|
104 | |
---|
105 | ! units is a flux, including radab |
---|
106 | common/radiation/radab_phy_i(maxnlay),radab_s(maxnlay), |
---|
107 | & radab_phy_s(maxnlay),radtr_i(0:maxnlay), |
---|
108 | & radab_alg_i(maxnlay),radtr_s(0:maxnlay), |
---|
109 | & radab_i(maxnlay) |
---|
110 | |
---|
111 | ! par, energy available at the top of each layer |
---|
112 | common/biology/chla_i(maxnlay), par(maxnlay) |
---|
113 | |
---|
114 | ! Remapping |
---|
115 | INTEGER snswi,icsuswi,icboswi,snicswi |
---|
116 | INTEGER snind,icsuind,icboind,snicind,limsum |
---|
117 | |
---|
118 | |
---|
119 | !-- End of file 'thermo.com' |
---|
120 | !-----------------------------------------------------------------------------! |
---|