function cmpstat

****************************************************************
*****                                                      *****
*****  This script is the driver for calling: pltdmap      *****
*****  pltdmap compares the RMS and STD from 2 experiments *****
*****  (EXP0 and EXP1 from the corcmp.rc)                  *****
*****                                                      *****
****************************************************************

'reinit'
'run getenv "GEOSUTIL"'
             geosutil = result

n = 0

'getresource corcmp.rc EXP'n ;  exp.n = result
'getresource corcmp.rc DSC'n ; desc.n = result
 if( exp.n != NULL | desc.n  != NULL )
     length = strlen(desc.n)
     desc.n = substr(desc.n,2,length)
          n = n+1
 else
     say 'You must supply a CONTROL and COMPARISON experiment list'
     say 'in the file:  corcmp.rc'
     return
 endif

while( n >= 0 )
'getresource corcmp.rc EXP'n ;  exp.n = result
'getresource corcmp.rc DSC'n ; desc.n = result
 if( exp.n != NULL | desc.n  != NULL )
     length = strlen(desc.n)
     desc.n = substr(desc.n,2,length)
          n = n+1
 else
     ntot = n 
        n = -999
 endif
endwhile

* Experiment 0
* ------------
'!remove 'exp.0'/setup3'
'!make_globl_ctl1 'exp.0

'!remove setup3_'desc.0
'!/bin/ln -s 'exp.0'/setup3 setup3_'desc.0
'exec    setup3_'desc.0

'setlons'
'sety'
'setz'
'statmak h 'desc.0
'statmak u 'desc.0
'statmak v 'desc.0
'statmak t 'desc.0
'set z 1'
'statmak p 'desc.0

'closefiles'
'!remove TINC.txt'

* Experiment 1
* ------------
'!remove 'exp.1'/setup3'
'!make_globl_ctl1 'exp.1

'!remove setup3_'desc.1
'!/bin/ln -s 'exp.1'/setup3 setup3_'desc.1
'exec    setup3_'desc.1

'setlons'
'sety'
'setz'
'statmak h 'desc.1
'statmak u 'desc.1
'statmak v 'desc.1
'statmak t 'desc.1
'set z 1'
'statmak p 'desc.1


'run 'geosutil'/plots/fcst/pltdmap.gs 'desc.1' 'desc.0