| 1 | MODULE ice |
|---|
| 2 | !!====================================================================== |
|---|
| 3 | !! *** MODULE ice *** |
|---|
| 4 | !! Sea Ice physics: diagnostics variables of ice defined in memory |
|---|
| 5 | !!===================================================================== |
|---|
| 6 | #if defined key_ice_lim |
|---|
| 7 | !!---------------------------------------------------------------------- |
|---|
| 8 | !! 'key_ice_lim' : LIM sea-ice model |
|---|
| 9 | !!---------------------------------------------------------------------- |
|---|
| 10 | !! |
|---|
| 11 | !! ** Purpose : Define in memory diagnostics Sea Ice variables |
|---|
| 12 | !! |
|---|
| 13 | !! History : |
|---|
| 14 | !! 2.0 ! 03-08 (C. Ethe) F90: Free form and module |
|---|
| 15 | !!---------------------------------------------------------------------- |
|---|
| 16 | !! LIM 2.0, UCL-LODYC-IPSL (2002) |
|---|
| 17 | !!---------------------------------------------------------------------- |
|---|
| 18 | !! * Modules used |
|---|
| 19 | USE par_ice ! LIM sea-ice parameters |
|---|
| 20 | |
|---|
| 21 | IMPLICIT NONE |
|---|
| 22 | PUBLIC |
|---|
| 23 | |
|---|
| 24 | !! * Share Module variables |
|---|
| 25 | INTEGER , PUBLIC :: & !!! ** ice-dynamic namelist (namicedyn) ** |
|---|
| 26 | nlminn = 103 , & ! minimum index for computation of ice drift(NH) |
|---|
| 27 | nlmaxs = 45 , & ! maximum index for computation of ice drift(SH) |
|---|
| 28 | nbiter = 1 , & ! number of sub-time steps for relaxation |
|---|
| 29 | nbitdr = 250 ! maximum number of iterations for relaxation |
|---|
| 30 | |
|---|
| 31 | REAL(wp), PUBLIC :: & !!! ** ice-dynamic namelist (namicedyn) ** |
|---|
| 32 | epsd = 1.0e-20, & ! tolerance parameter for dynamic |
|---|
| 33 | alpha = 0.5 , & ! coefficient for semi-implicit coriolis |
|---|
| 34 | dm = 0.6e+03, & ! diffusion constant for dynamics |
|---|
| 35 | om = 0.5 , & ! relaxation constant |
|---|
| 36 | resl = 5.0e-05, & ! maximum value for the residual of relaxation |
|---|
| 37 | cw = 5.0e-03, & ! drag coefficient for oceanic stress |
|---|
| 38 | angvg = 0.0 , & ! turning angle for oceanic stress |
|---|
| 39 | pstar = 1.0e+04, & ! first bulk-rheology parameter |
|---|
| 40 | c_rhg = 20.0 , & ! second bulk-rhelogy parameter |
|---|
| 41 | etamn = 0.0e+07, & ! minimun value for viscosity |
|---|
| 42 | creepl = 2.0e-08, & ! creep limit |
|---|
| 43 | ecc = 2.0 , & ! eccentricity of the elliptical yield curve |
|---|
| 44 | ahi0 = 350.e0 ! sea-ice hor. eddy diffusivity coeff. (m2/s) |
|---|
| 45 | |
|---|
| 46 | INTEGER , PUBLIC :: & |
|---|
| 47 | nlmaxn = jpjm1 , & ! maximum index for computation of ice drift(NH) |
|---|
| 48 | #if defined key_mpp |
|---|
| 49 | nlmins = 2 ! minimum index for computation of ice drift(SH) |
|---|
| 50 | #else |
|---|
| 51 | nlmins = 3 ! minimum index for computation of ice drift(SH) |
|---|
| 52 | #endif |
|---|
| 53 | |
|---|
| 54 | REAL(wp), PUBLIC :: & |
|---|
| 55 | usecc2 , & ! = 1.0 / ( ecc * ecc ) |
|---|
| 56 | rhoco , & ! = rau0 * cw |
|---|
| 57 | sangvg, cangvg , & ! sin and cos of the turning angle for ocean stress |
|---|
| 58 | pstarh ! pstar / 2.0 |
|---|
| 59 | |
|---|
| 60 | REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: & |
|---|
| 61 | u_oce, v_oce, & ! surface ocean velocity used in ice dynamics |
|---|
| 62 | ahiu , ahiv , & ! hor. diffusivity coeff. at ocean U- and V-points (m2/s) |
|---|
| 63 | pahu , pahv , & ! ice hor. eddy diffusivity coef. at ocean U- and V-points |
|---|
| 64 | hsnm , hicm , & ! mean snow and ice thicknesses |
|---|
| 65 | ust2s ! friction velocity |
|---|
| 66 | |
|---|
| 67 | |
|---|
| 68 | REAL(wp), DIMENSION(jpi,jpj) :: & |
|---|
| 69 | firic , & ! IR flux over the ice (only used for outputs) |
|---|
| 70 | fcsic , & ! Sensible heat flux over the ice (only used for outputs) |
|---|
| 71 | fleic , & ! Latent heat flux over the ice (only used for outputs) |
|---|
| 72 | qlatic , & ! latent flux |
|---|
| 73 | rdvosif, & ! Variation of volume at surface (only used for outputs) |
|---|
| 74 | rdvobif, & ! Variation of ice volume at the bottom ice (only used for outputs) |
|---|
| 75 | fdvolif, & ! Total variation of ice volume (only used for outputs) |
|---|
| 76 | rdvonif, & ! Lateral Variation of ice volume (only used for outputs) |
|---|
| 77 | sist , & ! Sea-Ice Surface Temperature (Kelvin ??? degree ??? I don't know) |
|---|
| 78 | tfu , & ! Melting point temperature of sea water |
|---|
| 79 | hsnif , & ! Snow thickness |
|---|
| 80 | hicif , & ! Ice thickness |
|---|
| 81 | hicifp , & ! Ice production/melting |
|---|
| 82 | frld , & ! Leads fraction = 1-a/totalarea |
|---|
| 83 | phicif , & ! ice thickness at previous time |
|---|
| 84 | pfrld , & ! Leads fraction at previous time |
|---|
| 85 | qstoif , & ! Energy stored in the brine pockets |
|---|
| 86 | fbif , & ! Heat flux at the ice base |
|---|
| 87 | rdmsnif, & ! Variation of snow mass |
|---|
| 88 | rdmicif, & ! Variation of ice mass |
|---|
| 89 | qldif , & ! heat balance of the lead (or of the open ocean) |
|---|
| 90 | qcmif , & ! Energy needed to bring the ocean surface layer until its freezing |
|---|
| 91 | fdtcn , & ! net downward heat flux from the ice to the ocean |
|---|
| 92 | qdtcn , & ! energy from the ice to the ocean |
|---|
| 93 | ! point (at a factor 2) |
|---|
| 94 | thcm , & ! part of the solar energy used in the lead heat budget |
|---|
| 95 | fstric , & ! Solar flux transmitted trough the ice |
|---|
| 96 | ffltbif, & ! Array linked with the max heat contained in brine pockets (?) |
|---|
| 97 | fscmbq , & ! Linked with the solar flux below the ice (?) |
|---|
| 98 | fsbbq , & ! Also linked with the solar flux below the ice (?) |
|---|
| 99 | qfvbq , & ! Array used to store energy in case of toral lateral ablation (?) |
|---|
| 100 | dmgwi ! Variation of the mass of snow ice |
|---|
| 101 | |
|---|
| 102 | REAL(wp), DIMENSION(jpi,jpj) :: & |
|---|
| 103 | albege , & ! Albedo of the snow or ice (only for outputs) |
|---|
| 104 | albecn , & ! Albedo of the ocean (only for outputs) |
|---|
| 105 | tauc , & ! Cloud optical depth |
|---|
| 106 | sdvt ! u*^2/(Stress/density) |
|---|
| 107 | |
|---|
| 108 | |
|---|
| 109 | REAL(wp),DIMENSION(jpi,jpj) :: & |
|---|
| 110 | u_ice, v_ice, & ! two components of the ice velocity (m/s) |
|---|
| 111 | tio_u, tio_v ! two components of the ice-ocean stress (N/m2) |
|---|
| 112 | |
|---|
| 113 | REAL(wp),DIMENSION(jpi,jpj,nsmax) :: & |
|---|
| 114 | scal0 |
|---|
| 115 | |
|---|
| 116 | REAL(wp), DIMENSION(jpi,jpj,nlayersp1) :: & |
|---|
| 117 | tbif ! Temperature inside the ice/snow layer |
|---|
| 118 | |
|---|
| 119 | REAL(wp), DIMENSION(jpi,jpj,0:kmax+1) :: & |
|---|
| 120 | reslum ! Relative absorption of solar radiation in each ocean level |
|---|
| 121 | |
|---|
| 122 | REAL(wp), DIMENSION(jpi,jpj) :: & |
|---|
| 123 | sxice, syice, sxxice, syyice, sxyice, & ! moments for advection |
|---|
| 124 | sxsn, sysn, sxxsn, syysn, sxysn, & |
|---|
| 125 | sxa, sya, sxxa, syya, sxya, & |
|---|
| 126 | sxc0, syc0, sxxc0, syyc0, sxyc0, & |
|---|
| 127 | sxc1, syc1, sxxc1, syyc1, sxyc1, & |
|---|
| 128 | sxc2, syc2, sxxc2, syyc2, sxyc2, & |
|---|
| 129 | sxst, syst, sxxst, syyst, sxyst |
|---|
| 130 | |
|---|
| 131 | #else |
|---|
| 132 | !!---------------------------------------------------------------------- |
|---|
| 133 | !! Default option Empty module NO LIM sea-ice model |
|---|
| 134 | !!---------------------------------------------------------------------- |
|---|
| 135 | #endif |
|---|
| 136 | |
|---|
| 137 | !!====================================================================== |
|---|
| 138 | END MODULE ice |
|---|
