[4] | 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 | !-----------------------------------------------------------------------------! |
---|