| 1 | subroutine nlthermeq(nlon, nlayer, pplev, pplay) |
|---|
| 2 | c |
|---|
| 3 | c Compute the number of layers nlaylte (stored in common yomlw.h) |
|---|
| 4 | c over which local thermodynamic equilibrium |
|---|
| 5 | c radiation scheme should be run to be sure of covering at least to a |
|---|
| 6 | c height greater than (pressure lower than) p=pminte, set in nlteparams.h. |
|---|
| 7 | c The maximum layer needed is found for the worst possible case. |
|---|
| 8 | c Stephen Lewis 6/2000 |
|---|
| 9 | c Modified Y. Wanherdrick/ F. Forget 09/2000 |
|---|
| 10 | use dimphy |
|---|
| 11 | implicit none |
|---|
| 12 | c#include "dimradmars.h" |
|---|
| 13 | #include "nlteparams.h" |
|---|
| 14 | c#include "yomlw.h" |
|---|
| 15 | #include "clesphys.h" |
|---|
| 16 | |
|---|
| 17 | c |
|---|
| 18 | c Input: |
|---|
| 19 | integer nlon, nlayer |
|---|
| 20 | real pplev(nlon, nlayer+1) |
|---|
| 21 | real pplay(nlon, nlayer) |
|---|
| 22 | c |
|---|
| 23 | c Local: |
|---|
| 24 | integer igpmax, ismax |
|---|
| 25 | logical firstcall |
|---|
| 26 | data firstcall /.true./ |
|---|
| 27 | save firstcall, igpmax |
|---|
| 28 | |
|---|
| 29 | |
|---|
| 30 | INTEGER i,ix |
|---|
| 31 | real sxmax |
|---|
| 32 | |
|---|
| 33 | ccc |
|---|
| 34 | if(firstcall) then |
|---|
| 35 | c Find the location of maximum surface pressure. |
|---|
| 36 | c Location won't vary much so only do it at the start; |
|---|
| 37 | c with no topography location would vary, but this is only |
|---|
| 38 | c needed for an estimate so any point would do in that case. |
|---|
| 39 | ismax=1 |
|---|
| 40 | sxmax=pplev(1,1) |
|---|
| 41 | ix=1 |
|---|
| 42 | do i=1,nlon-1 |
|---|
| 43 | if(pplev(i,ix).gt.sxmax) then |
|---|
| 44 | sxmax=pplev(i,ix) |
|---|
| 45 | ismax=i+1 |
|---|
| 46 | endif |
|---|
| 47 | enddo |
|---|
| 48 | |
|---|
| 49 | |
|---|
| 50 | igpmax = ismax ! longitude/ latitude where pression is maximum |
|---|
| 51 | write(*, 10) ptrans |
|---|
| 52 | write(*, 20) zw |
|---|
| 53 | write(*, 30) pminte |
|---|
| 54 | firstcall = .false. |
|---|
| 55 | endif |
|---|
| 56 | |
|---|
| 57 | IF(callnlte .or. callnirco2) THEN |
|---|
| 58 | c Find first layer above pminte at this location |
|---|
| 59 | do nlaylte = nlayer, 1, -1 |
|---|
| 60 | if (pplay(igpmax, nlaylte).gt.pminte) go to 100 |
|---|
| 61 | enddo |
|---|
| 62 | ELSE |
|---|
| 63 | nlaylte=nlayer |
|---|
| 64 | END IF |
|---|
| 65 | |
|---|
| 66 | c write(*,*) 'LTE rad. calculations up to layer ', nlaylte |
|---|
| 67 | |
|---|
| 68 | c |
|---|
| 69 | 100 return |
|---|
| 70 | c |
|---|
| 71 | |
|---|
| 72 | 10 format(' nlthermeq: transition to NLTE centred at ',f6.2,'Pa') |
|---|
| 73 | 20 format(' half-width (scale heights) ',f6.2) |
|---|
| 74 | 30 format(' suggested LTE coverage at least ',f6.2,'Pa') |
|---|
| 75 | 40 format(' nlthermeq: purely NLTE contribution over (nlayer) ',f6.2) |
|---|
| 76 | |
|---|
| 77 | end |
|---|
