Ignore:
Timestamp:
Mar 5, 2014, 2:19:12 PM (10 years ago)
Author:
lguez
Message:

Converted to free source form files in libf/phylmd which were still in
fixed source form. The conversion was done using the polish mode of
the NAG Fortran Compiler.

In addition to converting to free source form, the processing of the
files also:

-- indented the code (including comments);

-- set Fortran keywords to uppercase, and set all other identifiers
to lower case;

-- added qualifiers to end statements (for example "end subroutine
conflx", instead of "end");

-- changed the terminating statements of all DO loops so that each
loop ends with an ENDDO statement (instead of a labeled continue).

-- replaced #include by include.

File:
1 moved

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phylmd/moy_undefSTD.F90

    r1988 r1992  
    1 !
     1
    22! $Id$
    3 !
    4       SUBROUTINE moy_undefSTD(itap,itapm1)
    5       USE netcdf
    6       USE dimphy
    7       USE phys_state_var_mod
    8       USE phys_cal_mod, only : mth_len
    9       IMPLICIT none
    10       include "clesphys.h"
    11 c
    12 c====================================================================
    13 c
    14 c I. Musat : 09.2004
    15 c
    16 c Moyenne - a des frequences differentes - des valeurs bien definies
    17 c        (.NE.missing_val) des variables interpolees a un niveau de
    18 c        pression.
    19 c 1) les variables de type "day" (nout=1) ou "mth" (nout=2) sont sommees
    20 c    tous les pas de temps de la physique
    21 c
    22 c 2) les variables de type "NMC" (nout=3) sont calculees a partir
    23 c    des valeurs instantannees toutes les 6 heures
    24 c
    25 c
    26 c NB: mettre "inst(X)" dans le write_hist*NMC.h !
    27 c====================================================================
    28 cym#include "dimensions.h"
    29 cym      integer jjmp1
    30 cym      parameter (jjmp1=jjm+1-1/jjm)
    31 cym#include "dimphy.h"
    32 c
    33 c
    34 c variables Input
    35 c    INTEGER nlevSTD, klevSTD, itap
    36 c    PARAMETER(klevSTD=17)
    37       INTEGER itap, itapm1
    38 c
    39 c variables locales
    40 c    INTEGER i, k, nout, n
    41 c    PARAMETER(nout=3) !nout=1 day/nout=2 mth/nout=3 NMC
    42       INTEGER i, k, n
    43 c    REAL dtime, freq_outNMC(nout), freq_moyNMC(nout)
    44 c    REAL freq_outNMC(nout), freq_calNMC(nout)
    45       REAL freq_moyNMC(nout)
    46 c
    47 c variables Output
    48 c    REAL tnondef(klon,klevSTD,nout)
    49 c    REAL tsumSTD(klon,klevSTD,nout)
    50 c
    51       REAL un_jour
    52       PARAMETER(un_jour=86400.)
    53       REAL missing_val
    54 c
    55       missing_val=nf90_fill_real
    56 c
    57       DO n=1, nout
    58        IF(freq_outNMC(n).LT.0) THEN
    59          freq_moyNMC(n)=(mth_len*un_jour)/freq_calNMC(n)
    60 c        print*,'moy_undefSTD n freq_out freq_moy =',
    61 c    $n,freq_moyNMC(n)
    62        ELSE
    63          freq_moyNMC(n)=freq_outNMC(n)/freq_calNMC(n)
    64        ENDIF         
    65 c
    66 c calcul 1 fois pas mois, 1 fois par jour ou toutes les 6h
    67 c
    68        IF(n.EQ.1.AND.itap.EQ.itapm1.OR.
    69      $n.GT.1.AND.MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.0) THEN
    70 c
    71 c      print*,'moy_undefSTD n itap itapm1',n,itap,itapm1
    72 c
    73         DO k=1, nlevSTD
    74          DO i=1, klon
    75           IF(tnondef(i,k,n).NE.(freq_moyNMC(n))) THEN
    76            tsumSTD(i,k,n)=tsumSTD(i,k,n)/
    77      $     (freq_moyNMC(n)-tnondef(i,k,n))
    78           usumSTD(i,k,n)=usumSTD(i,k,n)/
    79      $     (freq_moyNMC(n)-tnondef(i,k,n))
    80           vsumSTD(i,k,n)=vsumSTD(i,k,n)/
    81      $     (freq_moyNMC(n)-tnondef(i,k,n))
    82           wsumSTD(i,k,n)=wsumSTD(i,k,n)/
    83      $     (freq_moyNMC(n)-tnondef(i,k,n))
    84           phisumSTD(i,k,n)=phisumSTD(i,k,n)/
    85      $     (freq_moyNMC(n)-tnondef(i,k,n))
    86           qsumSTD(i,k,n)=qsumSTD(i,k,n)/
    87      $     (freq_moyNMC(n)-tnondef(i,k,n))
    88           rhsumSTD(i,k,n)=rhsumSTD(i,k,n)/
    89      $     (freq_moyNMC(n)-tnondef(i,k,n))
    90           uvsumSTD(i,k,n)=uvsumSTD(i,k,n)/
    91      $     (freq_moyNMC(n)-tnondef(i,k,n))
    92           vqsumSTD(i,k,n)=vqsumSTD(i,k,n)/
    93      $     (freq_moyNMC(n)-tnondef(i,k,n))
    94           vTsumSTD(i,k,n)=vTsumSTD(i,k,n)/
    95      $     (freq_moyNMC(n)-tnondef(i,k,n))
    96           wqsumSTD(i,k,n)=wqsumSTD(i,k,n)/
    97      $     (freq_moyNMC(n)-tnondef(i,k,n))
    98           vphisumSTD(i,k,n)=vphisumSTD(i,k,n)/
    99      $     (freq_moyNMC(n)-tnondef(i,k,n))
    100           wTsumSTD(i,k,n)=wTsumSTD(i,k,n)/
    101      $     (freq_moyNMC(n)-tnondef(i,k,n))
    102           u2sumSTD(i,k,n)=u2sumSTD(i,k,n)/
    103      $     (freq_moyNMC(n)-tnondef(i,k,n))
    104           v2sumSTD(i,k,n)=v2sumSTD(i,k,n)/
    105      $     (freq_moyNMC(n)-tnondef(i,k,n))
    106           T2sumSTD(i,k,n)=T2sumSTD(i,k,n)/
    107      $     (freq_moyNMC(n)-tnondef(i,k,n))
    108           O3sumSTD(i,k,n)=O3sumSTD(i,k,n)/
    109      $     (freq_moyNMC(n)-tnondef(i,k,n))
    110           O3daysumSTD(i,k,n)=O3daysumSTD(i,k,n)/
    111      $     (freq_moyNMC(n)-tnondef(i,k,n))
     3
     4SUBROUTINE moy_undefstd(itap, itapm1)
     5  USE netcdf
     6  USE dimphy
     7  USE phys_state_var_mod
     8  USE phys_cal_mod, ONLY: mth_len
     9  IMPLICIT NONE
     10  include "clesphys.h"
     11
     12  ! ====================================================================
     13
     14  ! I. Musat : 09.2004
     15
     16  ! Moyenne - a des frequences differentes - des valeurs bien definies
     17  ! (.NE.missing_val) des variables interpolees a un niveau de
     18  ! pression.
     19  ! 1) les variables de type "day" (nout=1) ou "mth" (nout=2) sont sommees
     20  ! tous les pas de temps de la physique
     21
     22  ! 2) les variables de type "NMC" (nout=3) sont calculees a partir
     23  ! des valeurs instantannees toutes les 6 heures
     24
     25
     26  ! NB: mettre "inst(X)" dans le write_hist*NMC.h !
     27  ! ====================================================================
     28  ! ym#include "dimensions.h"
     29  ! ym      integer jjmp1
     30  ! ym      parameter (jjmp1=jjm+1-1/jjm)
     31  ! ym#include "dimphy.h"
     32
     33
     34  ! variables Input
     35  ! INTEGER nlevSTD, klevSTD, itap
     36  ! PARAMETER(klevSTD=17)
     37  INTEGER itap, itapm1
     38
     39  ! variables locales
     40  ! INTEGER i, k, nout, n
     41  ! PARAMETER(nout=3) !nout=1 day/nout=2 mth/nout=3 NMC
     42  INTEGER i, k, n
     43  ! REAL dtime, freq_outNMC(nout), freq_moyNMC(nout)
     44  ! REAL freq_outNMC(nout), freq_calNMC(nout)
     45  REAL freq_moynmc(nout)
     46
     47  ! variables Output
     48  ! REAL tnondef(klon,klevSTD,nout)
     49  ! REAL tsumSTD(klon,klevSTD,nout)
     50
     51  REAL un_jour
     52  PARAMETER (un_jour=86400.)
     53  REAL missing_val
     54
     55  missing_val = nf90_fill_real
     56
     57  DO n = 1, nout
     58    IF (freq_outnmc(n)<0) THEN
     59      freq_moynmc(n) = (mth_len*un_jour)/freq_calnmc(n)
     60      ! print*,'moy_undefSTD n freq_out freq_moy =',
     61      ! $n,freq_moyNMC(n)
     62    ELSE
     63      freq_moynmc(n) = freq_outnmc(n)/freq_calnmc(n)
     64    END IF
     65
     66    ! calcul 1 fois pas mois, 1 fois par jour ou toutes les 6h
     67
     68    IF (n==1 .AND. itap==itapm1 .OR. n>1 .AND. mod(itap,nint(freq_outnmc(n)/ &
     69        dtime))==0) THEN
     70
     71      ! print*,'moy_undefSTD n itap itapm1',n,itap,itapm1
     72
     73      DO k = 1, nlevstd
     74        DO i = 1, klon
     75          IF (tnondef(i,k,n)/=(freq_moynmc(n))) THEN
     76            tsumstd(i, k, n) = tsumstd(i, k, n)/(freq_moynmc(n)-tnondef(i,k,n &
     77              ))
     78            usumstd(i, k, n) = usumstd(i, k, n)/(freq_moynmc(n)-tnondef(i,k,n &
     79              ))
     80            vsumstd(i, k, n) = vsumstd(i, k, n)/(freq_moynmc(n)-tnondef(i,k,n &
     81              ))
     82            wsumstd(i, k, n) = wsumstd(i, k, n)/(freq_moynmc(n)-tnondef(i,k,n &
     83              ))
     84            phisumstd(i, k, n) = phisumstd(i, k, n)/ &
     85              (freq_moynmc(n)-tnondef(i,k,n))
     86            qsumstd(i, k, n) = qsumstd(i, k, n)/(freq_moynmc(n)-tnondef(i,k,n &
     87              ))
     88            rhsumstd(i, k, n) = rhsumstd(i, k, n)/(freq_moynmc(n)-tnondef(i,k &
     89              ,n))
     90            uvsumstd(i, k, n) = uvsumstd(i, k, n)/(freq_moynmc(n)-tnondef(i,k &
     91              ,n))
     92            vqsumstd(i, k, n) = vqsumstd(i, k, n)/(freq_moynmc(n)-tnondef(i,k &
     93              ,n))
     94            vtsumstd(i, k, n) = vtsumstd(i, k, n)/(freq_moynmc(n)-tnondef(i,k &
     95              ,n))
     96            wqsumstd(i, k, n) = wqsumstd(i, k, n)/(freq_moynmc(n)-tnondef(i,k &
     97              ,n))
     98            vphisumstd(i, k, n) = vphisumstd(i, k, n)/ &
     99              (freq_moynmc(n)-tnondef(i,k,n))
     100            wtsumstd(i, k, n) = wtsumstd(i, k, n)/(freq_moynmc(n)-tnondef(i,k &
     101              ,n))
     102            u2sumstd(i, k, n) = u2sumstd(i, k, n)/(freq_moynmc(n)-tnondef(i,k &
     103              ,n))
     104            v2sumstd(i, k, n) = v2sumstd(i, k, n)/(freq_moynmc(n)-tnondef(i,k &
     105              ,n))
     106            t2sumstd(i, k, n) = t2sumstd(i, k, n)/(freq_moynmc(n)-tnondef(i,k &
     107              ,n))
     108            o3sumstd(i, k, n) = o3sumstd(i, k, n)/(freq_moynmc(n)-tnondef(i,k &
     109              ,n))
     110            o3daysumstd(i, k, n) = o3daysumstd(i, k, n)/ &
     111              (freq_moynmc(n)-tnondef(i,k,n))
    112112          ELSE
    113            tsumSTD(i,k,n)=missing_val
    114            usumSTD(i,k,n)=missing_val
    115            vsumSTD(i,k,n)=missing_val
    116            wsumSTD(i,k,n)=missing_val
    117            phisumSTD(i,k,n)=missing_val
    118            qsumSTD(i,k,n)=missing_val
    119            rhsumSTD(i,k,n)=missing_val
    120            uvsumSTD(i,k,n)=missing_val
    121            vqsumSTD(i,k,n)=missing_val
    122            vTsumSTD(i,k,n)=missing_val
    123            wqsumSTD(i,k,n)=missing_val
    124            vphisumSTD(i,k,n)=missing_val
    125            wTsumSTD(i,k,n)=missing_val
    126            u2sumSTD(i,k,n)=missing_val
    127            v2sumSTD(i,k,n)=missing_val
    128            T2sumSTD(i,k,n)=missing_val
    129            O3sumSTD(i,k,n)=missing_val
    130            O3daysumSTD(i,k,n)=missing_val
    131           ENDIF !tnondef(i,k,n).NE.(freq_moyNMC(n))
    132          ENDDO !i
    133         ENDDO !k
    134        ENDIF !MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.0
    135 c
    136       ENDDO !n
    137 c
    138       RETURN
    139       END 
     113            tsumstd(i, k, n) = missing_val
     114            usumstd(i, k, n) = missing_val
     115            vsumstd(i, k, n) = missing_val
     116            wsumstd(i, k, n) = missing_val
     117            phisumstd(i, k, n) = missing_val
     118            qsumstd(i, k, n) = missing_val
     119            rhsumstd(i, k, n) = missing_val
     120            uvsumstd(i, k, n) = missing_val
     121            vqsumstd(i, k, n) = missing_val
     122            vtsumstd(i, k, n) = missing_val
     123            wqsumstd(i, k, n) = missing_val
     124            vphisumstd(i, k, n) = missing_val
     125            wtsumstd(i, k, n) = missing_val
     126            u2sumstd(i, k, n) = missing_val
     127            v2sumstd(i, k, n) = missing_val
     128            t2sumstd(i, k, n) = missing_val
     129            o3sumstd(i, k, n) = missing_val
     130            o3daysumstd(i, k, n) = missing_val
     131          END IF !tnondef(i,k,n).NE.(freq_moyNMC(n))
     132        END DO !i
     133      END DO !k
     134    END IF !MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.0
     135
     136  END DO !n
     137
     138  RETURN
     139END SUBROUTINE moy_undefstd
Note: See TracChangeset for help on using the changeset viewer.