/[lmdze]/trunk/Sources/phylmd/Orography/orosetup.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/Orography/orosetup.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 168 by guez, Wed Apr 29 15:47:56 2015 UTC revision 169 by guez, Mon Sep 14 17:13:16 2015 UTC
# Line 10  SUBROUTINE orosetup(nlon, ktest, kkcrit, Line 10  SUBROUTINE orosetup(nlon, ktest, kkcrit,
10    
11    USE dimens_m    USE dimens_m
12    USE dimphy    USE dimphy
13      use nr_util, only: pi
14    USE suphec_m    USE suphec_m
15    USE yoegwd    USE yoegwd
16    
# Line 64  SUBROUTINE orosetup(nlon, ktest, kkcrit, Line 65  SUBROUTINE orosetup(nlon, ktest, kkcrit,
65    !old  zcons2=g**2/cpd    !old  zcons2=g**2/cpd
66    zcons2 = rg**2/rcpd    zcons2 = rg**2/rcpd
67    !old  zcons3=1.5*api    !old  zcons3=1.5*api
68    zcons3 = 1.5*rpi    zcons3 = 1.5*pi
69    
70    ! 2.    ! 2.
71    
# Line 203  SUBROUTINE orosetup(nlon, ktest, kkcrit, Line 204  SUBROUTINE orosetup(nlon, ktest, kkcrit,
204          zu = pulow(jl)          zu = pulow(jl)
205       END IF       END IF
206       zphi = atan(pvlow(jl)/zu)       zphi = atan(pvlow(jl)/zu)
207       ppsi(jl, klev+1) = ptheta(jl)*rpi/180. - zphi       ppsi(jl, klev+1) = ptheta(jl)*pi/180. - zphi
208       zb(jl) = 1. - 0.18*pgamma(jl) - 0.04*pgamma(jl)**2       zb(jl) = 1. - 0.18*pgamma(jl) - 0.04*pgamma(jl)**2
209       zc(jl) = 0.48*pgamma(jl) + 0.3*pgamma(jl)**2       zc(jl) = 0.48*pgamma(jl) + 0.3*pgamma(jl)**2
210       pd1(jl) = zb(jl) - (zb(jl)-zc(jl))*(sin(ppsi(jl, klev+1))**2)       pd1(jl) = zb(jl) - (zb(jl)-zc(jl))*(sin(ppsi(jl, klev+1))**2)
# Line 337  SUBROUTINE orosetup(nlon, ktest, kkcrit, Line 338  SUBROUTINE orosetup(nlon, ktest, kkcrit,
338             zrhop = prho(jl, jk+1)             zrhop = prho(jl, jk+1)
339             znup(jl) = znup(jl) + (zdelp/rg)*((zstabp/zrhop+zstabm/zrhom)/2.)/ &             znup(jl) = znup(jl) + (zdelp/rg)*((zstabp/zrhop+zstabm/zrhom)/2.)/ &
340                  zwind                  zwind
341             IF ((znum(jl)<=rpi/2.) .AND. (znup(jl)>rpi/2.) .AND. (kkcrith( &             IF ((znum(jl)<=pi/2.) .AND. (znup(jl)>pi/2.) .AND. (kkcrith( &
342                  jl)==klev)) kkcrith(jl) = jk                  jl)==klev)) kkcrith(jl) = jk
343    
344          END IF          END IF
# Line 362  SUBROUTINE orosetup(nlon, ktest, kkcrit, Line 363  SUBROUTINE orosetup(nlon, ktest, kkcrit,
363                zu = pum1(jl, jk)                zu = pum1(jl, jk)
364             END IF             END IF
365             zphi = atan(pvm1(jl, jk)/zu)             zphi = atan(pvm1(jl, jk)/zu)
366             ppsi(jl, jk) = ptheta(jl)*rpi/180. - zphi             ppsi(jl, jk) = ptheta(jl)*pi/180. - zphi
367          END IF          END IF
368       end DO       end DO
369    end DO    end DO
Legend:
Removed from v.168  
changed lines
  Added in v.169
  ViewVC Help
Powered by ViewVC 1.1.21