Changeset 847


Ignore:
Timestamp:
04/27/16 17:37:34 (6 years ago)
Author:
ymipsl
Message:

bugfix : automatic generation of 2D regular domain.

YM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • XIOS/trunk/src/node/domain.cpp

    r836 r847  
    545545      if (bounds_lat_start >= 90 ) bounds_lat_start=90 ; 
    546546      else if (bounds_lat_start <= -90 ) bounds_lat_start=-90 ; 
    547       else if (bounds_lat_start <= 90 && bounds_lat_start >= lon(0)) 
    548       { 
    549         if ( 90-bounds_lat_start <= 0.1*std::abs(latStepStart)) bounds_lat_start=90 ; 
    550       } 
    551       else if (bounds_lat_start >= -90 && bounds_lat_start <= lon(0)) 
    552       { 
    553         if ( -90 + bounds_lat_start <= 0.1*std::abs(latStepStart)) bounds_lat_start=-90 ; 
     547      else if (bounds_lat_start <= 90 && bounds_lat_start >= lat(0)) 
     548      { 
     549        if ( std::abs(90-bounds_lat_start) <= 0.1*std::abs(latStepStart)) bounds_lat_start=90 ; 
     550      } 
     551      else if (bounds_lat_start >= -90 && bounds_lat_start <= lat(0)) 
     552      { 
     553        if ( std::abs(-90 - bounds_lat_start) <= 0.1*std::abs(latStepStart)) bounds_lat_start=-90 ; 
    554554      } 
    555555    } 
     
    563563      if (bounds_lat_end >= 90 ) bounds_lat_end=90 ; 
    564564      else if (bounds_lat_end <= -90 ) bounds_lat_end=-90 ; 
    565       else if (bounds_lat_end <= 90 && bounds_lat_end >= lon(nj_glo-1)) 
    566       { 
    567         if ( 90-bounds_lat_end <= 0.1*std::abs(latStepEnd)) bounds_lat_end=90 ; 
    568       } 
    569       else if (bounds_lat_end >= -90 && bounds_lat_end <= lon(nj_glo-1)) 
    570       { 
    571         if ( -90 + bounds_lat_end <= 0.1*std::abs(latStepEnd)) bounds_lat_end=-90 ; 
     565      else if (bounds_lat_end <= 90 && bounds_lat_end >= lat(nj_glo-1)) 
     566      { 
     567        if ( std::abs(90-bounds_lat_end) <= 0.1*std::abs(latStepEnd)) bounds_lat_end=90 ; 
     568      } 
     569      else if (bounds_lat_end >= -90 && bounds_lat_end <= lat(nj_glo-1)) 
     570      { 
     571        if ( std::abs(-90 - bounds_lat_end) <= 0.1*std::abs(latStepEnd)) bounds_lat_end=-90 ; 
    572572      } 
    573573    } 
Note: See TracChangeset for help on using the changeset viewer.