| 1 | subroutine evap_generic(ngrid,nlayer,nq,dtime,pt,pq,pdq,pdt,igcm_generic_vap, & |
|---|
| 2 | igcm_generic_ice,dqevap,dtevap,qevap,tevap) |
|---|
| 3 | Use genericcommon_h |
|---|
| 4 | Use tracer_h |
|---|
| 5 | implicit none |
|---|
| 6 | !================================================================== |
|---|
| 7 | ! |
|---|
| 8 | ! Purpose |
|---|
| 9 | ! ------- |
|---|
| 10 | ! Evaporation of all generic tracers in the atmopshere. |
|---|
| 11 | ! |
|---|
| 12 | ! Authors |
|---|
| 13 | ! ------- |
|---|
| 14 | ! Adapted from evap.F by Lucas Teinturier (2022) |
|---|
| 15 | ! evap.F is adapted from the LMDTERRE code by B. Charnay (2010) |
|---|
| 16 | ! Original author Z. X. Li (1993) |
|---|
| 17 | ! |
|---|
| 18 | !================================================================== |
|---|
| 19 | !Inputs |
|---|
| 20 | integer, intent(in) :: ngrid,nlayer,nq |
|---|
| 21 | real, intent(in) :: pt(ngrid,nlayer) !Temperature (K) |
|---|
| 22 | real, intent(in) :: pq(ngrid,nlayer,nq) ! Tracers (kg/kg) |
|---|
| 23 | real, intent(in) :: pdt(ngrid,nlayer) !Tendency on temperature (K/s) |
|---|
| 24 | real, intent(in) :: pdq(ngrid,nlayer,nq) ! Tendency on tracer (kg/kg/s) |
|---|
| 25 | real, intent(in) :: dtime ! physical timestep |
|---|
| 26 | integer,intent(in) :: igcm_generic_vap ! index of the vapor generic tracer |
|---|
| 27 | integer,intent(in) :: igcm_generic_ice !index of the ice generic tracer |
|---|
| 28 | |
|---|
| 29 | ! Ouputs |
|---|
| 30 | real,intent(out) :: qevap(ngrid,nlayer,nq) |
|---|
| 31 | real,intent(out) :: tevap(ngrid,nlayer) |
|---|
| 32 | real,intent(out) :: dqevap(ngrid,nlayer) |
|---|
| 33 | real,intent(out) :: dtevap(ngrid,nlayer) |
|---|
| 34 | |
|---|
| 35 | ! Local variables |
|---|
| 36 | real zlvdcp |
|---|
| 37 | integer l,ig ! Used for loops on the layers and grid and the tracers respectively |
|---|
| 38 | ! Evaporate all the ice from the tracer |
|---|
| 39 | zlvdcp = RLVTT/cpp ! RLVTT is the latent heat of vaporization (comes from genericcommon_h attention) |
|---|
| 40 | DO l=1,nlayer |
|---|
| 41 | DO ig=1,ngrid |
|---|
| 42 | qevap(ig,l,igcm_generic_vap) = pq(ig,l,igcm_generic_vap) & |
|---|
| 43 | +pdq(ig,l,igcm_generic_vap)*dtime |
|---|
| 44 | qevap(ig,l,igcm_generic_ice) = pq(ig,l,igcm_generic_ice) & |
|---|
| 45 | +pdq(ig,l,igcm_generic_ice)*dtime |
|---|
| 46 | tevap(ig,l) = pt(ig,l)+pdt(ig,l)*dtime |
|---|
| 47 | dqevap(ig,l) = max(0.,qevap(ig,l,igcm_generic_ice)/dtime) |
|---|
| 48 | dtevap(ig,l) = -dqevap(ig,l)*zlvdcp |
|---|
| 49 | qevap(ig,l,igcm_generic_vap) = qevap(ig,l,igcm_generic_vap) & |
|---|
| 50 | +dqevap(ig,l)*dtime |
|---|
| 51 | qevap(ig,l,igcm_generic_ice) = 0. |
|---|
| 52 | tevap(ig,l) = tevap(ig,l) + dtevap(ig,l)*dtime |
|---|
| 53 | ENDDO ! ngrid |
|---|
| 54 | ENDDO ! nlayer |
|---|
| 55 | |
|---|
| 56 | |
|---|
| 57 | end subroutine evap_generic |
|---|
| 58 | |
|---|
