Changeset 2222 for LMDZ5/trunk/libf


Ignore:
Timestamp:
Mar 10, 2015, 10:58:04 AM (9 years ago)
Author:
Ehouarn Millour
Message:

Bug fix: Poles are single points on physics grid. True mesh area there is thus the sum of corresponding dynamics "polar meshes" areas.
Note that this also implies that some extra work is needed in physics to generate correct areas at "polar points" for outputs (especially for zoomed grids).
EM

Location:
LMDZ5/trunk/libf
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/dyn3d/gcm.F

    r2221 r2222  
    435435         zcvfi(ngridmx) = cv(ip1jm-iim)
    436436         CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,aire,airefi)
    437          WRITE(lunout,*)
    438      .       'GCM: WARNING!!! vitesse verticale nulle dans la physique'
     437         ! Poles are single points on physics grid
     438         airefi(1)=sum(aire(1:iim))
     439         airefi(ngridmx)=sum(aire(ip1jmp1-(iim+1)+1:ip1jmp1-1))
     440!         WRITE(lunout,*)
     441!     .       'GCM: WARNING!!! vitesse verticale nulle dans la physique'
    439442! Physics:
    440443#ifdef CPP_PHYS

  • LMDZ5/trunk/libf/dyn3dmem/gcm.F

    r2221 r2222  
    438438         zcvfi(ngridmx) = cv(ip1jm-iim)
    439439         CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,aire,airefi)
    440 
    441          WRITE(lunout,*)
    442      .       'GCM: WARNING!!! vitesse verticale nulle dans la physique'
     440         ! Poles are single points on physics grid
     441         airefi(1)=sum(aire(1:iim))
     442         airefi(ngridmx)=sum(aire(ip1jmp1-(iim+1)+1:ip1jmp1-1))
     443!         WRITE(lunout,*)
     444!     .       'GCM: WARNING!!! vitesse verticale nulle dans la physique'
    443445! Physics:
    444446#ifdef CPP_PHYS

  • LMDZ5/trunk/libf/dyn3dpar/gcm.F

    r2221 r2222  
    437437         zcvfi(ngridmx) = cv(ip1jm-iim)
    438438         CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,aire,airefi)
    439 
    440          WRITE(lunout,*)
    441      .       'GCM: WARNING!!! vitesse verticale nulle dans la physique'
     439         ! Poles are single points on physics grid
     440         airefi(1)=sum(aire(1:iim))
     441         airefi(ngridmx)=sum(aire(ip1jmp1-(iim+1)+1:ip1jmp1-1))
     442!         WRITE(lunout,*)
     443!     .       'GCM: WARNING!!! vitesse verticale nulle dans la physique'
    442444! Physics:
    443445#ifdef CPP_PHYS
Note: See TracChangeset for help on using the changeset viewer.