Changeset 4081 for TOOLS/MOSAIX


Ignore:
Timestamp:
10/19/18 13:40:07 (6 years ago)
Author:
omamce
Message:

O.M. : add output fields for diagnostics

Location:
TOOLS/MOSAIX
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • TOOLS/MOSAIX/CreateWeightsMask.bash

    r4037 r4081  
    5353# ============== 
    5454#OCE=ORCA2.3 
    55 #OCE=eORCA1.2 
    56 OCE=ORCA025 
    57  
    58 ATM=ICO30 
    59 #ATM=ICO40 
     55OCE=eORCA1.2 
     56#OCE=ORCA025 
     57 
     58#ATM=ICO30 
     59ATM=ICO40 
    6060#ATM=ICO450 
    6161#ATM=LMD9695 
    … …  
    6666# =========================== 
    6767ListOCEGRID="T U V" 
    68 ListOrder="1st" # 2nd" 
    69 ListNormalize="false" # true" 
     68ListOrder="1st 2nd" 
     69ListNormalize="false true" 
    7070 
    7171## =========================================================================== 
    … …  
    217217 
    218218time /usr/bin/time ${MPIRUN} ./interpol.exe --mask_src=yes --mask_dst=no 
    219  
    220219## 
    221220## Correct spurious values (extremes) 
    … …  
    237236cp dia_t${oce}_to_t${atm}_1storder_false.nc  dia_t${oce}_to_t${atm}_1storder_false_mask.nc 
    238237ncks --alphabetize --history --overwrite --variable OceFrac dia_t${oce}_to_t${atm}_1storder_false_mask.nc  ${ATM}_grid_maskFrom_${OCE}.nc 
    239 ncks --alphabetize --history --append    --variable aire    atm_grid.nc                                    ${ATM}_grid_maskFrom_${OCE}.nc 
    240 [[ ${atm} = *ico* ]] && ncks --alphabetize --history --append --variable bounds_lon,bounds_lat atm_grid.nc ${ATM}_grid_maskFrom_${OCE}.nc 
    241238 
    242239cat <<EOF > creation_masque.nco 
    … …  
    244241where (OceFrac <= 0.0 )  OceFrac=0 ; 
    245242EOF 
     243 
    246244ncap2 --history --overwrite --script-file creation_masque.nco dia_t${oce}_to_t${atm}_1storder_false_mask.nc tmp_dia_t${oce}_to_t${atm}_1storder_false_mask.nc ; mv tmp_dia_t${oce}_to_t${atm}_1storder_false_mask.nc dia_t${oce}_to_t${atm}_1storder_false_mask.nc 
    247245ncrename --history --variable OceFrac,OceMask dia_t${oce}_to_t${atm}_1storder_false_mask.nc 
    248246 
    249 ncks --alphabetize --history --append --variable OceMask dia_t${oce}_to_t${atm}_1storder_false_mask.nc   ${ATM}_grid_maskFrom_${OCE}.nc 
     247ncks --overwrite --history --variable OceMask dia_t${oce}_to_t${atm}_1storder_false_mask.nc tmp_OceMask.nc 
     248ncks --history --append tmp_OceMask.nc ${ATM}_grid_maskFrom_${OCE}.nc 
    250249rm dia_t${oce}_to_t${atm}_1storder_false_mask.nc 
     250 
     251# Change dimension names, and some attributes accordingly 
     252if [[ $(ncdump -h ${ATM}_grid_maskFrom_${OCE}.nc | grep domain_dst | wc -l) -gt 0 ]] ; then 
     253    case ${atm} in 
     254        ( *ico* )  
     255        ncrename --history --dimension cell_domain_dst,cell       ${ATM}_grid_maskFrom_${OCE}.nc 
     256        ;; 
     257        ( *lmd* )  
     258        ncrename --history --dimension x_domain_dst,x             ${ATM}_grid_maskFrom_${OCE}.nc 
     259        ncrename --history --dimension y_domain_dst,y             ${ATM}_grid_maskFrom_${OCE}.nc 
     260        ;; 
     261    esac 
     262    ncrename --history --dimension nvertex_domain_dst,nvertex        ${ATM}_grid_maskFrom_${OCE}.nc 
     263    ncrename --history --variable lat_domain_dst,lat                 ${ATM}_grid_maskFrom_${OCE}.nc 
     264    ncrename --history --variable lon_domain_dst,lon                 ${ATM}_grid_maskFrom_${OCE}.nc 
     265    ncrename --history --variable bounds_lat_domain_dst,bounds_lat   ${ATM}_grid_maskFrom_${OCE}.nc 
     266    ncrename --history --variable bounds_lon_domain_dst,bounds_lon   ${ATM}_grid_maskFrom_${OCE}.nc 
     267    ncatted  --history --attribute bounds,lat,m,c,"bounds_lat"       ${ATM}_grid_maskFrom_${OCE}.nc 
     268    ncatted  --history --attribute bounds,lon,m,c,"bounds_lon"       ${ATM}_grid_maskFrom_${OCE}.nc 
     269    ncatted  --history --attribute coordinates,OceFrac,m,c,"lat lon" ${ATM}_grid_maskFrom_${OCE}.nc 
     270    ncatted  --history --attribute coordinates,OceMask,m,c,"lat lon" ${ATM}_grid_maskFrom_${OCE}.nc 
     271fi 
     272 
     273ncks --history --alphabetize  --append    --variable aire    atm_grid.nc                                   ${ATM}_grid_maskFrom_${OCE}.nc 
     274[[ ${atm} = *ico* ]] && ncks --alphabetize --history --append --variable bounds_lon,bounds_lat atm_grid.nc ${ATM}_grid_maskFrom_${OCE}.nc 
     275 
    251276 
    252277## 
    … …  
    347372EOF 
    348373ncap2 --history --overwrite --script-file coastal.nco ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_coastal_maskFrom_${OCE}.nc 
     374cp ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}_n07.nc 
    349375ncks --history --append --variable AtmCoastal ${ATM}_coastal_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}.nc 
     376cp ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}_n08.nc 
    350377rm ${ATM}_coastal_maskFrom_${OCE}.nc 
    351378ncks --alphabetize --history --overwrite --fl_fmt=${FMT_OASIS} ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}_${FMT_OASIS}.nc 
     379cp ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}_n09.nc 
     380 
    352381## 
    353382## ATM towards NEMO points - runoff 
    … …  
    396425#ncatted --history -a coordinates,aire,m,c,"lat lon" ${ATM}_grid_maskFrom_${OCE}.nc 
    397426ncap2 --overwrite --history --script 'defdim("time_counter",1) ;' ${ATM}_grid_maskFrom_${OCE}.nc tmp_${ATM}_grid_maskFrom_${OCE}.nc ; mv tmp_${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}.nc 
    398 ncdump -h ${ATM}_grid_maskFrom_${OCE}.nc 
     427cp ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}_n10.nc 
    399428 
    400429ncks --alphabetize --history --overwrite --fl_fmt=${FMT_OASIS} ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}_${FMT_OASIS}.nc 
    401 ncdump -h ${ATM}_grid_maskFrom_${OCE}_${FMT_OASIS}.nc 
    402430 
    403431if [[ ${atm} = ico ]] ; then 
    … …  
    408436EOF 
    409437    ncap2 --overwrite --history --script-file add_time.nco ${ATM}_grid_maskFrom_${OCE}.nc tmp_${ATM}_grid_maskFrom_${OCE}.nc ; mv tmp_${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}.nc 
     438    cp ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}_n11.nc 
    410439    ncatted --history \ 
    411             -a coordinates,OceFrac,c,c,"time_counter cell"    \ 
    412             -a coordinates,OceMask,c,c,"time_counter cell"    \ 
    413             -a coordinates,AtmCoastal,c,c,"time_counter cell" \ 
    414             -a coordinates,aire,c,c,"cell" \ 
     440            -a coordinates,OceFrac,m,c,"time_counter lat lon"    \ 
     441            -a coordinates,OceMask,m,c,"time_counter lat lon"    \ 
     442            -a coordinates,AtmCoastal,m,c,"time_counter lat lon" \ 
     443            -a coordinates,aire,c,c,"lat lon" \ 
    415444            ${ATM}_grid_maskFrom_${OCE}.nc 
     445    cp ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}_n12.nc 
    416446fi 
    417447ncks --alphabetize --history --overwrite --fl_fmt=${FMT_OASIS} ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}_${FMT_OASIS}.nc 
    … …  
    426456    ncdump -h tmp_${ATM}_grid_maskFrom_${OCE}.nc 
    427457    mv tmp_${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}.nc 
     458    cp ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}_n13.nc 
    428459     
    429460    ncatted --history \ 
    … …  
    433464            -a coordinates,aire,m,c,"lat lon" \ 
    434465            ${ATM}_grid_maskFrom_${OCE}.nc 
     466    cp ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}_n14.nc 
    435467fi 
    436468ncks --alphabetize --history --overwrite --mk_rec time_counter   ${ATM}_grid_maskFrom_${OCE}.nc tmp_${ATM}_grid_maskFrom_${OCE}.nc ; mv tmp_${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}.nc 
     469cp ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}_n15.nc 
    437470ncks --alphabetize --history --overwrite --fl_fmt=${FMT_OASIS} ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}_${FMT_OASIS}.nc 
     471cp ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}_n16.nc 
    438472 
    439473## 
    … …  
    480514            ${InFile} 
    481515done 
    482  
    483  
     516cp ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}_n20.nc 
    484517## 
    485518## Update and complete weights file to fit OASIS requested format 

  • TOOLS/MOSAIX/iodef_atm_to_oce.xml

    r3665 r4081  
    8282        <variable name="SVN_Id"        type="string">$Id$</variable> 
    8383        <variable name="SVN_HeadURL"   type="string">$HeadURL$</variable> 
     84 
     85        <field field_ref="field01_src" operation="once" prec="8" /> 
     86        <field field_ref="field02_src" operation="once" prec="8" /> 
     87        <field field_ref="field03_src" operation="once" prec="8" /> 
     88        <field field_ref="field04_src" operation="once" prec="8" /> 
     89        <field field_ref="field05_src" operation="once" prec="8" /> 
     90        <field field_ref="field06_src" operation="once" prec="8" /> 
    8491         
    8592        <field name="field01_dst" field_ref="field01_src" domain_ref="domain_dst" operation="once" prec="8" /> 

  • TOOLS/MOSAIX/iodef_oce_to_atm.xml

    r3665 r4081  
    8282        <variable name="SVN_Id"        type="string">$Id$</variable> 
    8383        <variable name="SVN_HeadURL"   type="string">$HeadURL$</variable> 
     84 
     85        <field field_ref="mask_src"     operation="once" prec="8" /> 
     86        <field field_ref="field01_src"  operation="once" prec="8" /> 
     87        <field field_ref="field02_src"  operation="once" prec="8" /> 
     88        <field field_ref="field03_src"  operation="once" prec="8" /> 
     89        <field field_ref="field04_src"  operation="once" prec="8" /> 
     90        <field field_ref="field05_src"  operation="once" prec="8" /> 
     91        <field field_ref="field06_src"  operation="once" prec="8" /> 
    8492         
    8593        <field name="OceFrac"     field_ref="mask_src"    domain_ref="domain_dst" operation="once" prec="8" /> 
Note: See TracChangeset for help on using the changeset viewer.