source: trunk/LMDZ.COMMON/libf/evolution/compute_tend_mod.F90 @ 3984

Last change on this file since 3984 was 3977, checked in by jbclement, 3 weeks ago

PEM:

  • All operations computed by the PEM on the PCM data (averages, minima) are now performed by XIOS with two dedicated file. One is for daily operation outputs ("Xoutdaily4pem*.nc") and the other for yearly operation outputs ("Xoutyearly4pem*.nc").
  • Deletion of the reshaping tool "reshape_XIOS_output" used to convert XIOS outputs onto the PCM grid. Thus, the PEM is now able to read directly the format of XIOS outputs.
  • Addition of subroutines to convert data between a lon x lat array and a vector.

JBC
File size: 1.2 KB
RevLine 
[3149]1MODULE compute_tend_mod
[3076]2
3implicit none
4
5!=======================================================================
6contains
7!=======================================================================
8
[3498]9SUBROUTINE compute_tend(ngrid,nslope,min_ice,d_ice)
[2779]10
[3070]11implicit none
[2779]12
13!=======================================================================
14!
[3367]15! Compute the initial tendencies of the ice evolution based on the PCM data
[2779]16!
17!=======================================================================
18
[3977]19! Arguments
20! ---------
[3149]21integer,                         intent(in) :: ngrid   ! # of grid points
22integer,                         intent(in) :: nslope  ! # of subslopes
23real, dimension(ngrid,nslope,2), intent(in) :: min_ice ! Minima of ice at each point for the PCM years
[3977]24real, dimension(ngrid,nslope), intent(out) :: d_ice ! Difference between the minima = evolution of perennial ice
[2779]25
[3977]26! Code
27! ----
[3149]28! We compute the difference
[3498]29d_ice = min_ice(:,:,2) - min_ice(:,:,1)
[2779]30
[3149]31! If the difference is too small, then there is no evolution
[3498]32where (abs(d_ice) < 1.e-10) d_ice = 0.
[2779]33
[3367]34! If the minimum over the last year is 0, then we have no perennial ice
[3498]35where (abs(min_ice(:,:,2)) < 1.e-10) d_ice = 0.
[3367]36
[3149]37END SUBROUTINE compute_tend
[2779]38
[3149]39END MODULE compute_tend_mod
Note: See TracBrowser for help on using the repository browser.