Ignore:
Timestamp:
Sep 5, 2014, 5:42:24 PM (10 years ago)
Author:
Laurent Fairhead
Message:

Suite (et fin?) des modifications pour permettre un controle 'tout xml'
des fichiers de sorties XIOS.
Un nouveau paramètre logique est introduit, ok_all_xml, false par défaut, et lu
dans le run.def qui permet de faire du 'tout xml'

Follow-up modifications to ensure total xlm control over the output files
from XIOS.
A new logical parameter, ok_all_wml, is introduced. False by default, it is
read from the run.def file and, if true, will give over control to the
XIOS xml files

Location:
LMDZ5/trunk/libf/phylmd
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phylmd/clesphys.h

    r2110 r2114  
    1212       LOGICAL cycle_diurne,soil_model,new_oliq,ok_orodr,ok_orolf
    1313       LOGICAL ok_limitvrai
     14       LOGICAL ok_all_xml
    1415       INTEGER nbapp_rad, iflag_con,iflag_ener_conserv
    1516       REAL co2_ppm, co2_ppm0, solaire
    … …  
    115116     &     , iflag_rrtm, ok_strato,ok_hines                             &
    116117     &     , iflag_ice_thermo, ok_gwd_rando, NSW                        &
    117      &     , ok_conserv_q
     118     &     , ok_conserv_q, ok_all_xml
    118119     
    119120       save /clesphys/

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

    r2110 r2114  
    8686    logical,SAVE              :: ok_snow_omp
    8787    logical,SAVE              :: ok_newmicro_omp
     88    logical,SAVE              :: ok_all_xml_omp
    8889    logical,SAVE        :: ok_journe_omp, ok_mensuel_omp, ok_instan_omp, ok_hf_omp       
    8990    logical,SAVE        :: ok_LES_omp   
    … …  
    243244       CALL abort_gcm('conf_phys','t_coupl = 0.',1)
    244245    END IF
     246
     247    !
     248    !Config Key  = ok_all_xml
     249    !Config Desc = utiliser les xml pourles définitions des champs pour xios
     250    !Config Def  = .false.
     251    !Config Help =
     252    !
     253    ok_all_xml_omp = .false.
     254    call getin('ok_all_xml', ok_all_xml_omp)
     255    !
    245256
    246257    !
    … …  
    18491860    ! Martin
    18501861
     1862    ok_all_xml = ok_all_xml_omp
    18511863    ok_newmicro = ok_newmicro_omp
    18521864    ok_journe = ok_journe_omp
    … …  
    19781990    write(lunout,*)' Config veget = ', ok_veget,type_veget
    19791991    write(lunout,*)' Snow model SISVAT : ok_snow = ', ok_snow
     1992    write(lunout,*)' Config xml pour XIOS : ok_all_xml = ', ok_all_xml
    19801993    write(lunout,*)' Sortie journaliere = ', ok_journe
    19811994    write(lunout,*)' Sortie haute frequence = ', ok_hf

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

    r2097 r2114  
    204204  IMPLICIT NONE
    205205  include 'dimensions.h'
     206  include 'clesphys.h'
    206207   
    207208    character*(*), INTENT(IN) :: name
    … …  
    227228    IF((.NOT. is_using_mpi) .OR. is_mpi_root) THEN
    228229        ! ug Création du fichier
     230      IF (.not. ok_all_xml) THEN
    229231        CALL wxios_add_file(name, ffreq, lev)
     232      ENDIF
    230233    END IF
    231234#endif
    … …  
    603606
    604607#ifdef CPP_XIOS
    605       IF ( var%flag(iff)<=lev_files(iff) ) THEN
    606         CALL wxios_add_field_to_file(var%name, 2, iff, phys_out_filenames(iff), &
    607         var%description, var%unit, var%flag(iff), typeecrit)
    608         IF (prt_level >= 10) THEN
    609           WRITE(lunout,*) 'histdef2d: call wxios_add_field_to_file var%name iff: ', &
    610                           trim(var%name),iff
     608      IF (.not. ok_all_xml) THEN
     609        IF ( var%flag(iff)<=lev_files(iff) ) THEN
     610          CALL wxios_add_field_to_file(var%name, 2, iff, phys_out_filenames(iff), &
     611          var%description, var%unit, var%flag(iff), typeecrit)
     612          IF (prt_level >= 10) THEN
     613            WRITE(lunout,*) 'histdef2d: call wxios_add_field_to_file var%name iff: ', &
     614                            trim(var%name),iff
     615          ENDIF
    611616        ENDIF
    612617      ENDIF
    … …  
    688693
    689694#ifdef CPP_XIOS
    690       IF ( var%flag(iff)<=lev_files(iff) ) THEN
    691         CALL wxios_add_field_to_file(var%name, 3, iff, phys_out_filenames(iff), &
    692         var%description, var%unit, var%flag(iff), typeecrit)
    693         IF (prt_level >= 10) THEN
    694           WRITE(lunout,*) 'histdef3d: call wxios_add_field_to_file var%name iff: ', &
    695                           trim(var%name),iff
    696         ENDIF
    697       ENDIF
     695       IF (.not. ok_all_xml) THEN
     696         IF ( var%flag(iff)<=lev_files(iff) ) THEN
     697         CALL wxios_add_field_to_file(var%name, 3, iff, phys_out_filenames(iff), &
     698         var%description, var%unit, var%flag(iff), typeecrit)
     699           IF (prt_level >= 10) THEN
     700             WRITE(lunout,*) 'histdef3d: call wxios_add_field_to_file var%name iff: ', &
     701                             trim(var%name),iff
     702           ENDIF
     703         ENDIF
     704       ENDIF
    698705#endif
    699706#ifndef CPP_IOIPSL_NO_OUTPUT
    … …  
    890897  INCLUDE 'dimensions.h'
    891898  INCLUDE 'iniprint.h'
     899  include 'clesphys.h'
    892900
    893901    TYPE(ctrl_out), INTENT(IN) :: var
    … …  
    922930!$OMP MASTER
    923931      !Si phase de définition.... on définit
     932      IF (.not. ok_all_xml) THEN
     933      if (prt_level >= 10) then
     934      write(lunout,*)"histwrite2d_phy: .not.vars_defined ; time to define ", &
     935                     trim(var%name)
     936      endif
    924937      DO iff=iff_beg, iff_end
    925938         IF (clef_files(iff)) THEN
    … …  
    927940         ENDIF
    928941      ENDDO
     942      ENDIF
    929943!$OMP END MASTER
    930944  ELSE
    … …  
    944958! La boucle sur les fichiers:
    945959      firstx=.true.
    946       DO iff=iff_beg, iff_end
     960
     961      IF (ok_all_xml) THEN
     962#ifdef CPP_XIOS
     963          if (prt_level >= 10) then
     964             write(lunout,*)'Dans iophy histwrite2D,var%name ',&
     965                             trim(var%name)                       
     966          endif
     967          CALL xios_send_field(var%name, Field2d)
     968          if (prt_level >= 10) then
     969             write(lunout,*)'Dans iophy histwrite2D,var%name apres xios_send ',&
     970                             trim(var%name)                       
     971          endif
     972#else
     973        CALL abort_gcm ('iophy','cannot have ok_all_xml = .T. without CPP_XIOS defined' ,1)
     974#endif
     975      ELSE 
     976        DO iff=iff_beg, iff_end
    947977            IF (var%flag(iff) <= lev_files(iff) .AND. clef_files(iff)) THEN
    948978
    … …  
    10011031                deallocate(fieldok)
    10021032            ENDIF !levfiles
    1003       ENDDO ! of DO iff=iff_beg, iff_end
     1033        ENDDO ! of DO iff=iff_beg, iff_end
     1034      ENDIF
    10041035!$OMP END MASTER   
    10051036  ENDIF ! vars_defined
    … …  
    10261057  INCLUDE 'dimensions.h'
    10271058  INCLUDE 'iniprint.h'
     1059  include 'clesphys.h'
    10281060
    10291061    TYPE(ctrl_out), INTENT(IN) :: var
    … …  
    10781110! BOUCLE SUR LES FICHIERS
    10791111     firstx=.true.
     1112
     1113      IF (ok_all_xml) THEN
     1114#ifdef CPP_XIOS
     1115          if (prt_level >= 10) then
     1116             write(lunout,*)'Dans iophy histwrite3D,var%name ',&
     1117                             trim(var%name)                       
     1118          endif
     1119          CALL xios_send_field(var%name, Field3d)
     1120#else
     1121        CALL abort_gcm ('iophy','cannot have ok_all_xml = .T. without CPP_XIOS defined' ,1)
     1122#endif
     1123      ELSE 
     1124
     1125
    10801126     DO iff=iff_beg, iff_end
    10811127            IF (var%flag(iff) <= lev_files(iff) .AND. clef_files(iff)) THEN
    … …  
    11361182            ENDIF
    11371183      ENDDO
     1184      ENDIF
    11381185!$OMP END MASTER   
    11391186  ENDIF ! vars_defined

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

    r2110 r2114  
    236236    type_ecri_files(:)=type_ecri(:)
    237237
     238!    if (ok_all_xml) phys_out_filelevels = 999
     239
    238240    WRITE(lunout,*)'phys_out_lonmin=',phys_out_lonmin
    239241    WRITE(lunout,*)'phys_out_lonmax=',phys_out_lonmax
    … …  
    305307       zoutm(iff) = ecrit_files(iff) ! Frequence ou l on ecrit en seconde
    306308
    307     IF (clef_files(iff)) THEN
    308309
    309310#ifdef CPP_XIOS
    310311!!! Ouverture de chaque fichier XIOS !!!!!!!!!!!
    311     if (prt_level >= 10) then
    312       print*,'phys_output_open: call wxios_add_file with phys_out_filenames(iff)=',trim(phys_out_filenames(iff))                                                                       
    313     endif
    314     CALL wxios_add_file(phys_out_filenames(iff),chtimestep(iff),lev_files(iff)) 
     312    IF (.not. ok_all_xml) then
     313      if (prt_level >= 10) then
     314        print*,'phys_output_open: call wxios_add_file with phys_out_filenames(iff)=',trim(phys_out_filenames(iff))   
     315      endif
     316      CALL wxios_add_file(phys_out_filenames(iff),chtimestep(iff),lev_files(iff)) 
     317    ENDIF
    315318
    316319!!! Declaration des axes verticaux de chaque fichier:
    … …  
    334337#endif
    335338
     339        IF (clef_files(iff)) THEN
    336340!!!!!!!!!!!!!!!!! Traitement dans le cas ou l'on veut stocker sur un domaine limite !!
    337341!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    … …  
    404408#endif
    405409
     410     ENDIF ! clef_files
     411
    406412!CR: ajout d'une variable eau
    407413!      IF (nqtot>=3) THEN
     414
    408415       IF (nqtot>=nqo+1) THEN
    409416!            DO iq=3,nqtot 
    … …  
    413420                           tname(iiq),'Tracer '//ttext(iiq), "-",  &
    414421                           (/ '', '', '', '', '', '', '', '', '' /))
    415 
    416             o_dtr_vdf(iq-nqo) = ctrl_out((/ 5, 7, 7, 7, 10, 10, 11, 11, 11 /), &
     422            o_dtr_vdf(iq-nqo) = ctrl_out((/ 4, 7, 7, 7, 10, 10, 11, 11, 11 /), &
    417423                              'd'//trim(tname(iq))//'_vdf',  &
    418424                              'Tendance tracer '//ttext(iiq), "-" , &
    … …  
    491497      ENDIF
    492498
    493      ENDIF ! clef_files
    494499
    495500   ENDDO !  iff

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

    r2103 r2114  
    240240    USE phys_cal_mod, only : mth_len
    241241
     242
    242243    IMPLICIT NONE
     244
    243245
    244246    INCLUDE "temps.h"
    … …  
    12931295       ENDDO !nfiles
    12941296!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    1295 !       IF (nqtot.GE.3) THEN
    1296 !          DO iq=3,nqtot
    1297           IF (nqtot.GE.nqo+1) THEN
     1297        IF (nqtot.GE.nqo+1) THEN
    12981298            DO iq=nqo+1,nqtot
     1299
    12991300             CALL histwrite_phy(o_trac(iq-nqo), qx(:,:,iq))
    13001301             CALL histwrite_phy(o_dtr_vdf(iq-nqo),d_tr_cl(:,:,iq-nqo))
Note: See TracChangeset for help on using the changeset viewer.