Commit 7a99b3b4 authored by ldemaine's avatar ldemaine
Browse files

bug fix 2 (problème de varible)

parent 33757f08
......@@ -1503,14 +1503,14 @@ c******************************
endif ! of if (perturkey.eq.5)
c 4.3.5 Computation of high precision wind
write(*,*) 'temp', temp
if (hireskey.eq.1) then
call slope_wind(dset,xlon,xlat,zsurface,zradius,
& temp,temp_gcm,sheight,
& zonwind,merwind,old_zonwind,old_merwind,
& tl,temp_gcm,sheight,
& ul,vl,old_zonwind,old_merwind,
& upslope_wind,crossslope_wind)
c$$$ call slope_wind(dset,xlon,xlat,zsurface,z_0,theta_s,
......@@ -7755,7 +7755,7 @@ c$$$ real,intent(in) :: pratio
parameter (degtorad=pi/180.0d0)
open(1, file='Test_slope/TEST_4.txt', status = 'new')
open(1, file='Test_slope/TEST_5.txt', status = 'new')
write(1,*) ''
write(1,*) '|-----------------------------------------|'
......@@ -7879,7 +7879,8 @@ c$$$ & /(1.-(sigma(levhi_0)/sigma(levlow_0)))
write(1,*) '| Cd', Cd
upslope_wind = sign(abs(sqrt(Float / Cd)), Float/Cd)
upslope_wind = sign(sqrt(abs(Float / Cd)), Float/Cd)
c$$$ upslope_wind = sign(abs(sqrt(Float / Cd)), Float/Cd)
write(1,*) '| upslope_wind', upslope_wind
c$$$ do i=0,1000
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment