| 1 | SUBROUTINE LW_venus_ve( |
|---|
| 2 | S PPB, pt, psi, deltapsi, |
|---|
| 3 | S PCOOL, |
|---|
| 4 | S PTOPLW,PSOLLW,PSOLLWDN, |
|---|
| 5 | S ZFLNET) |
|---|
| 6 | |
|---|
| 7 | use dimphy |
|---|
| 8 | use cpdet_mod, only: cpdet |
|---|
| 9 | IMPLICIT none |
|---|
| 10 | |
|---|
| 11 | #include "YOMCST.h" |
|---|
| 12 | C |
|---|
| 13 | C ------------------------------------------------------------------ |
|---|
| 14 | C |
|---|
| 15 | C PURPOSE. |
|---|
| 16 | C -------- |
|---|
| 17 | C |
|---|
| 18 | c This routine uses the NER matrix |
|---|
| 19 | c (computed for a given cell and temp profile in radlwsw) |
|---|
| 20 | c to compute cooling rates and radiative fluxes. |
|---|
| 21 | c |
|---|
| 22 | C AUTHOR. |
|---|
| 23 | C ------- |
|---|
| 24 | C Sebastien Lebonnois |
|---|
| 25 | C |
|---|
| 26 | C MODIFICATIONS. |
|---|
| 27 | C -------------- |
|---|
| 28 | C ORIGINAL : 27/07/2005 |
|---|
| 29 | C version multimatrice (topographie, sommet nuages): 20/12/2006 |
|---|
| 30 | C ------------------------------------------------------------------ |
|---|
| 31 | C |
|---|
| 32 | C* ARGUMENTS: |
|---|
| 33 | C |
|---|
| 34 | c inputs |
|---|
| 35 | |
|---|
| 36 | REAL PPB(klev+1) ! inter-couches PRESSURE (bar) |
|---|
| 37 | REAL pt(klev) ! mid-layer temperature |
|---|
| 38 | real psi(0:klev+1,0:klev+1) ! NER in W/m**2 |
|---|
| 39 | real deltapsi(0:klev+1,0:klev+1) ! D NER / DT in W/m**2/K |
|---|
| 40 | C |
|---|
| 41 | c output |
|---|
| 42 | |
|---|
| 43 | REAL PCOOL(klev) ! LONGWAVE COOLING (K/s) within each layer |
|---|
| 44 | REAL PTOPLW ! LONGWAVE FLUX AT T.O.A. (net, + vers le haut) |
|---|
| 45 | REAL PSOLLW ! LONGWAVE FLUX AT SURFACE (net, + vers le haut) |
|---|
| 46 | REAL PSOLLWDN ! LONGWAVE FLUX AT SURFACE (down, + vers le bas) |
|---|
| 47 | REAL ZFLNET(klev+1) ! net thermal flux at ppb levels (+ vers le haut) |
|---|
| 48 | |
|---|
| 49 | C |
|---|
| 50 | C* LOCAL VARIABLES: |
|---|
| 51 | C |
|---|
| 52 | integer i,j,p |
|---|
| 53 | real zlnet(klev+1) ! net thermal flux (W/m**2) |
|---|
| 54 | real dzlnet(0:klev) ! Radiative budget (W/m**2) |
|---|
| 55 | real pdp(klev) ! epaisseur de la couche en pression (Pa) |
|---|
| 56 | |
|---|
| 57 | c -------------------------- |
|---|
| 58 | c Calculation of the fluxes |
|---|
| 59 | c -------------------------- |
|---|
| 60 | |
|---|
| 61 | c flux aux intercouches: |
|---|
| 62 | c zlnet(i+1) est le flux net traversant le plafond de la couche i (+ vers le haut) |
|---|
| 63 | do p=0,klev ! numero de la couche |
|---|
| 64 | zlnet(p+1) = 0.0 |
|---|
| 65 | do j=p+1,klev+1 |
|---|
| 66 | do i=0,p |
|---|
| 67 | zlnet(p+1) = zlnet(p+1)+ psi(i,j) |
|---|
| 68 | enddo |
|---|
| 69 | enddo |
|---|
| 70 | enddo |
|---|
| 71 | |
|---|
| 72 | c flux net au sol, + vers le haut: |
|---|
| 73 | PSOLLW = zlnet(1) |
|---|
| 74 | c flux vers le bas au sol, + vers le bas: |
|---|
| 75 | PSOLLWDN = 0.0 |
|---|
| 76 | do i=1,klev+1 |
|---|
| 77 | PSOLLWDN = PSOLLWDN+max(psi(i,0),0.0) |
|---|
| 78 | enddo |
|---|
| 79 | |
|---|
| 80 | c dfluxnet = radiative budget (W m-2) |
|---|
| 81 | do p=0,klev ! numero de la couche |
|---|
| 82 | dzlnet(p) = 0.0 |
|---|
| 83 | do j=0,klev+1 |
|---|
| 84 | dzlnet(p) = dzlnet(p)+psi(p,j) |
|---|
| 85 | enddo |
|---|
| 86 | enddo |
|---|
| 87 | |
|---|
| 88 | c -------------------------------------- |
|---|
| 89 | c Interpolation in the GCM vertical grid |
|---|
| 90 | c -------------------------------------- |
|---|
| 91 | |
|---|
| 92 | c Flux net |
|---|
| 93 | c -------- |
|---|
| 94 | |
|---|
| 95 | do j=1,klev+1 |
|---|
| 96 | ZFLNET(j) = zlnet(j) |
|---|
| 97 | enddo |
|---|
| 98 | PTOPLW = ZFLNET(klev+1) |
|---|
| 99 | |
|---|
| 100 | c Heating rates |
|---|
| 101 | c ------------- |
|---|
| 102 | |
|---|
| 103 | c cool (K/s) = dfluxnet (W/m2) ! positif quand nrj sort de la couche |
|---|
| 104 | c *g (m/s2) |
|---|
| 105 | c /(-dp) (epaisseur couche, en Pa=kg/m/s2) |
|---|
| 106 | c /cp (J/kg/K) |
|---|
| 107 | |
|---|
| 108 | do j=1,klev |
|---|
| 109 | pdp(j)=(PPB(j)-PPB(j+1))*1.e5 |
|---|
| 110 | enddo |
|---|
| 111 | |
|---|
| 112 | c calcul direct OU calcul par schema implicit |
|---|
| 113 | if (1.eq.1) then |
|---|
| 114 | do j=1,klev |
|---|
| 115 | ! ADAPTATION GCM POUR CP(T) |
|---|
| 116 | PCOOL(j) = dzlnet(j) *RG/cpdet(pt(j)) / pdp(j) |
|---|
| 117 | enddo |
|---|
| 118 | else |
|---|
| 119 | call lwi(klev,dzlnet,deltapsi,pdp,pt,PCOOL) |
|---|
| 120 | endif |
|---|
| 121 | c print*,dzlnet |
|---|
| 122 | c print*,pdp |
|---|
| 123 | c print*,PCOOL |
|---|
| 124 | |
|---|
| 125 | return |
|---|
| 126 | end |
|---|
| 127 | |
|---|
