Changeset 4038 for LMDZ6/trunk/libf


Ignore:
Timestamp:
Dec 1, 2021, 8:27:42 AM (3 years ago)
Author:
crisi
Message:

correction dans check_isotopes + suppression d'un call check_isotopes problematique dans advtrac

Location:
LMDZ6/trunk/libf/dyn3dmem
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/dyn3dmem/advtrac_loc.F

    r2622 r4038  
    156156     *                        pk,teta )
    157157
    158           !write(*,*) 'advtrac 162: apres appel vlspltgen_loc'
    159       if (ok_iso_verif) then
    160            call check_isotopes(q,ijb_u,ije_u,'advtrac 162')
    161       endif !if (ok_iso_verif) then
    162 
    163158#ifdef DEBUG_IO     
    164159          do j=1,nqtot
  • LMDZ6/trunk/libf/dyn3dmem/check_isotopes_loc.F

    r2281 r4038  
    2424        parameter (errmaxrel=1e-3)
    2525        parameter (qmin=1e-11)
    26         parameter (deltaDmax=200.0,deltaDmin=-999.9)
     26        parameter (deltaDmax=1000.0,deltaDmin=-999.0)
    2727        parameter (ridicule=1e-12)
    2828        real deltaD
     
    3030        if (ok_isotopes) then
    3131
    32         !write(*,*) 'check_isotopes 31: err_msg=',err_msg
     32!        write(*,*) 'check_isotopes 31: err_msg=',err_msg
    3333        ! verifier que rien n'est NaN
    3434        do ixt=1,ntraciso
     
    4545                  write(*,*) 'q,i,k,iq=',q(i,k,iq),i,k,iq
    4646                  write(*,*) 'borne=',borne
    47                   stop
     47                  call abort_gcm('check_isotopes_loc','plantage iso',0)
    4848                endif  !if ((x(ixt,i,j).gt.-borne).and.
    4949              enddo !DO i = ijb,ije
     
    5353        enddo !do ixt=1,ntraciso
    5454
    55         !write(*,*) 'check_isotopes 52'
     55!        write(*,*) 'check_isotopes 52'
    5656        ! verifier que l'eau normale est OK
    5757        if (use_iso(1)) then
     
    7171                  write(*,*) 'q,iq,i,k=',q(i,k,iq),iq,i,k
    7272                  write(*,*) 'q(i,k,phase)=',q(i,k,phase)
    73                   stop
     73                  call abort_gcm('check_isotopes_loc','plantage iso',0)
    7474              endif !if ((abs(q(i,k,phase)-q(i,k,iq)).gt.errmax).and.
    7575              ! bidouille pour éviter divergence:
     
    8181        endif !if (use_iso(1)) then
    8282       
    83         !write(*,*) 'check_isotopes 78'
     83!        write(*,*) 'check_isotopes 78'
    8484        ! verifier que HDO est raisonable
    8585        if (use_iso(2)) then
     
    9999                  write(*,*) 'q=',q(i,k,:)
    100100                  write(*,*) 'deltaD=',deltaD
    101                   stop
     101                  call abort_gcm('check_isotopes_loc','plantage iso',0)
    102102             endif !if ((deltaD.gt.deltaDmax).or.(deltaD.lt.deltaDmin)) then
    103103            endif !if (q(i,k,iq).gt.qmin) then
     
    108108        endif !if (use_iso(2)) then
    109109
    110         !write(*,*) 'check_isotopes 103'
     110!        write(*,*) 'check_isotopes 103'
    111111        ! verifier que O18 est raisonable
    112112        if (use_iso(3)) then
     
    126126                  write(*,*) 'xt=',q(i,k,:)
    127127                  write(*,*) 'deltaO18=',deltaD
    128                   stop
     128                  call abort_gcm('check_isotopes_loc','plantage iso',0)
    129129             endif !if ((deltaD.gt.deltaDmax).or.(deltaD.lt.deltaDmin)) then
    130130            endif !if (q(i,k,iq).gt.qmin) then
     
    136136
    137137
    138         !write(*,*) 'check_isotopes 129'
     138!        write(*,*) 'check_isotopes 129'
    139139        if (ok_isotrac) then
    140140
     
    158158                  write(*,*) 'ixt,ieau=',ixt,ieau
    159159                  write(*,*) 'q,iq,i,k,=',q(i,k,iq),iq,i,k
    160                   write(*,*) 'deltaD=',deltaD
    161                   stop
     160                  write(*,*) 'deltaD=',deltaD                 
     161                  call abort_gcm('check_isotopes_loc','plantage iso',0)
    162162                 endif !if ((deltaD.gt.deltaDmax).or.
    163163                endif !if (q(i,k,iq).gt.qmin) then
     
    190190                  write(*,*) 'i,k,=',i,k
    191191                  write(*,*) 'q(i,k,:)=',q(i,k,:)
    192                   stop
     192                  call abort_gcm('check_isotopes_loc','plantage iso',0)
    193193                 endif !if ((abs(q(i,k,phase)-q(i,k,iq)).gt.errmax).and.
    194194                 
     
    209209
    210210        endif ! if (ok_isotopes)
    211         !write(*,*) 'check_isotopes 198'
     211!        write(*,*) 'check_isotopes 198'
    212212       
    213213        end
Note: See TracChangeset for help on using the changeset viewer.