Changeset 17
- Timestamp:
- 03/21/07 14:18:07 (17 years ago)
- Location:
- trunk
- Files:
-
- 2 deleted
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/geomag_profile.sh
r10 r17 1 # !/bin/sh 1 # !/bin/sh -xv 2 2 # 3 3 # module : … … 6 6 # usage : 7 7 # online usage 8 # $ . geomag_profile.sh 8 # $ . geomag_profile.sh -d directory 9 9 # 10 10 # in ${HOME}/.profile, add the following line 11 # . geomag_profile.sh 11 # . geomag_profile.sh -d directory 12 # 13 # examples : 14 # for fplod, on aedon.locean-ipsl.upmc.fr 15 # $ cd /usr/home/fplod/incas/geomag/geomag_ws 16 # $ . geomag_profile.sh -d /usr/home/fplod/incas/geomag/geomag_ws/ 17 # for reee522 on rhodes.idris.fr 18 # $ cd ${HOME}/incas/geomag/geomag_ws/ 19 # $ . geomag_profile.sh -d ${HOME}/incas/geomag/geomag_ws/ 12 20 # 13 21 # original location … … 18 26 # ++ machine dependant 19 27 # ++ besoin de posix 28 # ++ pas de MANPATH defini par déut sur zeus 29 # fplod 2007-03-13T16:12:41Z zeus.locean-ipsl.upmc.fr (Linux) 30 # add creation of $GEOMAG_ID and $GEOMAG_OD if they don't exist 20 31 # reee522 2006-11-23T11:41:41Z rhodes (IRIX64) 21 32 # pas de set -o posix sur rhodes (IRIX64) … … 24 35 # creation 25 36 # 37 # 38 command=geomag_profile.sh 39 # 40 usage=" Usage : ${command} -d directory" 41 # 42 while [ ! -z "${1}" ] # ++ pb bash 43 do 44 case ${1} in 45 -d) # directory choosen by user 46 directory=${2} 47 shift 48 ;; 49 *) # other choice 50 echo "${usage}" 51 # nb : no exit because this file should be launched by login process 52 ;; 53 esac 54 shift # next flag 55 done 56 57 # check for ${directory} 58 if [ ! -d ${directory} ] 59 then 60 echo " eee : ${directory} not found" 61 # nb : no exit because this file should be launched by login process 62 fi 63 # 64 # check for permission on directory 65 if [ ! -x ${directory} ] 66 then 67 echo " eee : ${directory} not reachable" 68 # nb : no exit because this file should be launched by login process 69 fi 26 70 # 27 71 set -u … … 36 80 esac 37 81 # 38 command='geomag_profile.sh'39 82 # 40 GEOMAG=${HOME}/incas/geomag/ # ++ 41 # test if $GEOMAG available 42 if [ ! -d ${GEOMAG} ] 43 then 44 echo "${command} : eee : ${GEOMAG} directory not available" 45 # nb : no exit because this file should be launched by login process 46 # ++ how to worn user 47 fi 83 GEOMAG=${directory} 84 export GEOMAG 48 85 # 49 86 # add GEOMAG tools to PATH … … 82 119 GEOMAG_ID=${HOME}/geomag_d/ 83 120 export GEOMAG_ID 121 if [ ! -d ${GEOMAG_ID} ] 122 then 123 mkdir -p ${GEOMAG_ID} 124 echo "${command} : iii : creation of \${GEOMAG_ID}" 125 fi 84 126 GEOMAG_OD=${HOME}/geomag_d/ 85 127 export GEOMAG_OD 128 if [ ! -d ${GEOMAG_ID} ] 129 then 130 mkdir -p ${GEOMAG_OD} 131 echo "${command} : iii : creation of \${GEOMAG_OD}" 132 fi 86 133 # 87 134 # end 88 135 # nb : no exit because this file should be launched by login process 89 # ++ how to worn user -
trunk/initorca2_bab.pro
r2 r17 1 ; + 2 ; @file_comments 3 ; initialisation on ORCA2 grid 1 4 ; 5 ; @history 6 ; F. Vivier 2003 7 ; F. Pinsard adaptation to new SAXO 8 ; according to SAXO ./SRC/Tests/tst_initorca2.pro 2 9 ; 3 @common 10 ; @todo 11 ; parametrisation for ORCA025 12 ; 13 ; zminmesh and izmaxmesh may be set to 0,0 because we only work with 14 ; surface problem (no z) whatever the resolution of ORCA. this will safe lot 15 ; of memory 16 ; 17 ; use ncdf_getaxis to recover meshmask dimension 18 ; on x, retrieve -1 and -1 (bande de recouvrement haut et bas) 19 ; on y, for ORCA2, retrieve -2 (bande de recouvrement ++ et ++) 20 ; 21 ; reminder : ORCA2 = 182 x 149 and ORCA025 1442 x 1021 22 ; 23 ; for ORCA025, user /keepcont in triangule call. 24 ; 25 ; @version 26 ; $Id$ 27 ;- 28 ; 29 @cm_4mesh 30 @cm_4cal 31 @cm_4data 32 ; 4 33 ;--------------------------------------------------------- 5 ; bo rnes de la grille34 ; boundaries of the read grid regarding to the original grid 6 35 ;---------------------------------------------------------- 7 ;ixminmesh =2 8 ;ixmaxmesh =181 9 ; ; 10 ; iyminmesh =0 11 ; iymaxmesh =147 12 36 ; 13 37 ixminmesh =1 14 38 ixmaxmesh =180 15 16 ;ixminmesh =017 ;ixmaxmesh =17918 39 ; 19 40 iyminmesh =0 … … 26 47 time = 0 27 48 ;------------------------------------------------------ 28 ; lecture de la grille49 ; read the grid 29 50 ;------------------------------------------------------ 30 51 iodir = '/usr/work/sur/fvi/OPA/ORCA2/' 31 52 32 ;key_shift = 3133 key_shift = 034 53 key_stride = [1, 1, 1] 35 54 ncdf_meshlec, 'meshmask_bab.nc', glamboundary = [20, 380] 36 ;, glamboundary = [20, 380]37 ;,/checkdat38 55 ;------------------------------------------------------------- 39 56 domdef 40 57 ; 41 key_periodique = 1 42 triangles = triangule() 58 triangles_list = triangule() 43 59 ;---------------------------------------------------------- 44 ; bo rnes pour la lecture des donnees60 ; boundaries of the data regarding to the original grid 45 61 ;---------------------------------------------------------- 46 62 jpidta = jpiglo … … 53 69 izmindta = 0 54 70 izmaxdta = jpkdta-1 55 71 ;---------------------------------------------------------- 72 @updateold -
trunk/rseries_ncdf.pro
r13 r17 250 250 ; if we need to read more than one file, 251 251 ; first we read the first file 252 ;++ print,var,date1+century*1000000L, datelast[0]+century*1000000L 253 ;++READ, B, PROMPT='Enter Name: ' 252 254 res1 = read_ncdf(var, date1+century*1000000L, datelast[0]+century*1000000L $ 253 255 , filename = filename, /nostruct, _extra = ex) -
trunk/step1_diff.pro
r16 r17 13 13 ; 14 14 ; @examples 15 ; test (for header and delta) 15 ; test (for header and delta) 16 16 ; to compute difference between cond_sed in cond_sed_ORCA2.nc and 17 17 ; itself and write it in ginette.nc : 18 ; rhodes\@reee522$ cd /home/rech/eee/reee522/incas/geomag/geomag_ws 19 ; rhodes\@reee522$ idl 20 ; IDL> step1_diff, 'cond_sed_ORCA2.nc', 'cond_sed_ORCA2.nc', 'ginette.nc' 18 ; $ cd ${GEOMAG} 19 ; $ idl 20 ; IDL> file1=getenv('GEOMAG') + '/cond_sed_ORCA2.nc' 21 ; IDL> file2=getenv('GEOMAG_OD') + '/cond_sed_ORCA2.nc' 22 ; IDL> step1_diff, file1, file2, 'ginette.nc' 21 23 ; values of delta must be 0 everywhere 22 ; rhodes\@reee522$ rm ginette.nc24 ; $ rm ginette.nc 23 25 ; idem with BR file 24 ; IDL> step1_diff, 'Br_ORCA2.nc', 'Br_ORCA2.nc', 'ginette.nc' 26 ; IDL> file1=getenv('GEOMAG') + '/Br_ORCA2.nc' 27 ; IDL> file2=getenv('GEOMAG_OD') + '/Br_ORCA2.nc' 28 ; IDL> step1_diff, file1, file2, 'ginette.nc' 25 29 ; 26 30 ; real life … … 28 32 ; /home/rech/eee/reee522/geomag_d/cond_sed_ORCA2.nc and write it in 29 33 ; cond_sed_ORCA2_diff.nc : 30 ; rhodes\@reee522$ cd /home/rech/eee/reee522/incas/geomag/geomag_ws31 ; rhodes\@reee522$ idl32 ; IDL> step1_diff, 'cond_sed_ORCA2.nc', $33 ; '/home/rech/eee/reee522/geomag_d/cond_sed_ORCA2.nc', $34 ; 35 ; IDL> step1_diff, 'Br_ORCA2.nc', $36 ; '/home/rech/eee/reee522/geomag_d/Br_ORCA2.nc', $37 ; 34 ; $ cd ${GEOMAG} 35 ; $ idl 36 ; IDL> file1=getenv('GEOMAG') + '/cond_sed_ORCA2.nc' 37 ; IDL> file2=getenv('GEOMAG_OD') + '/cond_sed_ORCA2.nc' 38 ; IDL> step1_diff, file1, file2, 'cond_sed_ORCA2_diff.nc' 39 ; IDL> file1=getenv('GEOMAG') + '/Br_ORCA2.nc' 40 ; IDL> file2=getenv('GEOMAG_OD') + '/Br_ORCA2.nc' 41 ; IDL> step1_diff, file1, file2, 'Br_ORCA2_diff.nc' 38 42 ; 39 43 ; to see the difference, for example … … 41 45 ; IDL> xxx,'Br_ORCA2_diff.nc' 42 46 ; select xy on plt wigdet 43 ; 47 ; 44 48 ; @history 45 49 ; reee522 2006-12-18T10:38:49Z rhodes (IRIX64) … … 53 57 ; $Id$ 54 58 ; 55 PRO step1_diff, file1, file2,file359 PRO step1_diff, file1, file2, file3 56 60 ; 57 61 ncverbose=1 … … 102 106 checknonzero=where(delta NE 0.,count) 103 107 IF count EQ 0 THEN BEGIN 104 msg = ' 108 msg = 'iii : delta is zero everywhere' 105 109 PRINT, msg 106 110 ENDIF ELSE BEGIN 107 msg = ' 111 msg = 'iii : delta is not zero ' + STRING(count) + ' times' 108 112 PRINT, msg 109 113 ENDELSE
Note: See TracChangeset
for help on using the changeset viewer.