Ignore:
Timestamp:
Apr 3, 2023, 10:25:44 AM (15 months ago)
Author:
crisi
Message:

Bug corrections in LMDZiso, especially for water tagging

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmdiso/reevap.F90

    r4143 r4491  
    8585      do ixt=1,ntiso
    8686        call iso_verif_noNaN(xt_seri(ixt,i,k), &
    87      &     'physiq 2417: apres evap tot')
     87     &     'reevap 2417: apres evap tot')
    8888      enddo
    8989      if (iso_eau.gt.0) then
    9090              call iso_verif_egalite_choix( &
    9191     &           xt_seri(iso_eau,i,k),q_seri(i,k), &
    92      &          'physiq 1891+, après reevap totale',errmax,errmaxrel)
     92     &          'reevap 1891+, après reevap totale',errmax,errmaxrel)
    9393              call iso_verif_egalite_choix( &
    9494     &           xtl_seri(iso_eau,i,k),ql_seri(i,k), &
    95      &          'physiq 2209+, après reevap totale',errmax,errmaxrel)
     95     &          'reevap 2209+, après reevap totale',errmax,errmaxrel)
    9696       endif !if (iso_eau.gt.0) then       
    9797      if ((iso_HDO.gt.0).and.(iso_O18.gt.0)) then 
     
    100100     &              xt_seri(iso_HDO,i,k)/q_seri(i,k), &
    101101     &              xt_seri(iso_O18,i,k)/q_seri(i,k), &
    102      &              'physiq 2315: apres reevap totale').eq.1) then
     102     &              'reevap 2315: apres reevap totale').eq.1) then
    103103                  write(*,*) 'i,k,q_seri(i,k)=',i,k,q_seri(i,k)
    104104                  write(*,*) 'd_q_eva(i,k)=',d_q_eva(i,k)
     
    111111#ifdef ISOTRAC     
    112112             call iso_verif_traceur(xt_seri(1,i,k), &
    113      &           'physiq 2165')
     113     &           'reevap 2165a')
    114114             call iso_verif_traceur_pbidouille(xt_seri(1,i,k), &
    115      &           'physiq 2165b')
     115     &           'reevap 2165b')
    116116#endif               
    117117         
     
    146146      do ixt=1,ntiso
    147147      call iso_verif_noNaN(xt_seri(ixt,i,k), &
    148      &     'physiq 2417: apres evap tot')
     148     &     'reevap 2417: apres evap tot')
    149149      enddo
    150150      if (iso_eau.gt.0) then
    151151              call iso_verif_egalite_choix( &
    152152     &           xt_seri(iso_eau,i,k),q_seri(i,k), &
    153      &          'physiq 1891, après réévap totale',errmax,errmaxrel)
     153     &          'reevap 1891, après réévap totale',errmax,errmaxrel)
    154154              call iso_verif_egalite_choix( &
    155155     &           xtl_seri(iso_eau,i,k),ql_seri(i,k), &
    156      &          'physiq 2209, après réévap totale',errmax,errmaxrel)
     156     &          'reevap 2209, après réévap totale',errmax,errmaxrel)
    157157              call iso_verif_egalite_choix( &
    158158     &           xts_seri(iso_eau,i,k),qs_seri(i,k), &
    159      &          'physiq 2209b, après réévap totale',errmax,errmaxrel)
     159     &          'reevap 2209b, après réévap totale',errmax,errmaxrel)
    160160       endif !if (iso_eau.gt.0) then
    161161     
     
    165165     &              xt_seri(iso_HDO,i,k)/q_seri(i,k), &
    166166     &              xt_seri(iso_O18,i,k)/q_seri(i,k), &
    167      &              'physiq 2408: apres reevap totale').eq.1) then
     167     &              'reevap 2408: apres reevap totale').eq.1) then
    168168                  write(*,*) 'i,k,q_seri(i,k)=',i,k,q_seri(i,k)                       
    169169                  stop
     
    173173#ifdef ISOTRAC     
    174174             call iso_verif_traceur(xt_seri(1,i,k), &
    175      &           'physiq 2165')
     175     &           'reevap 2165c')
    176176             call iso_verif_traceur_pbidouille(xt_seri(1,i,k), &
    177      &           'physiq 2165b')
     177     &           'reevap 2165d')
    178178#endif                 
    179179#endif                 
Note: See TracChangeset for help on using the changeset viewer.