Changeset 5789 for LMDZ6


Ignore:
Timestamp:
Jul 28, 2025, 1:12:04 PM (6 days ago)
Author:
fhourdin
Message:

Bug fix for initialisation of tnondef.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/ini_undefSTD.f90

    r5685 r5789  
    5959    ! initialisation variables en debut du mois, de la journee ou des 6h
    6060
    61     IF ( clef_files(n) .and. mod(itap,nint(freq_outnmc(n)/phys_tstep))==1) THEN
    62       ! print*,'ini_undefSTD n itap',n,itap
     61    !IF ( clef_files(n) .and. mod(itap,nint(freq_outnmc(n)/phys_tstep))==1) THEN
     62    IF ( mod(itap,nint(freq_outnmc(n)/phys_tstep))==1) THEN
    6363      DO k = 1, nlevstd
    6464        DO i = 1, klon
    6565          tnondef(i, k, n) = 0.
    66           tsumstd(i, k, n) = 0.
    67           usumstd(i, k, n) = 0.
    68           vsumstd(i, k, n) = 0.
    69           wsumstd(i, k, n) = 0.
    70           phisumstd(i, k, n) = 0.
    71           qsumstd(i, k, n) = 0.
    72           rhsumstd(i, k, n) = 0.
    73           uvsumstd(i, k, n) = 0.
    74           vqsumstd(i, k, n) = 0.
    75           vtsumstd(i, k, n) = 0.
    76           wqsumstd(i, k, n) = 0.
    77           vphisumstd(i, k, n) = 0.
    78           wtsumstd(i, k, n) = 0.
    79           u2sumstd(i, k, n) = 0.
    80           v2sumstd(i, k, n) = 0.
    81           t2sumstd(i, k, n) = 0.
    82           o3sumstd(i, k, n) = 0.
    83           o3daysumstd(i, k, n) = 0.
    8466        END DO !i
    8567      END DO !k
     68      IF ( clef_files(n) ) THEN
     69        ! print*,'ini_undefSTD n itap',n,itap
     70        DO k = 1, nlevstd
     71          DO i = 1, klon
     72            tnondef(i, k, n) = 0.
     73            tsumstd(i, k, n) = 0.
     74            usumstd(i, k, n) = 0.
     75            vsumstd(i, k, n) = 0.
     76            wsumstd(i, k, n) = 0.
     77            phisumstd(i, k, n) = 0.
     78            qsumstd(i, k, n) = 0.
     79            rhsumstd(i, k, n) = 0.
     80            uvsumstd(i, k, n) = 0.
     81            vqsumstd(i, k, n) = 0.
     82            vtsumstd(i, k, n) = 0.
     83            wqsumstd(i, k, n) = 0.
     84            vphisumstd(i, k, n) = 0.
     85            wtsumstd(i, k, n) = 0.
     86            u2sumstd(i, k, n) = 0.
     87            v2sumstd(i, k, n) = 0.
     88            t2sumstd(i, k, n) = 0.
     89            o3sumstd(i, k, n) = 0.
     90            o3daysumstd(i, k, n) = 0.
     91          END DO !i
     92        END DO !k
     93      END IF !
    8694    END IF !
    8795  END DO !n
Note: See TracChangeset for help on using the changeset viewer.