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