Ignore:
Timestamp:
03/23/18 10:55:23 (6 years ago)
Author:
omamce
Message:

O.M. : correct OIceFrc variable name.
More documentation

File:
1 edited

Legend:

Unmodified
Added
Removed
  • TOOLS/MOSAIX/CreateRestartOce4Oasis.bash

    r3665 r3686  
    104104## Extract variables 
    105105## =========================================================================== 
    106 ncks --overwrite --fl_fmt=${FL_FMT} --history -d time_counter,0,0 -v ${IceVar} ${IceFile} sstoce_fields.nc 
    107 ncks --overwrite --fl_fmt=${FL_FMT} --history -d time_counter,0,0 -v tos       ${OceFile} oce_fields.nc 
     106ncks --overwrite --fl_fmt=${FL_FMT} --history -d time_counter,0,0 --variable ${IceVar} ${IceFile} sstoce_fields.nc 
     107ncks --overwrite --fl_fmt=${FL_FMT} --history -d time_counter,0,0 --variable tos       ${OceFile} oce_fields.nc 
    108108ncks --append --fl_fmt=${FL_FMT} --history oce_fields.nc   sstoce_fields.nc 
    109 ncwa --overwrite --fl_fmt=${FL_FMT} --history -a time_counter sstoce_fields.nc sstoce_fields_notime.nc # Remove time dimension 
     109ncwa --overwrite --fl_fmt=${FL_FMT} --history --average time_counter sstoce_fields.nc sstoce_fields_notime.nc # Remove time dimension 
    110110ncatted --history --attribute history,global,d,c,""  sstoce_fields_notime.nc           # Clean attributes 
    111111 
     
    126126 
    127127cat <<EOF > create_sstoce.nco 
    128 OIceFrac       = double ( ${IceVar}(:,:) ) ;  
     128OIceFrc        = double ( ${IceVar}(:,:) ) ;  
    129129//tos.simple_fill_miss(tos)               ; 
    130 O_SSTSST[y,x]  = double ( (tos   (:,:) + 273.15d) * (1.0d-OIceFrac) ) ; 
    131 //O_AlbIce[y,x]  = double ( ialb   (:,:)            * OIceFrac(:,:)   ) ; 
    132 //O_TepIce[y,x]  = double ( (tsice (:,:) + 273.15d) * OIceFrac(:,:)   ) ; 
    133 O_AlbIce[y,x]  = double (   0.65d            * OIceFrac(:,:)   ) ; 
    134 O_TepIce[y,x]  = double ( (-20.0d + 273.15d) * OIceFrac(:,:)   ) ; 
     130O_SSTSST[y,x]  = double ( (tos   (:,:) + 273.15d) * (1.0d-OIceFrc) ) ; 
     131//O_AlbIce[y,x]  = double ( ialb   (:,:)            * OIceFrc(:,:)   ) ; 
     132//O_TepIce[y,x]  = double ( (tsice (:,:) + 273.15d) * OIceFrc(:,:)   ) ; 
     133O_AlbIce[y,x]  = double (   0.65d            * OIceFrc(:,:)   ) ; 
     134O_TepIce[y,x]  = double ( (-20.0d + 273.15d) * OIceFrc(:,:)   ) ; 
    135135O_OCurx1[y,x]  = 0.0d ; 
    136136O_OCury1[y,x]  = 0.0d ; 
     
    138138EOF 
    139139 
    140 ncap2 --overwrite --fl_fmt=${FL_FMT} --history -S create_sstoce.nco sstoce_fields_notime.nc tmp_sstoc.nc 
     140ncap2 --overwrite --fl_fmt=${FL_FMT} --history --script-file create_sstoce.nco sstoce_fields_notime.nc tmp_sstoc.nc 
    141141 
    142 ncks --fl_fmt=${FL_FMT} --overwrite --history -v OIceFrac,O_SSTSST,O_AlbIce,O_TepIce,O_OCurx1,O_OCury1,O_OCurz1 tmp_sstoc.nc sstoc.nc 
     142ncks --fl_fmt=${FL_FMT} --overwrite --history --variable OIceFrc,O_SSTSST,O_AlbIce,O_TepIce,O_OCurx1,O_OCury1,O_OCurz1 tmp_sstoc.nc sstoc.nc 
    143143 
    144144ncatted --history --attribute comment,O_SSTSST,o,c,"SST weighted by fraction of open ocean"          sstoc.nc 
Note: See TracChangeset for help on using the changeset viewer.