temperature

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
sara.oceanography
Posts: 6
Joined: Wed May 30, 2018 3:18 pm
Location: Ghent university

temperature

#1 Post by sara.oceanography » Wed Jun 27, 2018 6:16 am

Hello everybody,
I used the Ecmwf data to create the forcing file, but I got unusual values in the coastal region and the temperature is high in the regions, I have read the previous posts and I understand I should extrapolate the data because the land-sea mask doesn't match with the grid file,
I did some attempts with cdo and ncl to extrapolate them but I did not succeed,
I am very very new to modeling, I would be grateful if you could tell me the easiest way to do that.

I am looking forward to any suggestions.
Best Regards,
Sara

User avatar
kate
Posts: 3700
Joined: Wed Jul 02, 2003 5:29 pm
Location: IMS/UAF, USA

Re: temperature

#2 Post by kate » Wed Jun 27, 2018 5:37 pm

My previous reply on this included the link to sosie. An example of using it (for MERRA):

Code: Select all

#!/bin/bash

SOSIE=/home/raphael/TOOLS/sosie-rutgers/bin
dirin=/t3/workdir/raphael/MERRA/2014
dirout=/t3/workdir/raphael/MERRA/2014/drowned

ybeg=2014
yend=2014

# build the list of variables to proceed
list_var=$( ls $dirin | grep $ybeg | sed -e "s/MERRA_//" -e "s/_$ybeg.nc//" | awk '{ print $1 }' )

for year in $( seq $ybeg $yend ) ; do

    for var in $list_var ; do

        filein=MERRA_${var}_$year.nc
        timevar=$( ncdump -h $dirin/$filein | grep time | head -1 | awk '{ print $1 }' )
        ncvar=$( ncdump -h $dirin/$filein | grep float | tail -1 | sed -e "s/(/ /" | awk '{ print $2 }' )

        echo Processing variable $ncvar in file $filein

        $SOSIE/mask_drown_field_roms.x -D -i $dirin/$filein -v $ncvar -t $timevar \
                                     -m $HOME/STORAGE/MERRA/mask_MERRA.nc -o $dirout/drowned_$filein



    done

done

Post Reply