Ignore:
Timestamp:
07/09/12 15:23:38 (12 years ago)
Author:
ymipsl
Message:

Update on 3D dynamic

YM

File:
1 edited

Legend:

Unmodified
Added
Removed

  • codes/icosagcm/trunk/src/etat0_jablonowsky06.f90

    r12 r15  
    140140  REAL(rstd) :: r2 
    141141  REAL(rstd) :: utot 
     142  REAL(rstd) :: lonx,latx 
    142143  
    143144    DO l=1,llm 
    … …  
    157158     
    158159    CALL lonlat2xyz(Pi/9,2*Pi/9,V0) 
     160 
    159161    u(:,:)=1e10       
    160162    DO l=1,llm 
    … …  
    165167          CALL xyz2lonlat(xyz_e(ij+u_right,:)/radius,lon,lat) 
    166168          CALL cross_product2(V0,xyz_e(ij+u_right,:)/radius,ep) 
    167           r2=(asin(sum(ep*ep)))**2 
     169          r2=(asin(sqrt(sum(ep*ep))))**2 
    168170          utot=u0*cos(etav(l))**1.5*sin(2*lat)**2 + up0*exp(-r2/0.01) 
    169171 
    … …  
    184186          CALL xyz2lonlat(xyz_e(ij+u_lup,:)/radius,lon,lat) 
    185187          CALL cross_product2(V0,xyz_e(ij+u_lup,:)/radius,ep) 
    186           r2=(asin(sum(ep*ep)))**2 
     188          r2=(asin(sqrt(sum(ep*ep))))**2 
    187189          utot=u0*cos(etav(l))**1.5*sin(2*lat)**2 + up0*exp(-r2/0.01) 
    188190          ulon(1) = -sin(lon) * utot 
    … …  
    205207          CALL xyz2lonlat(xyz_e(ij+u_ldown,:)/radius,lon,lat) 
    206208          CALL cross_product2(V0,xyz_e(ij+u_ldown,:)/radius,ep) 
    207           r2=(asin(sum(ep*ep)))**2 
     209          r2=(asin(sqrt(sum(ep*ep))))**2 
    208210          utot=u0*cos(etav(l))**1.5*sin(2*lat)**2 + up0*exp(-r2/0.01) 
    209211          ulon(1) = -sin(lon) * utot 
Note: See TracChangeset for help on using the changeset viewer.