Changes between Version 7 and Version 8 of Documentation/UserGuide/LoopIGCM
- Timestamp:
- 2020-02-28T14:18:21+01:00 (4 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
Documentation/UserGuide/LoopIGCM
v7 v8 5 5 Let'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. 6 6 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 13 7 == How it works to use cyclic forcing == 14 8 By 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: 16 10 17 11 {{{ … … 20 14 }}} 21 15 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 haveto be used in orchidee_ol.card and sechiba.card to copy the correct forcing files.16 libIGCM 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. 23 17 24 18 Now 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. … … 43 37 44 38 === COMP/sechiba.card === 45 The PFTmap file must be copied using the variable !CyclicYear _p1which 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:39 The 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: 46 40 47 41 {{{ … … 69 63 }}} 70 64 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.65 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. 72 66 73 67 === PARAM/run.def ===