Changeset 543


Ignore:
Timestamp:
Jun 29, 2004, 6:06:31 PM (20 years ago)
Author:
lmdzadmin
Message:

On met des valeurs pour les traceurs par defaut si le fichier traceur.def
n'existe pas
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/dyn3d/iniadvtrac.F

    r524 r543  
    125125     s     iostat=ierr)
    126126      if(ierr.eq.0) then
    127        print*,'ouverture de traceur.def ok'
    128        read(90,*) nq
    129        print*,'nombre de traceurs ',nq
    130        if (nq.gt.nqmx) then
    131        print*,'nombre de traceurs trop important'
    132        print*,'verifier traceur.def'
    133        stop
    134        endif
     127        print*,'ouverture de traceur.def ok'
     128        read(90,*) nq
     129        print*,'nombre de traceurs ',nq
     130        if (nq.gt.nqmx) then
     131          print*,'nombre de traceurs trop important'
     132          print*,'verifier traceur.def'
     133          stop
     134        endif
    135135C
    136        do iq=1,nq
    137          read(90,999) hadv(iq),vadv(iq),tnom(iq)
    138        end do
    139        close(90) 
    140        PRINT*,'lecture de traceur.def :'   
    141        do iq=1,nq
    142          write(*,*) hadv(iq),vadv(iq),tnom(iq)
    143        end do       
     136        do iq=1,nq
     137          read(90,999) hadv(iq),vadv(iq),tnom(iq)
     138        end do
     139        close(90) 
     140        PRINT*,'lecture de traceur.def :'   
     141        do iq=1,nq
     142          write(*,*) hadv(iq),vadv(iq),tnom(iq)
     143        end do       
    144144      else
    145        print*,'pb ouverture traceur.def'
    146        print*,'ATTENTION on prend des valeurs par defaut'
    147       endif
     145        print*,'pb ouverture traceur.def'
     146        print*,'ATTENTION on prend des valeurs par defaut'
     147        nq = 4
     148        hadv(1) = 14
     149        vadv(1) = 14
     150        tnom(1) = 'H2Ov'
     151        hadv(2) = 10
     152        vadv(2) = 10
     153        tnom(2) = 'H2Ol'
     154        hadv(3) = 10
     155        vadv(3) = 10
     156        tnom(3) = 'RN'
     157        hadv(4) = 10
     158        vadv(4) = 10
     159        tnom(4) = 'PB'
     160      ENDIF
     161      PRINT*,'Valeur de traceur.def :'
     162      do iq=1,nq
     163        write(*,*) hadv(iq),vadv(iq),tnom(iq)
     164      end do       
     165
    148166#endif
    149167c a partir du nom court du traceur et du schema d'advection au detemine le nom long.
Note: See TracChangeset for help on using the changeset viewer.