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

Last change on this file since 3977 was 3977, checked in by jbclement, 42 hours 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
Line 
1MODULE compute_tend_mod
2
3implicit none
4
5!=======================================================================
6contains
7!=======================================================================
8
9SUBROUTINE compute_tend(ngrid,nslope,min_ice,d_ice)
10
11implicit none
12
13!=======================================================================
14!
15! Compute the initial tendencies of the ice evolution based on the PCM data
16!
17!=======================================================================
18
19! Arguments
20! ---------
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
24real, dimension(ngrid,nslope), intent(out) :: d_ice ! Difference between the minima = evolution of perennial ice
25
26! Code
27! ----
28! We compute the difference
29d_ice = min_ice(:,:,2) - min_ice(:,:,1)
30
31! If the difference is too small, then there is no evolution
32where (abs(d_ice) < 1.e-10) d_ice = 0.
33
34! If the minimum over the last year is 0, then we have no perennial ice
35where (abs(min_ice(:,:,2)) < 1.e-10) d_ice = 0.
36
37END SUBROUTINE compute_tend
38
39END MODULE compute_tend_mod
Note: See TracBrowser for help on using the repository browser.