Changes between Version 7 and Version 8 of Documentation/UserGuide/LoopIGCM


Ignore:
Timestamp:
2020-02-28T14:18:21+01:00 (4 years ago)
Author:
maignan
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • Documentation/UserGuide/LoopIGCM

    v7 v8  
    55Let's say that you want to run a 200Y long simulation, but you only want to use a few years of forcing data (we won't worry about why you want to do this at the moment). This is pretty easy to do with libIGCM by slightly modifying Josefine's instructions for running with analytical spinup, which I've mostly copied here. 
    66 
    7 NOTE: You have to loop over multiples of four years.  If you do not, something strange happens in the leap years and only 11 out of the 12 months are printed in the history files sometimes.  I talked to Josefine about it and we couldn't figure out why, so we're just using this as a fix. I have no idea how this works for things like FLUXNET validation. 
    8  
    9 NOTE2: Because of the previous note, take care to make sure your leap years line up in your looping years and your "real" years. I'll explain this more below. 
    10  
    11 NOTE3: You need libIGCM tag libIGCM_v2.0 (or trunk since revision 819). 
    12  
    137== How it works to use cyclic forcing == 
    148By adding the 2 variables '''!CyclicBegin''' and '''!CyclicEnd''' to config.card the cyclic functionality can be used.  
    15 !CyclicBegin is the first year of the forcing file to be used and !CyclicEnd is the last year. So if I wanted to loop over years 1950-1953 (this is a multiple of four!), I would add the following to config.card: 
     9!CyclicBegin is the first year of the forcing file to be used and !CyclicEnd is the last year. So if I wanted to loop over years 1950-1953, I would add the following to config.card: 
    1610 
    1711{{{ 
     
    2014}}} 
    2115 
    22 libIGCM will use these variables to calculate the 2 new variables '''!CyclicYear''' and '''!CyclicYear_p1'''. !CyclicYear is the year within the forcing period [!CyclicBegin, !CyclicEnd]. !CyclicYear_p1 is the next coming cyclic year in the same interval (to be used for the PFTmap). These 2 new variables have to be used in orchidee_ol.card and sechiba.card to copy the correct forcing files. 
     16libIGCM will use these variables to calculate the new variable '''!CyclicYear'''. !CyclicYear is the year within the forcing period [!CyclicBegin, !CyclicEnd]. This new variable ha to be used in orchidee_ol.card and sechiba.card to copy the correct forcing files. 
    2317 
    2418Now let's say that I want to run for 10 total years (notice this is not a multiple of the number of cyclic years. If you are doing analytical spinup, it has to be, but not if you're just doing cyclic forcing. Again, we won't ask why you want to do this). We need to change config.card. 
     
    4337 
    4438=== COMP/sechiba.card === 
    45 The PFTmap file must be copied using the variable !CyclicYear_p1 which is the next coming forcing year. In sechiba.card, remove PFTmap from section [!InitialStateFiles]. Notice that you have to remove ",\" in the end of the line with floodplains.nc. This section becomes: 
     39The PFTmap file must be copied using the variable !CyclicYear which is the next coming forcing year. In sechiba.card, remove PFTmap from section [!InitialStateFiles]. Notice that you have to remove ",\" in the end of the line with floodplains.nc. This section becomes: 
    4640 
    4741{{{ 
     
    6963}}} 
    7064 
    71 Notice that you need to remove the additional "List=  ()" line. No specific treatment in sechiba.driver will now be done for the PFTmap.nc file. It is not necessary to change in sechiba.driver. For the first year without restart file the PFTmap the will be used corresponds to the next coming year. For all next coming years the correct PFTmap is used. It is not possible to come around this problem while running with a period length one year. 
     65Notice that you need to remove the additional "List=  ()" line. No specific treatment in sechiba.driver will now be done for the PFTmap.nc file. It is not necessary to change in sechiba.driver.  
    7266 
    7367=== PARAM/run.def ===