Memo about IsoMATSIRO †

↑

How to run test IsoMATSIRO experiment †

See http://hydro.iis.u-tokyo.ac.jp/~kei/IsoMAT_RN.html for detail (in Japanese).

If you are interested in using IsoMATSIRO, please ask me to give you the address for the source code. IsoMATSIRO is free only for research purpose.

↑

Installation (vertical 1D simulation of IPILPS) †

  1. Download the source code from designated address and extract it. ("agcm5.6" directory is created) 
    $ wget http://xxxx.xx.xxxx/xxxx.tar.gz
$ tar xvfz xxxx.tar.gz
  2. cd agcm5.6 ; setenv AGCMDIR $PWD 
    $ cd agcm5.6
$ setenv AGCMDIR $PWD
  3. Edit Mkinclude.single and link as Mkinclude. 
    $ vi Mkinclude.sigle (select PROJECT "offimat", SYSTEM "Altix", RTYPE "single".)
$ ln -sf Mkinclude.single Mkinclude
  4. cd $AGCMDIR/src/proj/offimat and linke Mkinclude.1D as Mkinclude 
    $ cd src/proj/offimat
$ ln -sf Mkinclude.1D Mkinclude
  5. cd $AGCMDIR/src and make 
    $ cd ../../
$ make dirs
$ make clean
$ make lib ($AGCMDIR/lib/Altix/libagcm5v00ml1soffimat.single.a is created)
$ make gcm ($AGCMDIR/bin/Altix/agcm5.v00ml1soffimat.single.PE1 is created)
↑

Run †

  1. cd $AGCMDIR/sh_kei and Edit ipilps_test.sh (modify AGCMDIR and LOC) 
    $ cd $AGCMDIR/sh_kei
$ vi ipilps_test.sh
  setenv AGCMDIR  /home/kei/Desktop/agcm5.6
  #setenv LOC 012138 ; Tumbarumba
  #setenv LOC 148055 ; Munich
  setenv LOC 300087 ; Manaus
  2. Edit gt2grd_dble.sh if necessary. 
    $ vi gt2grd_dble.sh
  make FC=ifort FFLAGS="-O -assume byterecl -convert big_endian" ${PROG}
  3. Run ipilps_test.sh 
    $ ./ipilps_test.sh
↑

Check output in GrADS †

  • cd $AGCMDIR/out/ and validate the results.
    • ipilps_test_kei/${LOC}_grd : validation data
    • ipipls_test/${LOC}_grd : your results
    • isotopic variables : IGPRCT, ietflx, irunoff, etc.
    • non-isotopic variables : GPRCT, etflx, runoff, etc.

for example

$ cd $AGCMDIR/out
$ grads -l
> open ipilps_test/300087_grd/IGPRCT.ctl
> open ipilps_test_kei/300087_grd/IGPRCT.ctl
> open ipilps_test/300087_grd/GPRCT.ctl
> open ipilps_test_kei/300087_grd/IGPRCT.ctl
> set t 1 3000
> d var1.1/var.3
> d var1.2/var.4
  ## two lines should be perfectly matched each other.