| 1 | SUBROUTINE heating(dist,rmu0,fract,sol_htg,swnet) |
|---|
| 2 | |
|---|
| 3 | |
|---|
| 4 | c======================================================================= |
|---|
| 5 | c |
|---|
| 6 | c Object: Computation of the solar heating rate |
|---|
| 7 | c SOL_HTG(klon,klev) |
|---|
| 8 | c |
|---|
| 9 | c Arguments: |
|---|
| 10 | c ---------- |
|---|
| 11 | c |
|---|
| 12 | c Input: |
|---|
| 13 | c ------ |
|---|
| 14 | c |
|---|
| 15 | c dist-----input-R- distance astronomique terre-soleil |
|---|
| 16 | c rmu0-----input-R- cosinus de l'angle zenithal |
|---|
| 17 | c fract----input-R- duree d'ensoleillement normalisee |
|---|
| 18 | c p(klon,nl) pressure (level) |
|---|
| 19 | c |
|---|
| 20 | c Output: |
|---|
| 21 | c ------- |
|---|
| 22 | c sol_htg-----output-R- echauffement atmospherique (visible) (K/s) |
|---|
| 23 | c swnet-------output-R- flux solaire net (+ vers le bas) (W/m2) |
|---|
| 24 | c |
|---|
| 25 | c======================================================================= |
|---|
| 26 | c----------------------------------------------------------------------- |
|---|
| 27 | c Declarations: |
|---|
| 28 | c ------------- |
|---|
| 29 | |
|---|
| 30 | use dimphy |
|---|
| 31 | IMPLICIT NONE |
|---|
| 32 | #include "dimensions.h" |
|---|
| 33 | |
|---|
| 34 | INTEGER NLEVEL,NLAYER,NSPECV |
|---|
| 35 | PARAMETER(NLAYER=llm,NLEVEL=NLAYER+1) |
|---|
| 36 | PARAMETER (NSPECV=24) |
|---|
| 37 | c |
|---|
| 38 | c ASTUCE POUR EVITER klon... EN ATTENDANT MIEUX |
|---|
| 39 | INTEGER ngrid |
|---|
| 40 | PARAMETER (ngrid=(jjm-1)*iim+2) ! = klon |
|---|
| 41 | c |
|---|
| 42 | |
|---|
| 43 | c Arguments: |
|---|
| 44 | c ---------- |
|---|
| 45 | |
|---|
| 46 | |
|---|
| 47 | real dist, rmu0(klon), fract(klon) |
|---|
| 48 | |
|---|
| 49 | real sol_htg(klon,klev) |
|---|
| 50 | real swnet(klon,klev+1) |
|---|
| 51 | |
|---|
| 52 | c Local: |
|---|
| 53 | c ------ |
|---|
| 54 | |
|---|
| 55 | INTEGER I,J,IG,K,IPRINT,ilat,nq |
|---|
| 56 | |
|---|
| 57 | c COMMONS for interface with local subroutines: |
|---|
| 58 | c --------------------------------------------- |
|---|
| 59 | |
|---|
| 60 | REAL UBARI,UBARV,UBAR0 |
|---|
| 61 | REAL CH4(NLEVEL),XN2(NLEVEL),H2(NLEVEL),AR(NLEVEL) |
|---|
| 62 | REAL XMU(NLEVEL),GAS1(NLAYER),COLDEN(NLAYER) |
|---|
| 63 | REAL FNETV(ngrid,NLEVEL),FUPV(ngrid,NLEVEL,NSPECV) |
|---|
| 64 | REAL FDV(ngrid,NLEVEL,NSPECV),FMNETV(ngrid,NLEVEL) |
|---|
| 65 | REAL CSUBP,RSFI,RSFV,F0PI |
|---|
| 66 | |
|---|
| 67 | COMMON /UBARED/ UBARI,UBARV,UBAR0 |
|---|
| 68 | |
|---|
| 69 | COMMON /GASS/ CH4,XN2 |
|---|
| 70 | & ,H2,AR |
|---|
| 71 | & ,XMU,GAS1 |
|---|
| 72 | & ,COLDEN |
|---|
| 73 | |
|---|
| 74 | COMMON /FLUXvV/ FNETV, |
|---|
| 75 | & FUPV, |
|---|
| 76 | & FDV, |
|---|
| 77 | & FMNETV |
|---|
| 78 | |
|---|
| 79 | COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI |
|---|
| 80 | |
|---|
| 81 | |
|---|
| 82 | c================================================================== |
|---|
| 83 | |
|---|
| 84 | fnetv = 0.0 |
|---|
| 85 | sol_htg= 0.0 |
|---|
| 86 | swnet = 0.0 |
|---|
| 87 | c pour sorties dans gfluxv... |
|---|
| 88 | iprint = 0 |
|---|
| 89 | |
|---|
| 90 | DO ig=1,klon |
|---|
| 91 | IF(fract(ig).LT.1.e-5) THEN |
|---|
| 92 | DO j=1,nlayer |
|---|
| 93 | sol_htg(ig,j)=0. |
|---|
| 94 | ENDDO |
|---|
| 95 | ELSE |
|---|
| 96 | ubar0=rmu0(ig) |
|---|
| 97 | |
|---|
| 98 | CALL sfluxv(iprint,ig,dist) ! #3 |
|---|
| 99 | |
|---|
| 100 | fnetv(ig,:) = fnetv(ig,:) *fract(ig) ! >0 vers le haut |
|---|
| 101 | c >0 vers le bas + conversion en W/m2: |
|---|
| 102 | swnet(ig,:) = -1.e-3*fnetv(ig,:) |
|---|
| 103 | |
|---|
| 104 | DO j=1,nlayer |
|---|
| 105 | sol_htg(ig,j)= ! K/s |
|---|
| 106 | s (fnetv(ig,j+1)-fnetv(ig,j)) |
|---|
| 107 | s /(colden(j)*csubp) |
|---|
| 108 | ENDDO |
|---|
| 109 | ENDIF |
|---|
| 110 | ENDDO |
|---|
| 111 | |
|---|
| 112 | RETURN |
|---|
| 113 | 191 FORMAT(F8.2,1P10E10.2) |
|---|
| 114 | 192 FORMAT(a8,1P10E10.2) |
|---|
| 115 | END |
|---|
