Changeset 64


Ignore:
Timestamp:
08/01/11 17:33:41 (13 years ago)
Author:
jbrlod
Message:

make kt the effective time step

Location:
trunk
Files:
9 edited

Legend:

Unmodified
Added
Removed

  • trunk/scripts/run_batch_sinobad.sh

    r55 r64  
    33#PBS -o ../scripts/script_outputs/sinobad.out 
    44#PBS -q short 
    5 #PBS -l ncpus=16,nodes=1:x64 
     5#PBS -l ncpus=24,nodes=1:x64 
    66 
    77#run reference experience 

  • trunk/scripts/sinobad.i

    r49 r64  
    2020set_modeltime 0 
    2121FORWARD 
    22 #xwriteout 50 obs_48 ../data_out/exp_T.nc 
    23 #savestate gcx_dynspg_flt 1 ijk 0 A 3 ../data_out/gcx_dynspg_flt_jb.dat 
    24 #savestate gcb_dynspg_flt 1 ijk 0 A 3 ../data_out/gcb_dynspg_flt_jb.dat 
    25 #savestate gcx2 1 ijk 0 A 3 ../data_out/gcx2_jb.dat 
    2622 
    27 #goto FINRUN 
     23 
    2824 
    2925##OBSERVATION 
    3026xwriteout 50 obs_48 ../data_out/exp_T.nc 
    31 savestate sshb  1   ij   50    A       1       ../data_out/sshb_obs_48.dat    
    32 savestate tb  1   ijk   50    A       1       ../data_out/tb_obs_48.dat    
    33 savestate sb  1   ijk   50    A       1       ../data_out/sb_obs_48.dat    
    34 savestate ub  1   ijk   50    A       1       ../data_out/ub_obs_48.dat    
    35 savestate vb  1   ijk   50    A       1       ../data_out/vb_obs_48.dat    
     27 
     28goto FINRUN 
     29 
     30savestate sshb  1   ij   50    A       1       ../data_out/sshb_obs_50at    
     31savestate tb  1   ijk   50    A       1       ../data_out/tb_obs_50.dat    
     32savestate sb  1   ijk   50    A       1       ../data_out/sb_obs_50.dat    
     33savestate ub  1   ijk   50    A       1       ../data_out/ub_obs_50.dat    
     34savestate vb  1   ijk   50    A       1       ../data_out/vb_obs_50.dat    
    3635 
    3736loadobs sshb  1   ij   50    A       1       ../data_out/sshb_obs_48.dat D    

  • trunk/src/.TERM_xterm

    r52 r64  
    1111 lun. juin 27 15:44:44 CEST 2011 ... ariane.locean-ipsl.upmc.fr  
    1212 mer. juin 29 16:45:11 CEST 2011 ... ariane.locean-ipsl.upmc.fr  
     13 lun. août  1 16:18:44 CEST 2011 ... ariane  

  • trunk/src/sbcgyre_emp1.h

    r1 r64  
    1717  else{ 
    1818    //      ! current day (in hours) since january the 1st of the current year 
    19     double  ztime = double(Yt-TU) * rdt / (rmmss * rhhmm)  //       !  total incrementation (in hours) 
     19    /*    double  ztime = double(Yt-TU) * rdt / (rmmss * rhhmm)  //       !  total incrementation (in hours) 
    2020      - (nyear  - 1) * rjjhh * raajj;  //           !  minus years since beginning of experiment (in hours) 
    21      
     21    */ 
     22    double  ztime = kt * rdt / (rmmss * rhhmm)  //       !  total incrementation (in hours) 
     23      - (nyear  - 1) * rjjhh * raajj;  //           !  minus years since beginning of experiment (in hours) 
    2224    double  ztimemax1 = ((5.*30.)+21.)* 24.;//                     ! 21th june     at 24h in hours 
    2325    double  ztimemin1 = ztimemax1 + rjjhh * raajj / 2.;//            ! 21th december        in hours 

  • trunk/src/sbcgyre_qsr_qns.h

    r1 r64  
    2424  else{ 
    2525    //      ! current day (in hours) since january the 1st of the current year 
    26     double  ztime = double(Yt-TU) * rdt / (rmmss * rhhmm)  //       !  total incrementation (in hours) 
     26    /*    double  ztime = double(Yt-TU) * rdt / (rmmss * rhhmm)  //       !  total incrementation (in hours) 
     27          - (nyear  - 1) * rjjhh * raajj;  //           !  minus years since beginning of experiment (in hours) 
     28    */ 
     29    double  ztime = kt * rdt / (rmmss * rhhmm)  //       !  total incrementation (in hours) 
    2730      - (nyear  - 1) * rjjhh * raajj;  //           !  minus years since beginning of experiment (in hours) 
    28      
    2931    double  ztimemax1 = ((5.*30.)+21.)* 24.;//                     ! 21th june     at 24h in hours 
    3032    double  ztimemin1 = ztimemax1 + rjjhh * raajj / 2.;//            ! 21th december        in hours 

  • trunk/src/sbcgyre_utau_vtau.h

    r9 r64  
    1919    //      !accumulates days of previous months of this year 
    2020    //      ! day (in hours) since january the 1st 
    21     double  ztime = float(Yt-TU)*rdt/(rmmss*rhhmm)//  ! incrementation in hour 
     21    /*    double  ztime = float(Yt-TU)*rdt/(rmmss*rhhmm)//  ! incrementation in hour 
     22          - (nyear - 1) * rjjhh * raajj;//         !  - nber of hours the precedent years */ 
     23    double  ztime =  kt*rdt/(rmmss*rhhmm)//  ! incrementation in hour 
    2224      - (nyear - 1) * rjjhh * raajj;//         !  - nber of hours the precedent years 
    2325    double  ztimemax = ((5.*30.)+21.)* 24.;//               ! 21th june     in hours 

  • trunk/src/sinobad.d

    r61 r64  
    11311131insert_fct  arg xchangesavemode 
    11321132 
     1133#ifdef FILTER 
    11331134insert_fct  arg load_eb 
     1135#endif 
    11341136 
    11351137insert_fct arg xwriteout 
    11361138insert_fct arg xwritegrad 
    11371139insert_fct arg xinitnc 
     1140insert_fct arg init_euler 
    11381141//#insert_fct  arg load_shape_func 
    11391142//#insert_fct  arg load_stdev_pca 

  • trunk/src/sinobad.h

    r61 r64  
    5050int nksr; 
    5151void xtraqsr_init();// init of solar radiation penetration 
    52  
     52double Ytold; //pour voir se pas de temps 
     53double ktinit; //premier pas ocean 
    5354#   define  rdttra( k ) rdt 
    5455double r2dt; 
    … …  
    133134  xdisplay(); 
    134135  xtraqsr_init(); 
    135    
     136  ktinit=0; //A mettre dans la namelist ???? 
    136137} 
    137138//____________________________________________________________________________ 
    … …  
    147148  fprintf(p,"%23.16e\n",YTotalCost); 
    148149  fclose(p); 
     150  Ytold=Yt; 
     151  kt=ktinit; 
    149152} 
    150153//_____________________________________________________________________________ 
    … …  
    164167void forward_before (int ctrp) 
    165168{       // permet d'intervenir si besoin avant le forward 
    166   kt=Yt-TU; //oceanic time step 
    167169  
     170 
     171  kt=ktinit + (Yt-TU); 
     172  printf("kt=%lf\n",kt); 
    168173    if(Yt==TU+1 && neuler==0) 
    169174      r2dt=rdt; 
    … …  
    174179void forward_after (int ctrp) 
    175180{                       // permet d'intervenir si besoin aprÚs le forward 
    176   printf("neuler=%d\n",neuler); 
     181  //printf("neuler=%d\n",neuler); 
    177182} 
    178183//_____________________________________________________________________________ 
    … …  
    959964   
    960965} 
     966 
     967void init_euler(int argc, char *argv[]){ 
     968  if (argc!=2) 
     969    { 
     970      printf("\nwrong number of argument in init_euler : specify 0 or 1\n"); 
     971      exit(1); 
     972    } 
     973  else 
     974    { 
     975      int n; 
     976      n=atoi(argv[1]); 
     977      if(n!=0 && n!=1) 
     978        { 
     979          printf("\nwrong argument in init_euler : specify 0 or 1 in command init_euler\n"); 
     980          exit(1); 
     981        } 
     982      neuler=n; 
     983    } 
     984} 
Note: See TracChangeset for help on using the changeset viewer.