Ignore:
Timestamp:
10/20/14 23:42:26 (10 years ago)
Author:
dubos
Message:

Partial etat0 cleanup (removed calls to xyz2lonlat)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • codes/icosagcm/trunk/src/etat0_dcmip2.f90

    r186 r286  
    8686          DO i=ii_begin,ii_end 
    8787             ij=(j-1)*iim+i 
    88              CALL comp_all(xyz_i(ij,:), ps(ij),phis(ij),temp(ij,l), dummy1,dummy2) 
     88             CALL comp_all(lon_i(ij),lat_i(ij), ps(ij),phis(ij),temp(ij,l), dummy1,dummy2) 
    8989          END DO 
    9090       END DO 
     
    100100          DO i=ii_begin,ii_end 
    101101             ij=(j-1)*iim+i 
    102              CALL comp_all(xyz_e(ij+u_right,:), dummy1,dummy2,dummy3, ulon(ij+u_right,l),ulat(ij+u_right,l)) 
    103              CALL comp_all(xyz_e(ij+u_lup,:), dummy1,dummy2,dummy3, ulon(ij+u_lup,l),ulat(ij+u_lup,l)) 
    104              CALL comp_all(xyz_e(ij+u_ldown,:), dummy1,dummy2,dummy3, ulon(ij+u_ldown,l),ulat(ij+u_ldown,l)) 
     102             CALL comp_all(lon_e(ij+u_right), lat_e(ij+u_right), & 
     103                  dummy1,dummy2,dummy3, ulon(ij+u_right,l),ulat(ij+u_right,l)) 
     104             CALL comp_all(lon_e(ij+u_lup), lat_e(ij+u_lup), & 
     105                  dummy1,dummy2,dummy3, ulon(ij+u_lup,l),ulat(ij+u_lup,l)) 
     106             CALL comp_all(lon_e(ij+u_ldown), lat_e(ij+u_ldown), & 
     107                  dummy1,dummy2,dummy3, ulon(ij+u_ldown,l),ulat(ij+u_ldown,l)) 
    105108          END DO 
    106109       END DO 
     
    112115    CONTAINS 
    113116       
    114       SUBROUTINE comp_all(xyz, psj,phisj,tempj, ulonj,ulatj) 
     117      SUBROUTINE comp_all(lon,lat, psj,phisj,tempj, ulonj,ulatj) 
    115118        USE dcmip_initial_conditions_test_1_2_3 
    116         REAL(rstd), INTENT(IN) :: xyz(3) 
     119        REAL(rstd), INTENT(IN) :: lon, lat 
    117120        REAL(rstd), INTENT(OUT) :: psj,phisj,tempj,ulonj,ulatj 
    118         REAL :: lon,lat, dummy 
     121        REAL :: dummy 
    119122        dummy=0. 
    120         CALL xyz2lonlat(xyz,lon,lat) 
    121123        SELECT CASE (icase) 
    122124        CASE(0) 
Note: See TracChangeset for help on using the changeset viewer.